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These  Investigations  were  carried  out  under  the  authority 
contained  In  Contract  No.  DA-44- 009  ENG  4420,  "Analytical 


Analysis,  and  Appl  ication  of 


Compensating  Equations  to  the  General  Block  Trlangulatlon 


Adjustment  .Program." 

_  The  Investigation  was  conducted  under  the  d I  rect  i  on  of 

Professor  Charles  1.  MlTler  and  Mr.  E.  Phillip  Gladding,  Instructor 


.Department  of.  CjvJi  .Engineer ing.  The  Investigation  was  conducted, 
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Z  lad  M,  Ellas,  and  Frank  S.  Greatorex,  Research  Assistants  In 
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.  , ,  --  V:r  -  i  ABSTRACT  • 


. The  objective  of  the  activities  reported  Is  to  effect  Improved 

Z^ZJIil^accuracy  Trv-the  supplied  General  Block  Tr langulatlon  digital  computer  pro- 
: ""  gram  through  Incorporation  In  the  program  means  of  error  adjustment  and  .  — r 
compensation.  The  first  volume  of  the  report  presents: 

1.  -The  nature  of  random  and  systematic  errors  and  the  basic 
techniques  for  treating  their  effects  as  applicable  to  the  analytical 
photogrammetrlc  problem; 

2.  the  basic  least  squares  method  and  Its  Incorporation  In  the 
computer  program; 

3.  complete  mathematical  description  of  the  program; 

4.  studies  of  the  nature  and  effects  of  the  Important  error 
sources:  lens  and  camera  errors,  atmospheric  refraction,  film  distortion; 

5.  the  study  of  various  techniques  for  the  solution  of  simultaneous 

equations; 

6.  operating  instructions; 

7.  the  results,  conclusions,  and  resulting  recommendations  of  test 
runs  of  the  final  computer  program. 

The  second  volume  contains  the  appendices  which  consist  of  the  complete 
flow  charts  representing  the  original  and  final  programs. 
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.  ..  -  SUMMARY 


This  report  describes  research  directed  toward  the  general 
problem  of  adjustment  and  compensation  of  errors  In  the  analytical  .  : 
aerial  photogramnetr ic  triangulation  problem."  The  primary  object Ive 
is  the  Incorporation  In  the  General  Block  Trlangjlatlon  and  Adjustment 
Program,  supplied  by  the  government,  the  necessary  routines  for 
effective  adjustment  of  random  errors.  The  report  includes  studies  of 
the  particular  error  sources  In  the  photogrammetrlc  "system. 

The  general  method  of  Incorporating  the  least  squares  adjustment 
of  random  errors  is  presented  along  with  a  detailed  mathematical 
description  of  the  final  program.  The  program  has  been  tested  using 
real  and  fictitious  data.  While  the  program  is  effective  In  adjusting 
for  random  errors,  the  systematic  error  content  of  the  available 
data  is  dominant.  Considerable  Improvement  In  the  compensation  of  lens 
and  film  errors  is  required.  A  procedure  is  recommended  for  compen¬ 
sating  for  atmospheric  refraction  effects. 

It  has  been  shown  that: 

1.  a  general  routine  can  be  used  to  compensate  for  both 
vertical  and  non-vert  leal  photography; 

2.  very  high  altitudes  present  no  additional  problems; 

3.  four  fiducial  are  not  sufficient  for  f 1 1m  d tstort Ion 
correction  -  more  control  is  required; 

4.  the  Jordan  Diagonal izatlon  method  for  solving  simultaneous 
equations  reduces  the  operating  time; 

5.  analytical  or  numerical  representations  of  lens  distortion  may 
be  used  to  implement  radial  lens  distortion  correction. 

vi  i 


1.  INTRODUCTION  _ 

A.  Protect  Objectives 

The  primary-objectives  ofthls-contract  are  to:  (l)  study 
and  evaluate  the  errors  associated  with  analytical  photogrammetry 
for  the  purpose  of,  (2)  incorporating  a  means  of  analytically 
correcting  and/or  compensating  for  those  errors  Into  the  ERDL 
General -Block  Tr langulatlon  and  Adjustment  Program. 

This  report  presents  the  results  of  the  random  and  systematic 
error  studies  and  the  resulting  programming  Incorporated  Into  the 
subject  computer  program,  with  an  evaluation  of  the  effectiveness 
of  the  resultant  program. 

B.  Sources  and  Treatment  of  Errors 
1 .  Scope 

A  typical  complaint  associated  with  technical  reports 
Is  that  the  author  and  reader  assign  different  Interpretations 
to  fundamental  terms  and  concepts.  Such  a  situation  exists  in 
the  consideration  of  errors.  The  following  paragraphs  of 
this  section  are  Intended  to  give  the  reader  a  brief  review 
of  basic  concepts  Involved  in  the  analysis  and  treatment  of 
errors. 

There  are  three  character  1st  leal ly  different  sources 
of  variance  between  the  actual  value  of  a  physical  quantity 
and  the  value  obtained  through  measurement.  A  brief  discussion 
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of  each  of  these  three  follows. 

2.  Blunders 

Blunders  are  mistakes  made  during  the  measuring  process 
or  during  subsequent  handling  of  the  data.  Blunders  are 
usually  human  mistakes.  There  is  no  general  approach  or 
mathematical  method  for  detection,  compensation,  or  elim¬ 
ination  of  blunders.  Blunders  must  be  detected  by  some 
independent  checking  method  and  e! imlnated  before  applying 
systematic  and  random  error  considerations.  Blunder-checking 
techniques  can  only  be  determined  by  analysis  of  the  Involved 
measurement  problem  and  system. 

3.  Systematic  Errors 

Systematic  errors  are  those  errors  contained  In 
measured  quantities  which  are  deterministically  caused  by 
the  conditions  of  measurement;  that  Is,  the  conditions  at 
the  time  and  place  of  the  measurement  which  deterministically 
alter  the  value  Intended  to  represent,  a  measured  quantity  and 
result  In  systematic  errors.  This  Is  a  relatively  succinct 
and  rigorous  definition.  If  a  certain  condition  affects  the 
results  of  a  measurement,  the  condition  not  being  part  of 
that  intended  to  be  measured,  then  the  component  resulting 
from  the  condition  is  an  error.  These  conditions  are 
commonly  called  error  sources.  The  relationship  between  cause, 
the  error  source  or  condition,  and  the  effect  --  the  error  -- 
must  be  known.  The  effect  must  be  a  deterministic  function 
of  the  condition.  This  implies  that  the  value  of  the  error 


may  be  computed.  This  Is  true  If  the  source  condition  can 
be  measured  and  if  the  relationship  between  the. condition 
and  its  effect -I s  known .  If,  however,  one  of  these  requ I re-  - 
ments  cannot  be  satisfied; but  the  effect  Itself  can  be 
measured,  then  the  effect  of  the  source  condition  Is  deter¬ 
mined  with  equivalent  utility. 

4.  Random  Errors 

Random  errors  are  those  errors  which  obey  the  laws -of 
chance.  A  common  assumption  Is  that  random  errors  are 
Gaussianly  distributed.  The  Important  general  Implications 
of  this  assumption  are: 

(1)  Small  errors  are  more  likely  than 
large  errors. 

(2)  Positive  and  negative  errors  are 
equally  likely. 

Unlike  systematic  errors,  random  errors  cannot  be  deter¬ 
mined  per  se.  Consideration  of  random  errors  Is  based  upon 
the  assumption  that  they  are  distributed  Gaussianly.  Rigorous 
considerations  of  random  error  propagation  and  the  adjustment 
of  related  measured  quantities  which  contain  random  errors 
are  results  of  the  Gaussian  distribution  expression,  The 
resulting  techniques  and  procedures  are  part  of  what  is  norma 
known  as  the  Method  of  Least  Squares. 

5.  .Practical  Error  Conditions 

Before  any  mathematical  error  analysis  can  be  performed 
on  measured  data,  all  blunders  must  have  been  removed.  The 


measurement  data  ts  then  assumed  to  be  blunder-free.  The 
next  step  Is  the  removal  of  systematic  errors.  Each  piece 
of  measured  date  Is  analyzed  end  corrected  for  every  appl T- • 
cable  systematic  error  source.  A  condition  for  rigorous 
appl (cation  of  the  random  error  considerations  Is  that  all 
systematic-errors  have  been  removed.  One  must,  therefore, 
be  able  to  determine  the  values  of  al 1  systematic  errors;  -- 
but  this  Is  not  always  possible.  Results  of  random  error  - 
considerations  are  not  always  what  they  should  be  because 
of  the  unavoidable .presence  of  systematic  errors.  Once  all 
the  systematic  errors  that  can  be  determined  have  been 
removed,  then  random  error  adjustment  may  be  applied,  remem¬ 
bering  that  It  Is  no  better  than  the  degree  to  which  the 
assumptions  have  been  met.  The  final  result  of  the  random 
error  adjustment  of  the  measured  quantities  Is  based,  unlike 
systematic  error  correction,  on  the  relationship  between  the 
measured  quantities.  The  adjustment  criteria,  a  result 
Of  the  Gaussian  expression,  Is  that  the  sum  of  the  squares 
of  the  weighted  differences  between  the  observed  and  adjusted 
values  for  all  measured  quantities  Is  minimized.  The  weight 
associated  with  each  measurement  Is  Inversely  proportional 
to  the  square  of  the  standard  deviation  for  that  particular 
measurement.  This,  again,  Is  a  result  of  the  Gaussian 
assumption. 
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11.  RANDOM  ERROR  ADJUSTMENT  IN  ANALYTICAL  PHOTOGRAMMETRV 


A.  Mode  of  Attack  1  ;  ' 

Accuracy  In  analytical  photogrammetry  Is  related  to  two  con¬ 
siderations:  first,  the  underlying  geometrical  representation  of  the 
probjem 5 -second,  the  accuracy  of  the  Input  data.  The  geometrical-  . 
representation  assumes  that  1  Ight  rays  are  straight  lines- and  pass 
without  deviation  through  a  geometrical  point.  Atmospheric  refraction 
and  lens  distort  Ion  Introduce  errors  In  this  assumption. 

Film  shrinkage,  Instrumental  as  well  as  biased  errors  In  measuring 
the  coordinates  of  an  Image  point  and  other  Information  such  as  ground 
control  point  coordinates  Introduce  errors  In  the  Input  data.  The 
determination  of  ground  point  coordinates  from  the  locations  of 
Its  two  images  on  two  photographs  will  be  affected  by  two  classes 
of  errors.  One  Is  present  In  the  elements  of  exterior  orientation 
and  depends  on  the  factors  discussed  above.  The  other  Is  Introduced 
through  the  discrepancy  between  the  measured  photo  Image  coordinates 
and  the  theoretically  correct  ones. 

When  they  are  systematic,  the  errors  can  be  eliminated  by 
a  correction,  but  only  to  a  certain  degree.  What  Is  left  are 
unknown  random  errors  and  residuals  from  systematic  errors.  The 
only  way  which  Is  seen  to  study  them  and  their  effects  Is  to  study 
their  probabilistic  behavior  and  consider  all  of  them  as  random 
errors.  We  will  consider,  therefore,  that  we  have  a  set  of  random 
errors  in  the  Input  data.  In  case  of  redundant  data,  the  adjustment 
will  be  carried  according  to  the  method  of  least  squares. 


B.t  The  Dodae  Method  ~~ 

;  :  I  nth  Is  method,  residuals*  are  not  cons  I dered.  For  example , 
when  the  condition  for  rays  Oj I ^  and  0^ I 2  to  Intersect  Is  expressed, 
the  positions  of  Images  lj  and  .1,  in.  their  respective  photographs 
Is  assumed  to  be  exactly  known. -More  clearly,  In  the  Iterative  process 
of  solution  only  the  exterior  orientation  of  a  bundle  of  rays  Is 
altered  at  each  cycle,  but  the  Interior  orientation  is  fixed.  As  the 
primary  sou_rces  of  errors  are  In  the  location  of  the  theoretical 
Image  points,  a  rational  adjustment  by  the  method  of  least  squares 
should  take  Into  account  the  residuals  of  the  photo  Image  coordinates. 


In  the  following,  these  residuals  and  those  of  ground  control  data 
will  be  formed  on  the  condition  that  the  sum  square  of  the  weighted 
residuals  is  minimum. 

C.  Remarks  on  the  Notation  and  the  Derivation 

A  nondeflned  notation  has  the  same  meaning  as  In  the  Dodge 
report:  "Sequential  Presentation  of  Analytical  Aerotr (angulation 
by  the  Revised  Direct  Geodetic  Restraint  Method."  In  the  following 
derivation  of  the  completed  equations,  the  same  order  of  presentation 

The  word  "residuals"  is  used  here  and  subsequently  to  indicate  the 
measurement  errors. 
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as  In  the  above  will  be  followed  and  the  same  equation  numbers  will 
be  used.  The  reader  Is  assumed  to  be  familiar  with  the  notation  and  .7 
the  derivation  of  the  equations  of  the  Dodge  report. 

t  :  -  The  ixfMftls.  - ; :  -  _ 

Thesoareoftwoklnds:  -  -  - 

1 .  residuals  In  a "photograph  2  per  Image  —  ;  ■ 

. .  dl  and  dl 

—  .  X  ...  y  _  '  _ _ 

2.  residuals  on  the  ground  --  1,  2,  or  3  per 

control  ground  point!  :  d<$>  and/or  d  -  -  - 

and/or  dh 

Thenumber  of_ residuals  _ln_ a  given  equat Ion  depends  on  the 
number  of  Images  and  the  type  of  the  ground  point  Involved.  See 
Chart  1 . 

D.  Hethod  of  -Incorporation  of  Least  Souares  Adlustment 

The  need  to  compensate  for  random  errors  Initiated  an  examination 
of  the  least  squares  method.  A  set  of  equations  Is  said  to  be  redundant 
when  the  number  of  equations  (m)  exceeds  the  number  of  unknowns  (n). 

In  theory,  there  should  be  one,  and  only  one,  explicit  solution. 
However,  because  of  measuring  errors  In  the  measured  data,  one  could 
arrange  the  redundant  set  Into  S  ■  ^m_nj i  subsets  of  n  equations 
with  n  unknowns,  the  result  being  S  Independent  solution  sets,  with 
each  solution  vector  being  different.  To  establish  a  means  of  obtaining 
a  single  solution  vector  for  the  redundant  set,  the  method  of  least 
squares  was  introduced.  The  classic  approach  Is  to  assign  to  each 
equation  a  residual  term  V  equal  to  the  difference  between  the  equation 
reel  and  observed  values;  that  Is,  V  *  E  -  AX  -  BY  -  ...  -  FZ,  and 
then  to  minimize  the  weighted  sum  square  of  these  residual  terms.  When 
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practiced,  this  then  would  give  an  adjusted  solution  vector  which 
ts  said  to  be  the  most  probable  solution  vector  for  the  redundant 
i.  ....  set.  The  In  t  roducti  on  of  we  I  ghts  Into  t  hi  s  method  Is  done  by  mu  1 1 1  -  - 

plying  the  equations  by  a  set  of  predetermined  numbers.  In  essence, 
this  can  be  interpreted  to  mean  that  certain  equations  are  stronger, 
l.e.,  have  less  measurement  error  than  other  equations  within  the 
redundant  set.  Thlscan  be  proved  to  be  equivalent  to  writing  an 

equation  a  certain  number  of  times  more  than  a  weaker  equation.  To  . 

-summarize,  this  approach,  then,  gives  greater  significance  to  those 
_  _  equations  known  to  have  less  measuring  error,  by  minimizing  the 

weighted  sum  square  of  the  equation  residual  errors.  This  classical 
approach  Is  the  simplest  case  of  least  squares  adjustment  provided 
that  an  appropriate  selection  of  weights  for  the  redundant  set  of 
equations  is  possible.  The  question  of  weight  selection  proves  to 
be  the  heart  of  the  problem.  The  word  weight  Is  associated  with  a 
measured  quantity  and  Is  a  number  Inversely  proportional  to  the  square 
of  the  standard  deviation  or  variance  of  this  measured  quantity. 

To  Illustrate  this,  consider  the  following  set  of  3  equations 
with  2  unknowns,  x  and  y: 

a,,*  +  a12*  "  bl 
a2,x  +  a22y  -  bg 

a3,x  -  a32y  -  b3 

Assume  that  the  coefficients  of  the  equations  are  known  to  be  without 
error,  but  that  the  constant  terms  bj ,  b^,  bj  are  determined  through 
direct  measurements  having  standard  deviations  of  Jj,  Jj,  and  J3  units 
respectively.  The  classical  least  squares  method  would  solve  the 
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following  set  of  equations s 


allK  *  “ijV  •  Vl  ■  Sl 


:  -  a2!x  +  ^  -  V  “2 

a3*x  +  w  Vb3 

and  the  solution  xQ,  y'  satisfies  the  following  condition: 

of  all  the  possible  pairs  of  2  numbers  x  »  y  the 
solution  1  s  that  pair  which  when  .subs  1 1  toted  jn  the 
1.-  following  expressions 

*1  w  aiiao  *  at2^o  "  bl  - 

v2  ■  a2i*,  *  Wo  '  b2 
v3  ■  a3l  V+ Vo  •  b3 


makes 


f y*  v2  V2  N 

/  JL  +  JL  +  _i_  )  a  minimum 

A  j?  4  4  I 


In  the  case  where  both  the  coefficients  and  the  constant  terms  are 
measured  quantities  or  functions  of  measured  quantities,  the  problem 
becomes  more  complicated. 


This  Is  the  case  In  a  problem  such  as  the  "General  Trlangulatlon 
and  Adjustment  Program"  where  errors  are  Introduced  into  the  coefficients 
and  Into  the  constant  terms  through  the  control  data.  Accordingly, 
a  method  which  would  use  as  criteria,  (1)  the  minimization  of  the 
measurement  errors  themselves,  and  (2)  would  weight  each  measured 
value  as  a  function  of  the  expected  magnitude  of  Its  measurement  error, 
was  desired.  This  method  was  found  in  a  report  (l)  presented  by  D.C.  Brown. 


Before  proceeding  Into  the  main  of  the  report,  It  will  be  profitable 
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to  review  briefly  this  new  least  squares  method. 

Let  mQ  ■  the  measured  (observed)  values  —  ... 

m  :  «  the  true  value  of  the  measured  Items 
v  *  the  measuring  errors 
thus,  m  ■  mQ  +  v  • 

Let  “  xQ  ■  the  approximated  parameter  value 
x  ■  the  true  parameter  value 

_ _ r  ■  the  parameter  residual  _ 

then,  x  ■  xQ  +  r 

Accordingly,  the  conditional  equations  can  be  expressed  as: 

»  f  (xQ  +  r,  mQ  +  v) 

when  the  conditional  equations  are  approximated  by  the  first 


two  terms  of  a  Taylor  expansion  where: 

£  ■  f  (x  ,  m  ) 

O  o 

rV\ 

a  « 

m 

b  - 

i)X 

the  linearized  equations  take  the  form  of 

vrbi2r2+"--,at!vi  +  ,i2v2  +  -+£r° 

Now  the  primary  approach  is  to  minimize  the  weighted  sum  of  the 


squares  of  the  measurement  errors  or 
2  A  2  , 

S  ■  WjVj  +  Wj|V2  +  .  .  . 

where  wj  * 

d  o' ■  the  unit  variance 

and  d  j4  »  the  observation  variance. 

For  convenience  with  matrix  notation,  let6°  be  the  diagonal  matrix 
whose  elements  are  the  Inverse  of  the  Wj 
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A  •  the  matrix  of  the  "a"  elements 

B  ■  the  matrix  of  the  "b"  elements 

v  ■  the  column  matrix  of  the  measurement  errors 

r  ■  the -Column  matrix  of  the  unknown  parameter 

residuals  _ .  _  _  _  . 

£  -  the  column  matrix  of  the  1 Ineartzed  conditional 
-  -—-equation -values .  .  -  -  -  -------  -  --------- — - 

Then  equation  (7)  can  be  written  as 

Br  +  Av  +  £  a  0  '  — '  --  : 

and 

S  -  V*  (d°)  V 

At  this  stage,  Lagrange  mul tipi ler 's  are  Introduced  to  aid  In  the 
obtaining  of  a  solution;  thus,  let  -  the  row  matrix  of  Lagrange 
multipliers  and  rewrite  (ll)  as: 

S  -  V*  (d  °)  V  -  2  A.  (Br  +  Av  +&  ) 

The  conditions  that  S  be  minimum  are  obtained  by  equating  to  zero 
the  partial  derivatives  of  S  with  respect  to  the  measurement  errors 
V  and  the  parameters  residuals  r. 

This  leads  to: 

(d°f  V  -  AlA  -  0 
and  8*^-0 
solving  (13)  f°r  V  gives 
V  *  <i  °A*  A.. 

and  (15)  substituted  Into  (10)  gives 
Br  +  (A  6  °At)  X  +  £  •  0 
Then,  solving  (16)  for  X  gives 

X  .  -  (A<5  °At)  (Br  +  %-) 
and  (17),  substituted  into  (14)  gives  finally 
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: :T  ,  -  __  .  -  'T. IV.  ~  ~  - 

(Bt(A  ff  V)  B)  r  +  Bl  (A  ^  V)  *-  «  0 
_  To  check  the  matrix  dimension  validity,  let  : 

p  ■  the  number  of  unknown  parameters 

EQ.  ■  the  number  of  redundant  conditional  equations 

Q,  ■  the  number  of  measurement  errors  . 

Then ,  the  Individual  matrix  dimensions  are:  _____  . 

B  -  (Ed  x  P)  ..  . .  ...”  ;  _  . 

A  -  (Ed  x  Q)  T_ 

a  °  -  (Q  x  Q) 

_  r  -  (P  x  I)  _ 

t  -  (EQ  x  1)  _  c 

Thus, 

(A  i‘  °At)  -  (EQ  x  Q)  (Q  x  Q)  (Q  x  EQ)  -  (EQ  x  EQ) 

B*  (A  ./V)"'  -  (  P  X  EQ)  (EQ  x  Q)  -  (P  x  EQ) 

Therefore, 

-1  _ 

(8*  (A  J °At)  B)  r  -  (P  x  EQ)  (EQ  x  P)  (P  x  l)  -  (  P  x  l) 

and 

B*  (A  rf  V)  k.  -  (P  x  EQ)  (EQ  x  1)  -  (P  x  0 

Thus,  the  above  method  allows  one  to  minimize  the  weighted  sum 
square  of  the  physical  errors  associated  with  the  measured  (observed) 
data  without  Increasing  the  number  of -unknowns  over  that  defined 
by  the  classic  least  squares  method.  Nor  does  It,  in  any  way,  alter 
the  EROl  selected  approach  for  solving  the  set  of  normal  equations. 

The  basic  difference  In  computer  mathematical  logic  Is  in  the 
superposition  of  the  redundant  equations  In  forming  the  normal  equations. 
With  the  classic  least  squares  method,  each  Individual  linearized 


equation  could  be  normalized  and  then  Its  contribution  superimposed 
add  it Ivel y  Into  the  norma  I  equations .  However .with  this  new  approach , 
all  linearized  equations  associated  with  a  particular  point  must  be 
normalized  as. a  set.  Once  tills  Is^Ona,  the  normal  i  zed  set1  s 


contribution  can  be  superimposed  additive  I  y  Into  the  normal  equat  Ions.,. 


-  j  • -  In  the  "General  Block  Trlangulatlon  and  Adjustment  Program",  :  7 

\-  ,  - the  coefficients  of  the  equations  derived  by  H.  F.  Dodge. represent 

"  [  the  elements  of  the  aforementioned  "B"  matrix.  Thus,  the  work  required 

i  "  %f  this  laboratory  was  to  derive  additions  to  these  equations  which 


then  would  supply  the  elements  of  the  "A"  matrix,  l.e.,  the  matrix 
of  the  coeff  I  dents  of  the  unknown  measurement  errors .  “  7 


111.  LENS  AND  CAMERA  ERRORS 


A.  Focal  Length  - 

_ The  computer  program,  by  Its  mathematical  .organization,  requires 

the  use  of  a  constant  focal  length  value.  Therefore,  variable  focal 
length  techniques  to  compensate  for.  symmetric  radial  lens  distortion 
cannot  be  emp_[oyed_ here._The_symmetr I c_  «idla2  distortion  curves 


presently  being  supplied  by  camera  cal lbrat Ions -are-plots  of:  distortion 


vs.  radial  distance.  This  requires  the  use  of  a  constant  arbitrary 
focal  length  to  construct  the  curve.  By  def In  1 1 Ion,  this  focal  length 


Is  called  Lite  calibrated  focal  length.  In  regards  to  a  small  focal 
length  change,  It  has  been  demonstrated  that  the  resulting  propagated 
error  Is  very  smal I . 


B.  Asymmetric  and  Tangential  lens  Distortions 


While  it  has  become  apparent  that  asymmetric  and  tangential 
distortions  can  be  simulated  and  be  made  to  behave  as  systematic 
errors,  their  distortion  pattern  orientations  on  the  photograph  could 
be  a  matter  of  chance  (such  as  the  prism  effect  from  removable  filters) 
In  actual  practice,  there  appear  to  be  too  many  variables  within 
the  system  to  allow  for  an  economical  analytic  treatment  of  these 
error  sources  for  each  Job.  Accordingly,  It  would  seem  that  the 
proper  means  of  avoiding  these  sources  (following  past  operational 
procedures)  lies  in  the  "acceptance  tolerances"  defined  for  the 
optics  and  camera  manufacturers. 


C.  Principal  Point 


Three  different  determinations  of  principal  points  are  of  Interest 
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They  are!  (l)  Indicated  principal  point,  (2)  autocoi 1 Imated 
principal  point,  and  (3)  point  of  symmetry. 

__  The  Indicated  principal  point  Is  defined  as  that  point 

formed  at  the  Intersection  of  two  lines  Joining  opposite  fldu- - 

clai  marks.  :  ..  ..  -- 

-  ^  The  autocoi 1 imated  principal  point  Is  _that  point  on  the  phd.to  _ 
plane  at  which  the  zeroth  ray  (during  cal  I  brat  .{on  testing)  is 

Imaged.  T  ~  '  '  ’ . . 

- --;7ZThe  point  of  symmetry  Is  that  statistically  determined  point 
from  which  the_symmetrjc_radlal  distortion  curve  ]s  orlglrled. 

Examination  of  a  typical  planlgon  lens  calibration  test  shows 
the  point  of  symmetry  to  be  +15  microns  In  x  and  -9  microns  In  y 
from  the  Indicated  principal  point,  and  the  autocoi I Imated  principal 
point  to  be  +8  microns  in  x  and  -6  microns  In  y  from  the  Indicated 
principal  point. 

While  It  Is  not  clear  that  the  discrepancies  In  principal  point 
location  can  be  neglected,  it  Is  worthwhile  to  point  out  that  it  has 
been  demonstrated  that  a  small  principal  point  displacement  (5  microns 
In  x  and  8  microns  In  y)  has  a  negligible  effect  on  the  final  results 
(ground  point  coordinates). 

Since  the  lens  calibration  curve  Is  referenced  from  the  point 
of  symmetry,  the  calculations  of  image  corrections  for  radial 
symmetric  lens  distortion  will  also  have  to  be  referenced  from  this 
point.  The  images  themselves  will,  as  is  normal,  be  measured  relative 
to  the  Indication  principal  point.  These  coordinates  should  then 
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be  translated  to  be  referenced  from  the  autocol 1 Imated  principal 
point  as  this  point  Is,  optically  speaking,  the  image  of  the  ray 
passing  through  the  perspective  center  and  perpendicular  to  the 

film  plane.  -  -  -  —  —  - _ - _ _ _ _ 

D,  Radial  Symmetric  Lens  Distortion  '  — - 

-1.  Investigation  - 

:  The  Image  coordinate  shift  due  to  this  systematic  error  can 

be  direct ly  correlated  to  the  d I s tort  Ion  curves  suppl led  by_the-_-_ 
camera  calibrator.  This  curve  is  a  plot  of  distortion  versus  the 
readlal  distance  of  the  Image  from  the  "point  of  symmetry ."The 
main  point  of  concern,  analytically  speaking,  Is  how  to  best 
represent  this  curve  Internally  within  the  computer  program. 

The  most  obvious  approach  is  to  determine  the  coefficients 
of  an  n  order  polynomial  which  will  mathematically  simulate  the 
graphic  curve. 

If  d  ■  lens  distort  Ion 

X  *  the  distance  in  the  x  direction  of  the  Image 
5  from  the  "point  of  symmetry" 

Y  *  the  distance  In  the  y  direction  of  the  Image 
s  from  the  "point  of  symmetry" 

then  r  ,  the  radial  distance  from  the  "point  of  symmetry 
is 

l/2 

r  «  (X  2  +  Y  2) 
s  s 

and  the  equation  for  the  distortion  could  be  of  the  form: 
d  ■  Ar  +  Br  +  Cr2  +  Dr^  +  Er^+ 


V  Page  1 8 

Thecomputatlon  of  the  coefficients  (A. . F)  would  be 

Input  parameters  to  the  program,  and -they  would  have  to  be  determined 

only  after  each  camera  calibration  test,.:  ..  '  _  . 

The  primary  objection  to  the  above  treatment  is  that,  for 
certain  classes  of  lenses  (such  as  for  the  aviator) ,  the  lens  curve 
Is  not-smooth  nor  simple,  and  the  use  of  a  single  polynomial 
equation  couldronly  roughly  approximate  the  curve.- 

While  the  present  lens  being  employed  In  the  photogrammetrfc 
system  (the  Planlgon)  has  a  smooth  and  simple  curve,  It  Is  felt 
thal  the  correction  method  selected  should  be  general  enough  to  treat 
the  non-smooth,  non-simple  errors  associated  with  lenser;  which  might 
be  Introduced  Into  the  system  In  the  future.  A  good  working  solution 
to  treat  this  problem  can  be  found  In  Interpolation  methods. 

This  laboratory  has  selected  Newton's  Interpolation  formula 
with  divided  difference,  making  use  of  four  known  points  on  the 
curve  (for  each  Interpolation)  and  using  a  third  degree  polynomial 
to  calculate  the  distortion  for  points  r  microns  from  the  "point 
of  symmetry." 

The  procedure  for  this  method  Is  to  first  select  and  Input 
Into  the  computer  the  radial  distance  and  related  distortion  for 
those  points  on  the  calibration  curve  which  are  necessary  to  define 
the  curve  (i.e.,  maxlmums,  mlnlmums,  points  of  inflection,  etc.). 
Then,  when  the  distortion  for  a  photograph  image  is  desired,  the 
computer  program  selects  those  four  points  (nominally  two  curve 
control  points  on  either  side  of  the  Image)  which  will  best  define 
that  segment  of  the  calibration  curve  pertinent  to  the  Image  In 
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question.  A  third  degree  curve  is  fitted  to  these  four  points  and 
the  image  distortion  interpolated  from  this  curve. 


Let  d  and  r_»  the  distortion  and  radial  distance 

-- v  ^  v  of  the  Image  In  question 


d  ,  r.  -  the  distortion  and  radial  distance  of 
the  curve  control  points  selected  by 
.  the  program 

- - - -  2'  r2  ,  . 


then  d  -  f  (r) 

-  f(r0)  +  (r  -  rQ)  f  (r^r,)  +  (r  -  r,)  f  (r^.r, , 
+  (r  -  rQ)  (r  -  r,)  (r  -  rj  f  (r^r,  .^.rj) 

where  f(rQ)  -  do;  f(r,)  -  d,;  f(r2)  -  d.;  f(r.)  - 
The  divided  difference  Interrelationship  Is! 


,(ro!\ 

f(r,)v  m  f ^ro* rl *r2* 

><Vr2><^ 

Hrl)‘C  >ff(rl’r2'r3) 

jrf!r2,r3) 

f(r}K 


f(ro’rrr2 


thus  giving  that; 


Now,  by  definition,  pigs  (+)  distortion  Is  radially  outwards 
and  minus  (-)  distortion  is  radially  Inwards  relative  to  the  “■ 
"point  of  symmetry."  Thus,,  to  find  the  correctimage  coordinates, 
one  roust  adjust  the  distorted  measured  coordinates  Inwards  for 
plus  (+)  distortion  and  outwards  for  a  minus  (-)  distortion. 


Accordingly,  let:' 


corrected  X  coordinate 
corrected  V  coord Inate 


Thus,  XT*  *--(X  -  AX) 

C  s 


If  Xs  1  s_  (+) ,  / Jchen __X  I  s  (•« ) 

If  Xs  Is  (-),  then  X£  Is  (-) 

If  ,YS  ls.(+),  then  Yf  .js  : 

if  Y  is  (-),  then  Y  Is  (-) 
s  c 


2.  Discussion 

The  use  of  the  Newton  method  gives  a  general  approach  to  this 
problem.  It  allows  one  to  select  the  curve  control  points  which 
best  define  the  curve,  irregardless  of-  the  6r  Increment  between 
the  point,  and  provides  a  means  of  operating  with  lens  curves  which 
cannot  be  defined  explicitly  by  a  polynomial  equation. 

However,  If  the  lens  curve  is  a  simple  curve  and  can  be 
represented  by  a  polynomial,  it  Is  obvious  that  the  use  of  the 
Newton  method,  Instead  of  the  polynomial  approach,  would  require 
more  computation  time  and  thus  would  be  lacking  in  systems  efficiency. 

3-  Conclusions 

1.  If  the  lens  distortion  curve  is  a  simple,  a  smooth 
curve,  it  Is  more  efficient  to  use  an  explicit  poly¬ 
nomial  correction  equation  technique. 

2.  If  the  lens  distortion  curve  Is  a  non-smooth  and/or 
non-simple,  or  If  a  general  (all  lens)  correction 
technique  is  desired,  then  the  "Newton  Interpolation" 
technique  is  a  logical  method  to  be  used. 
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IV.  ATMOS PHEft I C  REFRACT  ION  COMPENSATE ON  - 


A,  Introduction 


Because  of  the  special  variation  of  the  refractive  index 
of  the  atmosphere,  a  light  ray  travels  along  a  curved  path.  The 


central  projection  of  a  photograph  does  not,  therefore,  represent 
correctly  the  original  bundle  of  light  rays  which  traveled  to 
the  camera  perspectlve-centerfrom  the  ground  pojnts  who-e  Images 
appear  on  the  photograph. 


b.  Ihfi.,.Pxgi>Um 

The  problem  of  atmospheric  refraction  correction  Is  that  of 
determining  the  displacement  of  an  Image  point  on  a  photograph 
so  that  the  straight  line  Joining  the  perspective  center  and 
the  ground  point  passes  through  the  displaced  image  point. 


The  usual  assumption  that  the  refractive  Index  Is  constant 
on  concentric  spheres  will  be  made.  The  center  of  these  spheres 
can  be  taken  as  the  center  of  a  spherical  earth  or  as  the 
center  of  curvature  of  the  ellipsoid  representing  the  geold  at 
the  photographed  area. 


In  the  following  only  the  significant  results  will  bo  presented 
The  derivations  of  the  formulae  to  follow  and  an  otherwise  more 
detailed  treatment  can  be  found  in  the  Phase  II  Interim  Technical 
Report  of  this  contract. 

C.  Apparent  Ground  Displacement 


In  Figure  1 : 
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A 

H 

B 

A' 


is  a  ground  point 
Is  the  exposure  station 
Is  the  nadir  '  ' 


)S  the  point  of  Intersection  of  the  tangent  to  the  light 
ray  at  H  with  the  sphere  (S)  whose  refractive  Index  Is  equal 
to  that  of  A 
Is  the  center  of  S 

The  sphere  (S)  will  be  taken  as  a  datum  for  altitudes.  — 


Let 

H 

h 

n 

n 


'H 


be  the  altitude  BN  of  the  perspective  center  H; 
be  the  altitude  of  a  point  on  the  light  ray; 
be  the  refractive  index  at  altitude  h; 
be  the  refractive  Index  at  the  exposure  station; 
be  the  angle  HA»,  HB; 

AQ  be  the  angle  CA1,  CA; 

j,  ,  be  the  angle  CA1,  A'H; 

(X  m  cos  ; 

R 


radius  of  the  sphere  (5) 


V 


_H_ 

R 


no.  i 
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The  apparent  ground  d Isplacement  Is  arc  AA'  and  is  given 


by. 


AA'  =  L,  -LK 


(!) 


where 


.  Hfp*9-  In  w  ^  ,  A  .5--^'  v  C  .  L  3 x,3c« ) 

i-r  — • Ta—  (^r  ><HC2.+  t  a  z.  «(•  n  s/ 


t?-) 


3H  d  c'  ..  C  '  i  "  'kJifr^u  C/  ) 

L*.a  -o - 11“  VC*  «.  *4  H  *' 


(3  J 


The  corresponding  angular  displacement  A0  (See  Figure  l)  Is  then 

A  *v  "  L*- 

A&  -  ^ 

The  quantities  C, .  C^...,  Cj,C£,...  appearing  In  (2)  and  (3)  are 
defined  as  fol lows! 


let 


*A 


N 


«  s  I-  U 


u  Is  then  a  function  of  h.  Ck  and  C-  for  k  -  1,2,3...  are  defined  by: 


<V  c«.  w  s  ^ 


(6J 


C^«  C^CHl 


_l  r‘ 

H  *■  / 


U1,  {■>  d  l> 


(0) 


It  should  be  noted  that  equations  (1),  (2),  and  (3)  are  the 
result  of  truncating  infinite  convergent  series.  Details  on  this 
question  and  the  associated  truncation  error  are  to  be  found  ln[i). 


0.  Correction  of  Image  Coordinates 

The  photographic  coordinates  x,y,z  of  an  image  point  I  are 
referenced  to  a  rectangular  right  handed  system  of  axes  haying  . 
the  perspective  center  H  as  the  origin.  The  photograph  is  here_-- 
assumed  to  lie  between  the  perspective  center  and  the  ground.  In 
Figure  3,  line  HP  Is  perpendicular  to  the  photograph,  P  i s  the 
principal  point  and  distance  HP  Is  the  focal- length  f.  Note  that 
the" z “coordinate  of  anyTmage  point  Is  negative  and  equal  to  ~f. 

Let  Q,  be  the  Intersection  of  the  vertical  (normal  to  the 

sphere  (s))  passing  through  H  and  the  plane  of  the  photograph, 

and  let  I  be  the  Image  point  of  a  ground  point  A.  As  was  mentioned 

earlier,  the  (curved)  light  ray  HIA  Is  In  a  plane  containing  HU, 

The  straight  line  HA  Intersects  the  photographic  plane  at  point 

I  which  is  the  corrected  image  point.  See  Figures  2  and  3-  The 

correction  of  the  Image  point  I  consists  of  determining  the 

components  A x  and  ay  of  the  vector  I  T  so  that  denoting  by  x 

c  c 

and  y  the  coordinates  of  I  will  have: 

'c  c 

Xc  »  X  -  A  X 
Yc,  -  V  "Ay 

Let  npn2,n3  135  t*ie  direction  cosines  with  respect  to  the  photo¬ 
graphic  axes  x,y,z  of  the  vertical  QH  directed  upwards.  Let 
f~2  5  2 

r  ■  distance  HI  ■  lx  +y  +z  and  let 
S  c  2.C.OS  (4a  -  1?  )  StN  ^ 

S,N  -VH  -  ft*} 
ax  and  ny  are  given  by 

(n,t  •  A<x)rS 
Z.-  nsrS 


AX  «• 


w 


Approximate  formulae  giving  II  In  the  case  of  a  vertical  photo- 

c 

graph  can  be  found  in  Cl]  where  their  accuracy  is  discussed. 

I  I 

E.  Determination  of  C ^ ,  Cg,  ...  ,  Cj,  C2,  ... 

These  quantities  are  defined  by: 


(<*) 


where  k  ■  1,2,... 
nH 


and  u  *  I 


n 
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let  M  ■  n-1 
then 

M-M 


U  -  T^r  -  (h-mh)  (i-H+n2 —  - - 

where  the  subscript  H  refers  to  the  exposure  station.  We  now  make 
the  approximation 
u  »  M-Mh. 

The  corresponding  error  Is  discussed  In  f^].  M  Is  related 

to  the  specific  mass  d  of  the  atmosphere  by  the  relation 

M  „  „ 

■j  =  constant . 

Two  atmospheres  will  be  considered. 

I .  Ideal  I  zed  Atmosphere 

The  Idealized  atmosphere  will  be  such  that  the  specific 
mass  d  varies  exponentially  with  altitude  according  to  the 

equation: 

,  .  -ah 

d  -  dQe 

where 

dQ  ■  specific  mass  at  the  datum 
a  ■  suitable  constant. 

Hwlll  then  satisfy  the  relation 

u  u  -ah 

M  ■  H  e 
o 

where 

M  *  value  of  M  at  the  datum* 
o 

Integration  of  (5)  gives 


Cx(h',  *-  ^  [t- 
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-sir,  -  aU 

c  (  l  +  &H  +.•  ■■) 


1H.' 


_ >  -an  ,  .  ...  (.am*'  (  a  H I  *• 

<3^  (HV  ST  r^irN®  [  I'  *  Cu  1  ~  +  — — • 


>-] 


( X) 


Xu)  - 


^(LK^C- >-•****<  VCM 


It  may  be  Interesting  to  note  that  C.^j  and  are-related - 

by  the  relation 

C  .  Jl.  c  -  A 

uk+l  aH  k  k+l 

Integration  of  (6)  gives  C£(H)  In  terms  of  Ck(2H)  and  ck(H) . 

C,;  (II)  r.  HoC-rX^'  -  5-H„  C.,X")  (n) 

where  M^  ■  HQe  a". 

If  a  Is  determined  on  the  condition  that  d  agrees  with  the 
ARDC  Model  Atmosphere  at  sea  level  h**0  and  an  altitude  of 
10  km. ,  It  Is  found  that 
a  ■  0. 10860  km-' , 

And  at  sea  level  for  dry  air,  ARDC  Model  Atmosphere  pressure 
and  temperature,  and  for  the  Kodak  filter  25(A),  It  Is  found 
that 

M  -  276.73  x  10'6 
o 

In  the  proceeding,  the  altitude  was  measured  with  regard  to 
a  spherical  datum  on  which  the  ground  point  lies.  Thus,  In 
the  formulae  giving  Cj(H),  €2(11),  ...  ,Cj(H), 

H  Is  the  altitude  of  the  exposure  station  wl th  respect  to 
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the  around  point  and  Hq  ts  the  value  of  M  at  the  ground 
point.  Thus  when  the  exposure  station  and  ground  point  altitudes 
with  respect  to  sea  level,  say  Hg  and  HQ,  are  given,  we  would 

have:''"””"  ...  = 


H  -  H  -H 
e  o 


MQ»  276.73  x  10"6  x  e‘aH° 

_HU-  M  e‘aH  »  276.73  x  I0'6  e'aHe 

■  H  O 


2.  ARDC  Model!  Atmosphere 


In  13  ( 


The  density  profile  of  the  ARDC  Model  Atmosphere.  Is  given 
M 


.  Since 


constant ,  we  have 


M  -  M 


o  d. 


4-  Is  tabulated  In  1  Si  and  numerical  Integrations  to 
determine  C)t  ...  ,  Cj,  C£,  ...  are  possible.  However, 
a  direct  tabulation  of  these  quantities  is  not  convenient 
since  they  depend  not  only  on  H,  but  also  on  the  altitude 
of  the  ground  point  with  respect  to  sea  level.  In  the 
following,  we  will  define  quantities  as  functions  of  H 
only,  from  which  Cj,  C2,...,  Cj,  ...  can  be  determined  easily. 
Now  (5)  and  (6)  can  be  re-wrltten  respectively: 


H  *"  dU 


f 

M  *  (N-Mr 

» J  CM-MjV''  ^ 


Mh^  '  dfo  -  Mh  “ 


MV'dk 


ZM  J  Mb*"'  dh  «•  Mh*-  ~ 


0 
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If  sea  level  Is  taken  as  datum  and  H  and  H_  are  the 

.  .  o  _ e  _  _  ,  _ _ _ _ _ 

altitudes  of  the  ground  point  and  exposure  station  res¬ 
pect  1  veiy  with  respect  -  to  sea  1 evel ,  then  we  wl 1 1  have 

nk  -  (h  -h  )k 

e  o 

and  _  _____  _  ^  _ _ _  . 


Cm,,  -  i«o  n 


H, 


( u«.,  •>« ) ' 

yC 


j  He. .  r  tlr.  "" 

( II,  -  )  V  '  r.  H  H  '''  d  h  -  ?-  M  „  /  M  h  cl  In  <  H'ti  - 

-'ll 


where  the  variable  of  Integration  h  now  has  sea  level  as 
its  origin,  With  this  same  convention,  let 

.  it 


I.  (10 


-  / 


M  h  d  li 


/ 

*  0 


I'^iO  ■-  /  hM/-1  cU, 


where  values  of  Fj,  Fj,»  ^3  ant*  F j  f°r  H  4  70  km  can  be 
found  In  Table  1 . 

C.  and  C'  are  then  given  by 


CHt-H,rCu  =  -  H, 


d. 


u. 


(.«»-  iokc;  -  f;  cm0 


-Wj-*- 


•>1 


+ 


n,f 

^  cT  k 


X  ' 


(2.0  I 

Lii) 
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where  M  and  d  are  the  values  of  H  and  d  at  any  given 
point,  and  de  Is  the  value  of  d  at,  the  exposure  station.  “ 

Is  given  In  Table  2withdQ  referring  to  sea  level .  Then 
Mq  -  276.73  x  10  0 

for  the  Kodak  filter  25(A),  dry  air,  and  ARDC  pressure  and 

temperature  at  -sea  level .  If  a  different  value  of  ^r-  say-- 

M1.  but  still  the  same  altitude  should  be  used,  Fk  shguld 
0  M!  "  J%\ 

be  corrected  by  the  factor  and  F^  by  the  factor!  -jpl 


F.  Application  of  Equations 

The  following  data  should  be  known  or  estimated  In  order 
to  apply  the  atmospheric  refraction  correction. 

1.  Altitude  of  the  exposure  station  with  regard  to  sea  level  He> 

2.  Altitude  of  the  ground  point  with  regard  to  sea  level  Hq. 

3.  The  direction  cosines  n^,  n^,  n^  of  the  vertical  passing 
through  the  exposure  station  with  regard  to  the  photographic 
axes. 

4.  Photographic  coordinates  of  the  image  point  x,y,z  ■  -f. 

5.  Local  radius  of  the  earth  R0« 

Knowing  Hq  and  He,  the  constants  C(,  . ..,  Cj,  ...  are 

determined  for  either  the  idealized  or  the  ARDC  atmosphere  using 
equations  (15)  through  (19)  for  the  former  and  Table  1  with  equations 


(20)  and  (21)  for  the  latter. 


Cos  V, 


N  <  n  1‘  * 


In 


VH  Is  determined  by 


Page  34 


The  datum  sphere  passes  through  the  ground  point.  Its-  radius  Is 

-----  — —  ■  ■  --  ij  - H  -h  -  -V.  '  ■ ,  -  - r  . 

then  R0+h0  and  XH  *  T  "  F7T  *  '  - 

O  O  . .  -  ;r_ :: 

Q0  Is  determined  by  “  ~ 

0-'.  ^h.K>*mV.h  ■  S'N  .  ■'  _ 

then  wl th 


H  -  H  -H 
e  o 

c<  ---.-C  OS.  <|t  '  ~  %- 

Formulae  (2)  and  (3)  are  used  to  compute 


A  e  * 


I  ”  I-  2. 

“X*' 


~  3-. 

■  n. 


s  Is  then  computed  using  formula  (7),  use  of  (8)  6 (9)  produces 
the  corrections  A*  and  Ay,  and  finally 
xc  “  x-  Ax 

y c  ■  y-  ay 

Table  I  contains  F|,  F^,  F^  and  Fj  only.  For  all  practical 
purposes,  the  terms  In  C^,...  and  In  C£,  C^,...  can  be  deleted. 
An  upper  bound  of  the  corresponding  truncation  error  Is  given 
In  Cl],  p.  77- 

Beyond  a  certain  altitude  H^,  the  refractive  index  becomes 
close  enough  to  unity  so  that  refraction  can  he  neglected  In  the 
region  above  H^,  In  Figure  5  the  light  ray  Is  HH^A,  H^H  can  be 
considered  as  a  straight  line,  it  Is  seen  that  for  H  x  H^,  the 
displacement  AA1  on  the  ground  and  the  corresponding  angular 
displacement  ee  are  Independent  of  the  altitude  H.  Thus  AA1  ■  L ^ 
or  46  =  —  - L,k-  are  calculated  using  equations  (2)  and  (3)  as  If 

Is 
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the  exposure  station  were  at  H^.  The  angle  1$  obtained  from 
the  angle  by  the  relation  ...  .  — 


(R+H, } sin  VH  -  (R+H)sln  Vu. 


H~ 


But  when  using  (?)- to  compute  s,  shouid.be  used. 

G .  .  .  Incorporation  of  a  Compensa t i on  Routine  1 n  the  Obi ect  Program 
Some  of  the  Information  needed  to  make  the  corrections,  such 


as.  the__aj _tltud.es jof_the  camera  staUon  and  of  the  ground  point 
“and  the  orientation  of  the  vertical  with  regard  to  the  photographic 
axes  can  only  be  estimated  at  the  beginning  of  a  Job.  The 
solution  to  this  dl ff icul ty . Is  to  include  the  routine  for  atmospheric 
refractions  In  the  Iteration  cycles  of  the  program  whereby  the 
altitudes  and  orientations  obtained  from  the  resection  at  the 
end  of  a  cycle  are  used  to  compute  the  corrections  to  every  Image 
point  coordinates  before  these  are  processed  In  the  next  cycle. 
However,  the  Inclusion  of  the  Iteration  process  should  not  be 
executed  before  a  sufficient  convergence  has  been  achieved. 

In  case  of  near  vertical  photography  and  a  reasonable  estimate 
of  the  altitudes,  the  correction  can  be  made  only  once,  before 
the  Image  point  coordinates  are  processed.  In  that  case,  a  simpler 
formula  can  be  used  to  compute  the  Image  point,  displacement  provided 
some  conditions  on  the  magnitude  of  H  and  are  met.  See  [ll  . 

It  should  not  be  overlooked  that  the  atmospheric  refraction 
corrections  should  be  made  preferably  along  with  other  systematic 
error  corrections,  such  as  lens  distortion  and  film  shrinkage,  so 
that  the  remaining  errors  can  be  treated  as  random  errors  and  be 
adjusted  as  such  through  the  least  squares  method. 
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V,  FILM  DISTORTION 

Invest Iqat ion  _  -  • 

I .  Theoretical  Considerations 

The  correction  of  aerial  photographs  for  film  distortion  Is,' 

In  general,  one  of  Interpolation.  A  certain  number  of  points  In 
the  photographic  plate  are  used  as  control  or  known  points,  arid 
-from  the  Information,  tha.t.  they,  provide  the  Images,  of  other  points  _ 
will  be  displaced  to  locate  them  as  close  as  possible  to  the  position 
that  they  would  occupy  had  there  not  been  any  distortion.  As  In  any 
Interpolation  procedure,  success  depends  on  the  number  of  known  pieces 
of  Information  available  as  well  as  on  the  knowledge  of  the  general 
pattern  that  the  system  follows.  For  example,  If  measurements  are 
made  of  the  x  and  y  coordinates  of  points  on  a  parabola,  we  can,  from 
the  known  properties  of  such  a  curve,  make  the  following  analysis: 

The  general  analytical  expression  for  a  parabola  with  axis  of 
symmetry  parallel  to  the  y  axis  Is: 

Y  "  C,  +  V  +  C3xZ  (1) 

in  which  the  three  C's  are  arbitrary  constants.  If  In  the  above 
described  manner  the  coordinates  of  three  different  points  are  measured 
say  points  PjUj.Yj),  P2^xj*y2^’  the  following  relations  can 

written  by  using  (l) . 


Y|  *  C1  +  Vl  +  C3X12 

y2  =  C|  +  C2*2  +  C3*/ 


(2) 


This  provides  a  set  of  three  equations  and  three  unknowns  which  usual ly 
have  a  unique  solution  for  the  C's.  If  the  measurements  were  carried  I n 

such  a  fashion  that  no  errors  are  made,  that  Is,  they  are  exact  (a  _ 

physical  impossibility),  the  three  C's,  when  used  In  (l)  will  give  the 
exact  formula  of  the  parabola.  However,  when  the  measurements  are  not 
exact,  the  C 's_wl  )1  only  give  an  approximate  expression  for  the  curve. 

If  more  than  three  points  on. the  parabola  aj^e  measured,  say  n,  a  set 
of  n  equations  In  three  unknowns  Is  obtained  In  which  any  three  of  them 
will  be  solvable  for  the  C's.  There  are  ni/  3'(n-3)!  different 
possjble  groups  of  three  so  that  the  same  number  of  parabolas  can  be 
formed  which  will,  In  general,  be  different  from  each  other.  To  determine 
from  them  what  parabola  to  take,  statistical  methods  should  be  used.  The 
most  suitable  approach  Is  that  of  minimizing  the  sum  of  the  squares  of 
the  res  1  dual s J 

If,  however,  we  did  not  know  that  the  curve  Is  a  parabola,  the 
problem  could  become  much  more  complicated  because  either  a  guess  of  .the 
type  of  curve  will  have  to  be  made  or  else  enough  points  will  have  to  be 
measured  In  order  to  fit  either  a  Fourier  or  an  algebraic  polynomial. 

When  this  Is  the  case,  a  study  of  the  physical  characteristics  of  the  prob 
lem  can,  In  many  Instances,  be  very  helpful  in  determining  the  exact 
function  or  at  least  an  approximation  to  the  exact  function. 


I.  See  F,  8.  Hildebrand,  Introduction  to  Numerical  Analysis.  McGraw- 
Hill,  1 956 ,  for  least  Squares  Methods. 
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The  above  example  Is  one  of  one- dimensional  Interpolation;  that  is,  the 
unknown  functl on  depends  on  only  one  varlabl e .  F i I m  d I s tor 1 1 on  cor  recti  on  I s , 
however,  a  two-dimensional  interpolation  problem  since  the  distortion 

(or  correction)  at  a  point  can,  In  general,  depend  on  two  Independent  _ 

variables  (usually  x  and  y).  Two-dimensional  Interpolation  Is  basically 
analogous  to  the  one-dlmenslonal- case- except  that- lengthier  computations  are  - 
-Involved,  and- 1 n.-genej'a I  more  data  Is  required  for  a  satisfactory  solution.  - 
One  further  complication  Is  present  In  the  film  distortion  correction 
problem.  At  any  one  point,  the  distortion  Is  not  a  scalar  but  rather  a 
vector;  that  Is, It  has  a  dl rect Ion  as  well  as  a  magnitude  associated  with 
It.  In  the  analysis  that  follows,  it  will  be  shown  that  the  vector  nature 
of  the  distortion  does  not  complicate  the  problem  appreciably. 

As  stated  earlier,  the  problem  of  compensating  for  film  distortion  can 
be  considered  as  one  of  defining  a  vector  field.  The  vector  field  provides 
the  size  and  the  direction  of  the  correction  that  should  be  applied  to 
every  point  on  the  photograph.  The  same  Image  results  as  existed  when  the 
film  was  exposed . 

The  first  problem  Is  to  determine  suitable  coordinate  axes.  Obviously, 
whatever  form  the  vector  field  has,  It  will  be  expressible  in  any  coordinate 
system.  The  selection  of  a  coordinate  system  Is  fairly  arbitrary,  and  the  system 
more  suited  to  each  phase  of  study  should  be  used.  Once  the  axes  have  been 

selected,  every  point  P  (l  1,2, . n)  In  a  photograph  can  be  located 

by  the  vector  V^k  which,  In  turn,  can  be  expressed  In  terms  of  Its 
and  components: 


kx 


+  V 


ky 


(3) 


Figure  I 

where  1  and  J  are  unit  vectors  in  the  x  and  y  direction  respectively. 

The  absolute  location  of  In  the  X,  Y  plane  Is  of  no  Importance; 
the  only  thing  of  Importance  Is  the  relative  location  of  all  the  points. 
If  the  convention  Is  adopted  that  P^  he  the  points  on  the  photograph 


and  P^  the  corresponding  points  on  the  master  grid  that  was  used  to 
make  the  prints,  we  can  see  that  If  by  definition, 


E.  -  v.  1  -  7. 

k  k  k 


<*■) 


a  transformation  of  the  set  of  points  P^  of  the  form: 

V  ■  \  ♦  rk  <s> 

will  obviously  define  a  set  of  points  with  the  exact  same  geometrical 
relations  as  those  of  P^.  The  set  of  vectors  E^  can  thus  be  said  to 
represent  the  discrete  value  of  a  vector  field  E  ■  F  (X,Y)  at  the  points 
(Xh,Yk).  This  vector  field  satisfies  the  necessary  conditions  at  all 
points  whose  coordinates  are  known  both  on  the  film  and  on  the  grid.  It 
should  be  noted,  however,  that  an  Infinite  number  of  vector  fields  take 


--  - 
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take  values  at  P^( ,  and  unless  an  infinite  number  of  points  are  known,  the 
true  E  cannot  be  singled  out,  unless  other  Information  is  given.  We  do  know, 


however.  that  F  (X,Y)  mgst  be  contlnuous  un jess  rips  or  tears  are  present 
In  the  film.  - L — 


In  an  empirical,  study  I  Ike  the  one  at  hand,  E  will  hopeful  ly  be 
determined  by  making  sufficient  observations  at  close  enough  Intervals  so. 


that  an  Interpolation-will  give  it  within  the  accuracy  desired.  As  pointed  out 
earlier,  the  interpolation  problem  Is  not  an  easy  one  because  It  Is  a  two- 
dimensional  problem  and  o iso  because  no  information  Is  yet  available  as  to 
the  shape  of  the  distortion  field. 

The  representation  of  E  requires  a  four-dimensional  space.  There  are 
two  independent  variables,  X  and  Y,  and  two  dependent  ones,  le  (magnitude 
of  E)  and  (angle  of  E  with  the  x-axis).  It  Is  found  convenient  to  break 
the  problem  down  into  two  parts,  each  dealing  with  a  three-dimensional  space 
in  the  following  manner: 
from  equation  (4) 


letting: 


\ p  Vk  -  vk  -  \  3>  -  <V +  YkJ) 

Ek  ■  (xk'  *  v  +  (Yk' "  vj 


exk  “  Xk  ■  xk 

e  ,  *  Y.  '  -  Y. 
yk  k  k 


(6) 


(7) 


(8) 


equation  (*')  can  be  rewritten: 

Ek  exk*  +  eyk^ 

In  general,  E  Is  a  function  of  X  and  Y,  The  X  and  Y  components  of  E  will 
thus  be  functions  of  X  and  Y: 
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J  -  F(X,Y)  f  (X , Y )  I  +  g(X,Y)j  •  <9) 

If  it  Is  desired  that  E  agrees  with  Vk  -  at  the  point  P^,  the  following 

relation  must  hold  _ 

Ek  -  F(Xk,Yk)  -  f(Xk.Vk)T  +  g(Xk,Yk)j  -  ej  +  eykj 

from  which: 


"exk 


do) 


--So,  once  f  and  g  are  known,  ex  and  ey  are  obtained  Identically. 

The  problem  has  now  reduced  to  the  determination  of  two  scalar 

point  functions  f(X,Y)  and  g(X,Y).  Each  of  these  functions  can  be  visualized 

as  a  surface  In  three-dimensional  space. 

There  being  no  Information  available  as  to  the  form  of  f  and  g,  a  form 
should  be  found  that  conforms  best  with  the  true  deformation.  This  can  only 

be  done  after  an  extensive  empirical  study. 

The  most  popular  Interpolation  formulas  deal  either  with  approximations 
by  means  of  polynomials  or  by  the  use  of  trigonometric  series. 

Expressing  f  and  g  as  algebraic  polynomials,  the  following  equations 

are  obtained: 


01) 


«<M>  *  cyo  *  V  *  V  '  cy3XY  +  -•  *  *  *v(M) 


where  the  C ' s  are  constants  which  will  usually  be  determined  from  the 
aval lable  data. 

At  the  points  Pk  where  the  ER  vector  Is  known,  the  following  relations 
can  be  written  by  combining  equations  (7),  (10),  and  (ll): 


*k  -  V  C„0  *  Cx.\  *  Vk  * 
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yo 


CylXk 


+  C  ,Y. 
y2  k 


-  (12) 


r  In  equa  t  Ion  (12)  the  only  unknowns  a  re  the  C 1  s.  Each  known  po  1  nt. 
gives  two  equations,  one  for  f  and  one  for  g.  The  maximum  size  of  the) 
In  ter  point  I  rig  pol  ynom  I  a  1  ,  as  well  as  the  maximum  number  of  cons  tan  ti> 


In  It,  depends  then  on  the  number  of  available  known  points. 


In  photogrammet r  Ic  practice,  only  fpu_r_  p_o!nts_.of  known  pos  1 1 1  on  fli  t- 
available  In.opch  photograph.  This  limits  the  number  of  constants  thin 
tan  bo  solved  for  to  four. 

tioiiKJ  special  cases  of-tlui  polynomial  Interpolai  Ion  are: 
a)  Uni  arm  I  Infill'  cot met Ion. 

l  oi  this  rase,  equal  lor.  (lO  .edutes  l«*i 

rCK.v)  -  tiM0  •  i  ,x 

il(X.Y)  -  t;yfl  *  f,Y  (13) 

It)  linear  nonet  Hon  lit  ai«l  Y 

This  Is  l he  tot  reel  Ion  usually  *|!tM!«il  In  phutnyrammeti  l« 

wru'ki  The  procedure  Is  It*  multiply  all  tl  1st  antes  In  the  x- 

•I, 

direction  hy  the  mtlu  whrte  df  h  the  distance  between  tin 

tl  ... 

related  fiducial  marks  tar  the  phot*Htr*ph  and  d^  Is  the  same 

distance  as  measured  In  the  master  grid.  In  the  notation  of  this 

report,  this  ease  become* • 

MM)  *  1  e„« 

«(«.*)-  cy6  *  ey,*  1,1,1 

c)  Linear  Interpolation 

This  differs  frn«  the  pftvleu*  cate  In  that  f  and  g 
are  'Inear  function*  ©f  both  x  and  V-  Thul, 

MM)  6  C*o  *  CM*  *  C,iY 
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9<XKY)  “  Cy„  *  V  +  V  (15) 

d)  Tewinkels's  Formula 

Hr.  Tewinkel  of  U.S.C.G.S.  suggests  a  correction  formula 
of  the  form: 

f(x.y)  »  c  ■  cx|x  .  cx2y  t  c 

x3 

g  (x ,  Y )  r-  c  *  C  ,X  ,  C  ,Y  >  c  ,XY  ( !  &) 

a  ’  yo  yl  y2  y3 

The  above  formulas  will  be  discussed  and  evaluated  later  in  this  report. 

Ms  a  concluding  remark  in  this  section,  it  should  be  noted  that, 

although  algebraic  polynomial  interpolation  is  the  more  popular  type  of 

correction,  there  is  yet  no  evidence  that  it  is  the  best  one  for  our 

purpose.  Furthermore,  the  above-mentioned  polynomials  are  not  the  only 

2  2  2  2 

possible  ones;  for  example  (X  +  Y  )  or  X  and  Y  can  be  used  instead  of 
XY  in  the  lost  formula.  A  trigonometric  or,  in  general,  a  tr anscendenta 1 
equation  might  give  better  results  than  an  algebraic  one.  The  decision  as 
to  which  would  be  the  more  efficient  type  of  equation  can  only  be  made  by 
a  statistical  evaluation  of  a  large  group  of  test  data. 

2.  Tests  on  the  Different  Interpolation  Methods 

Since  this  laboratory  has  no  facilities  for  the  preparation  of  data  for 
a  i i Im  distortion  study,  we  hod  to  rely  on  data  from  outside  sources. 

The  only  data  that  we  could  obtain  which  seemed  to  be  of  any  value 
for  the  project,  was  o  set  of  four  photographs  prepared  by  the  U.S.  Coast 
and  Geodetic  Survey  m  is  data  consisted  of  the  comparator  coordinates  of 
thirteen  grid  intersections  from  a  master  plate  and  the  contact  prints 
made  from  it  (on  aerograph ic  film).  The  time  of  measurements  ranged  between 
(  tve  days  to  forty-seven  days  uftei  development. 
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This  data  was  used  to  make  a  statistical  evaluation  of  the  inter-:! 
^elation  methods  discussed  In  the  previous  section..!  '/ 

A  set  of  programs  were  run  on  the  IBM  1620,  In  these  programs,  all 
four  methods  were  used  for  the  Interpolation,  and  the  constants  in  each 
formula  were  obtained  once  by  using  the  corner  flduclals  as  controls  and  once 
usjng  the  side  flducljjJs  as  controls.  After  obtaining  the  appropriate^ 
-—constants,  the  film  coordinates  were  used  In  the  formulas  and  the  corrected- 
Image  coordinates  were  compared  to  the  plate  coordinates.  The  sum  of  the 
squares  of  the  remaining  errors  In  the  X  and  Y  of  the  corrected  Image 
coordinates  gave  an  indication  of  the  effectiveness  of  the  different 
methods.  It  should  be  noted  here  that  theoretically,  in  the  first  three 
methods,  the  absolute  locations  of  the  origin  or  the  orientation  of  the 
coordinate  axes  do  not  affect  the  corrections.  In  Mr.  Tewlnkel's  method, 
however,  although  the  absolute  location  Is  not  significant,  the  orien¬ 
tation  of  the  coordinate  system  does  affect  the  corrections.  This  is 
evident  when  one  considers  that,  In  this  method,  the  curves  of  constant 
correction  are  hyperbolas  with  asymptotic  directions  parallel  to  the  X 
and  V  axes;  so  as  these  curves  of  constant  correction  depend  on  the 
orientation  of  the  axes,  the  corrections  also  depend  on  the  orientation  of 
the  oxe6.  It  should  also  be  noted  that  there  are  polynomials  whose 
corrections  would  depend  on  the  absolute  location  of  the  origin  and  the 
orientation  of  the  coordinate  system.  A  simple  example  of  this  case  is; 
f(X,Y)  -  C^X 
g(X, Y)  -  cyY 

Keeping  the  above  considerations  In  mind,  the  following  tests  were  made: 
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(0  methods  (a) 1 ,  (b) ,  and  (c)  once  with  the  corners  as  controls  and  once 
w  1th  the  sides  a  s"  controls;  (2)  method  (d)  with  brig  Ins  outs  I'de  the  photo- 
graph,  once  using  the  corners  as  controls  and  once  using  the  sides  as 
controls;  (3)  method  (d)  with  origin  at  the  center  of  the  photograph  and 
axes  almost  parallel  to  the  sides  of  the  photograph,  once  using  the 
corners  as  controls  and  once  using  the  sides  as  controls;  (4)  method 
(d)-wlth  origins  at  the  center  and  the  axes  almost  along  the  diagonals  of 
the  photograph,  once  using  the  corners  as  controls  and  once  using  the  sides  as 
controls. 

In  methods  (a),  (b) ,  and  (c) ,  as  there  were  more  equations  than  unknowns, 
normal  equations  were  formed  and  the  res  ltlng  set  of  equations  were  solved 
by  the  Jordan  diagonal Izatlon  method.  The  significance  of  using  two  different 
locations  for  the  origins  for  method  (d)  will  be  discussed  In  the  next 
section.  In  the  following  table,  the  results  of  these  runs  on  three  sets 
of  data  are  listed,  using  the  sum  of  the  squares  of  the  errors  In  X  and  Y 
of  the  corrected  Images  as  a  measure  of  the  efficiency  of  each  method. 


(a),  (b) .  and  (c)  refer  to  the  methods  described  In  the  last  section. 
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------  .  - 7  ‘  '  ... 

...  --1--. 

.  '  A.vi.  .  - --- 

.  -  - 

SUM  SQUARE 

OF  X  AHft  V  ERRORS  FOR 

•. '  . — .-  r-  — -  --  ------  — 

THIRTEEN  POINTS  IN  (MICRONS)1,  ,  1 

.  .  .  '  -----  .  - t " 

METHOO  of  interpolation 

DATA  H  I 

KB! 

OAT  A  *9 

.....  . 

CA)  with  corners  as  controls 

73  AA 

e  EAR 

TR3R 

cai  n/ith  Sioes  as  controls 

T&so 

8  5  at 

TOTS 

(.b)  WITH  CORNERS  as  CONTROLS 

5ASP 

THIS 

(.b)  v/lf  HT«rOE9~AS-  controls 

K213»9I 

“(C)  WITH  CORNERS  AS  CONTROLS 

SOTS 

t  (»iStT  : 

— Sots 

CCI  WITH  SIDES  AR  CONTROLS 

SHOE 

s  top 

SOSA-  “ 

(d)  WITH  ORIGIN  OUTSIDE  ANO  AXES 
PARALLEL  TO  THE  SIDES  OF  THE 
PHOTOGRAPH  ANO  USING  THE 

CORNERS  AS  CONTROLS, 

1351 

qis 

736 

<dl  WITH  ORIGIN  OUTSlOE  ANO  AXES 

parallel  to  the  sides  of  the 
photograph  USING  THE  SiOEA 

AS  CONTROLS. 

io< 

a.  i  » 107 

3,6  m  .o’ 

id)  WITH  ORIGIN  AT  CENTER  AND 

AXES  PARALLEL  TO  THE  SIDES 

ANO  CORNERS  AS  CONTROLS, 

1  as  1 

qis 

736 

Id)  WITH  ORIGIN  AT  CENTER  ANO 

AXES  PARALLEL  TO  THE  SI013 

ANO  SIOES  AS  CONTROLS 

«• 

A  H  (O 

E,  U  *  IO* 

%  A  «  Id'* 

Idl  WITH  ORIGIN  AT  CENTER  AND 
AXES  ACROSS  THE  DIAGONALS 

OF  THE  PHOTOGRAPH  AND  THE 
CORNERS  AS  CONTROLS 

»• 

1.3  *  to 

1.9  *  IO* 

3  »  lO* 

<d)  WITH  ORIGIN  At  CENTER  ANO 
AXES  ACROSS  THE  0 1 AOONALS 

OF  THE  PHOTO  ORAPH  AND  THE 
SIDES  AS  CONTROLS 

-AP07 

7.9  S9 

RS  36 

TABLE  SHOWING  THE  SQUARE  OF  fiRRoRS 

IN  X  AMD  V  OF  THE  CORRECTED  COORDINATES  OF  THIRTEEN 
POINTS  (.TWENTY- SIX  COORDINATED)  ,  FOR  D  IPFC  R.ENT  WgTHOOS 

OF  INTERPOLATION. 

CO),  Cb),  Cc),  ANO  <4>,  REFER  TO  PAETHOOS  0E%CRlSt6 
IN  THE  LAST  SECTION. 
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Discuss  1  on 

Although  three  sets  of  data  are  not  sufficient  for  a  reliable 
statistical  evaluation  of  the  different  methods  of  interpolation,  one  can 
still  expect  :.that-the  distortion  pattern  of  these -three-photographs  after  - 
correction  would  be  Indicative  of  the  efficiency  of  each  method.  Further¬ 
more,-  one  can  expect- a  Certain  method  to  be  preferable  to  another  by 
considering  some  theoretical  aspects.  1  _  ....  _ _  : . .  . 

Methods  (a)  ahd^(b)T "are  rather- Ineff  tc lent ,  and“the  sum  square-of  the 
errors  In  our  three  photographs  confirm  this  fact.  These  two  methods  gave  the 
largest  errors  In  every  case,  whether  using  the  corner  fiduclals  or  the 
side  fiducial s  as  controls.  Method  (c)  gave  better  results  than  methods 
(a)  and  (b).  Although  In  method  (c)  the  corner  fiduclals  gave  better  results 
than  the  side  fiduclals  In  all  three  cases,  the  difference  was  not  sig¬ 
nificant  enough  to  Indicate  that  corner  fiduclals  are  better  than  side 
fiduclals  for  this  method. 

Method  (d)  provided  more  problems  than  the  first  three  methods.  Although 
It  was  deduced  that  the  absolute  location  of  the  origin  should  not  affect 
the  results  of  the  Interpolation,  when  the  axes  were  parallel  to  the 
sides  of  the  photograph  and  the  side  fiduclals  were  used  as  controls,  the 
results  did  depend  on  the  location  of  the  origin.  It  It  should  be  noted 
that,  In  these  cases,  there  was  an  enormous  amount  of  errors  left  after 
correction.  Although  these  large  errors  were  not  all  due  to  round-off 
error,  there  was  evidence  that  there  was  a  large  amount  of  round-off  error 
In  the  case  where  the  origin  was  outside  the  photograph  with  the  axes 
parallel  to  the  sides  of  the  photograph  and  the  side  fiduclals  used  as 
controls.  This  wa9  evident  from  the  fact  that,  when  the  corrections  of  the 
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control  points  themselves  were  obtained,  they  had  errors  up  to  ten  microns. 
This  suggests  that  the  difference  of  the  results  due  to  change  of  origin 
was  due  to  a  different  amount  of  round-off  error  In  each  case. 

The  large  error  associated  with  the  above  case  was  due  to  using 
"weak"  points  as  controls  and  the  greatest  amounts  of  error  were  associated 
with  points  whose  XY  products  were  large.  We  noted  earlier  that,  when  using 
M  .  Tewlnkel's  formula,  the  curves  of  constant  correction  are  hyperbolas. 
Considering  the  X  correction,  C^o,  and  define  a  family  of 

hyperbolas  with  the  same  center  and  asymptotes.  Each  hyperbola  Is  defined 
as  a  function  of  f  or  the  x  correction,  and  the  magnitude  of  f  Increases 
as  the  hyperbolas  become  farther  from  the  center.  (Refer  to  Figure  2.) 

When  the  axes  were  rotated  forty-live  degrees  to  make  them  almost  along 
the  diagonals  of  the  photograph;  as  should  be  expected,  the  side  fiducial 
gave  good  results,  whereas  the  corner  flduclals  gave  ery  large,  errors  for 
points  with  large  XY  (origin  at  center  of  photograph).  This  Is,  of  course, 
exactly  the  reverse  of  the  case  where  the  n.es  are  parallel  to  the  sides. 
Note:  When  solving  for  the  constants  of  Tewlnkel's  formula,  It  Is 
theoretically  possible  that  the  determinant  of  the  matrix  Involved  will  be 
zero  and  no  solution  will  exist.  Such  a  case  will  arise  when  there  Is 
absolutely  no  error  in  the  coordinates  of  the  flduclals  on  the  film  with 
respect  to  the  fiduclals  on  the  grid  and  the  lines  joining  any  two  fiducials 
on  the  film  are  parallel  to  lines  joining  the  same  flduclals  on  the  grid. 
Although  such  a  case  is  theoretically  possible,  It  Is  highly  improbable 
that  it  will  occur,  as  there  will  always  be  some  kind  of  error  present  In 
the  coordinates  of  the  flduclals. 
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Figure  2:  This  diagram  shows  a  typical  distortion  correction  pattern 

when  Tewinkel's  method  is  used  with  side  flducials  Fj,  F^,  F j,  F(  as 
controls  and  X,  Y  axes  parallel  to  the  sides  of  the  photograph.  C  Is 
the  center  of  tne  family  of  hyperbolas  and  all  points  on  a  specific 
hyperbola  will  have  the  same  correction.  Pj  and  will  have  unnecessarily 


large  corrections. 
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Another  point  that  should  be  considered  is  that  if  the  determinant 
of  the  matrix  is  not  exactly  zero,  but  close  to  zero,  there  will  be  large 
errors  In  the  solutions.  Such  a  case  can,  in  general,  be  detected  by 
computing  the  remaining  errors  after  correction  In  the  coordinates  of  the 
flduclals  themselves.  Now,  in  all  three  cases  tried  the  center  of  the  family 
of  hyperbolas  was  close  to  the  center  of  the  photograph, and  so  the  side 
flduclals  or  our  controls  were  all  close  to  one  or  the  other  asymptote.  On 
the  other  hand,  It  was  observed  that  the  corner  points  of  the  photographs 
usually  had  the  same  order  of  absolute  magnitude  of  error  before  correction 
as  the  side  points. 

So,  when  these  side  flduclals  which  are  close  to  the  asymptotes  are  used 
to  obtain  the  family  of  hyperbolas,  a  point  which  is  far  from  both  asymptotes, 
o.q.,  a  corner  point  will  be  associated  with  a  very  large  f,  whereas  the 
true  f  might  be  the  same  order  of  magnitude  as  Tor  a  side  point.  Therefore, 
when  side  flduclals  ore  used  as  controls,  unnecessarily  large  amounts  of 
correction  are  made  on  all  points  whose  XV  are  large  and,  hence,  poor 
Interpolation  results  ore  obtained. 

On  the  other  hand,  when  coiner  flduclals  are  used  as  controls  and  the 
axes  ore  parallel  to  the  sides  of  the  photograph,  all  the  flduclals  cannot 
be  close  to  the  asymptotes;  therefore,  the  XV  of  the  corner  points  being 
usually  the  largest  of  any  point  on  the  photograph,  the  family  of  hyper¬ 
bolas  are  so  defined  that  no  point  on  the  photograph  will  be  associated  with 
an  unusually  large  correction  and  so  better  results  are  obtained. 


ifcvuiaiw n.'ir,'.  [<nn4. - 
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C.  Conclus Ions 

1.  The  data  available  was  Insufficient  to  draw  any  finalized 
conclusions. 

2.  As  film  distortion  Is  nearly  a  random  process,  there  is  a 
limit  to  the  accuracy  of  any  Interpolation  method  that  may  be 
used.  On  the  average,  this  accuracy  will  increase  as  the  number 
of  control  points  becomes  greater. 

3.  Defining  >  as  "Is  better  than",  the  data  available  has 
indicated  that: 

Tewlnkel's  formula  >  Linear  Interpolation  >  Linear  correction 
|n  X  and  Y  ^Uniform  linear  correction. 

h.  The  corner  flduclals  gave  better  results  than  side  flduclals 
In  all  tests  run  on  the  three  sets  of  data  available,  However, 

In  methods  (a),  (b) ,  and  (cl  there  was  not  any  significant 
difference  between  the  results  to  Indicate  0  general  pattern. 

5.  When  using  Tewlnkel's  formula,  to  assure  reliable  results  the 
axes  should  be  parallel  to  the  sides  of  the  photograph  when  the 
corner  flduclals  are  used,  and  the  axes  should  be  parallel  to  the 
diagonals  when  the  side  flduclals  are  used.  Furthermore,  there 
seems  to  he  less  round-ofr  error  when  the  center  of  the  photo¬ 
graph  Is  used  as  the  origin. 


Page  53 


THE  SOLUTION  OF  SIMULTANEOUS  EQUATIONS 


A .  I  ntrociuet  Ion 

A  general  study  of  different  methods  of  solutions  of  simultaneous 
equations  was  made  to  determine  whether  other  methods  might  be  more 
efficient  than  the  Stlfel  method  of  Iteration  for  the  solution  of  the 
normal  matrix  of  aerotrlangulatlon. 

As  a  result  of  this  study,  the  Jordan  diagonal Izatlon  method  has 
been  programmed  for  the  object  program  and  In  Its  final  form  con  be 
used  when  there  are  less  than  25  photographs  In  the  aerotrlangulatlon 
problem.  This  Is  discussed  further  In  Appendix  D. 

We  shall  use  the  following  notation  for  the  set  of  slmu! taneous 
equations  that  arc  to  bo  solved. 


1,1 

CM 

<T 

’•  A!,N-1 

A1  ,N 

2,1 

A„  o  • 
2,2 

“  A2,N-1 

A2,N 

N ,  1 

A 

, ,  A.  . 

A 

N,2 

N,N-1 

N,N 

lAl'lXl  » 


!  j 

Vi  AN,2  AN,N-1  AN,N  j  XNj  1  CN 

Aj  j's  are  the  equation  coefficients,  x^'s  are  the  unknowns,  and  Cj’s  are 
the  constants  associated  with  the  equations. 

B .  Direct  Methods 

I .  Jordan  D laeonal I zat ion  Method 
a.  Restrictions: 

Matrix  I A  t  must  have  non-zero  diagonal  elements.  The  program 


will  stop  If  any  of  the  diagonal  elements  become  zero  during  the 
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RM  JORDAN  DIAGONAL  I ZAT1 ON  METHOD 
RM  FOR  SOLUTION  OF  SIMULTANEOUS  EQUATIONS 

RM  650  BUFfcTRAN  SOURCE  PROGRAM 

DIMEN  A ( 27 * 28  > 

RM  N  MUST  BE  LESS  THAN  28 

12  READ.N 
M  =  N+1 

DO  1  I ~ 1 • N 
DO  1  J*1.M 
1  READ.  A ( 1 • J ) 

DO  4  1=1»N 
R  =  A(  1 .1 ) 

J=I+1 

DO  13  L  =  J.M 

13  A(I»L)*A(I«L)/R 
DO  4  K=1.N 

IF( K-I 18.4.8 

8  P  =  A  <  K  » I ) 

DO  14  L=J.M 

14  A(K»L)=A(K.L)-A( 1 »L)*P 

4  CONTINUE 

DO  9  J«1.N 

9  PUNCH.  A(J.M) 

GO  TO  12 

END 
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solution.  This  Is,  In  general,  quite  likely  to  happen, 

b.  Input  Format: 

First  card  contains  N,  the  number  of  equations,  then  the 
elements  of  matrix  A,  and  vector  C  are  read  In,  one  word 
per  card,  In  the  following  order: 


A|,l'  A 1 , 2 ’  ”•  Al,n-1 

A„  , ,  ,  C„,  ... 

2 ,  n- 1 ’  2,n  2’ 

C  . 


AH,r  A2,2’  “■ 

...  A  ,  A  , 
n,n-l  n,n 


n 

c.  Output  Format: 

Solutions  to  the  unknowns  Xf  are  punched  one  word  per  card 

In  the  following  order:  X,,  X  ,  ...  X  . 

12  n 

d,  Running  T Ime: 

The  program  takes  about  50  seconds  on  the  IBM  650  when  N, 
the  number  of  equations,  Is  four.  Forty-five  seconds  are 
required  on  the  IBM  1620  for  N  equal  to  eleven. 

2 .  Jordan  Diagonal Izat Ion  Method  for  Symmetric  Matrix 

This  Is  essentially  the  same  as  the  previous  program  except 
that  It  makes  use  of  the  symmetry  of  matrix  A  and  has  fewer 
operations . 

a.  Restrictions: 

Matrix  A  must  be  symmetric  and  must  have  non-zero 
diagonal  elements.  The  program  will  stop  If  a  diagonal  element 
becomes  zero  during  the  solution.  This  is,  In  general,  quite 
I Ikely  to  happen. 

b.  Input  Format: 

The  first  card  contains  N,  the  number  of  equations.  Then 
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RM  JORDAN  METHOD  FOR  SOLUTION  OF  A  SYMMETRIC  MATRIX. 
RM  650  BUFFTRAN  SOURCE  PROGRAM 

DIMEN  A( 24»25 ) »B( 25  > 

RM  N  MUST  BE  LESS  THAN  25 

12  READ  »N 
M»N+1 

DO  1  1  =  1. N 
DO  1  J=1*M 
1  READ.A(I.J) 

DO  4  I  =  1  *  N 
J»I  +  1 
R*A< 1  .1) 

DO  13  L  =  J  *M 
B  (  L  )  = A ( I »L) 

13  A( I *L ) = A( I .L) /R 
DO  4  K=1»N 

I  F  (  K-l >8.4.5 

8  P«A  (  K.  *  I  > 

DO  14  LB  J  *M 

14  A(K»L)=A(K»L>-A( I »L)*P 
GO  TO  4 

5  P-BUJ 

DO  7  L*K  »M 

7  A(K»L)-A(K*L)“A(  I  »L)*P 

4  CONTINUE 

DO  9  J=  1  »N 

9  PUNCH. A(J.M» 

GO  TO  12 
END 
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the  elements  of  matrix  A  and  vector  C  are  read  In  one 
word  per  card  In  the  following  order: 

A1,l*  A1 ,2*  A1,n-I’  Al,n’  CT  A2,|’  A2,2‘•• 

A2,n-I*  A2,n’  C2*  An,I’  An,2’  An,n-1* 

A  ,  C 
n ,  n  n 

c.  Output  Format: 

Solutions  to  the  unknown  are  punched  one  word  per 
card,  In  the  following  order:  X ^ ,  X^,  ...  X^. 

d.  Running  Time: 

The  program  takes  about  '0  seconds  on  the  IBM  650  when  the 
number  of  equations  Is  four.  No  other  estimate  Is  available, 
but  It  Is  faster  than  the  general  Jordan  Diagonal  I zat I  on  method 
discussed  previously. 


3. 


This  program  Is  not  very  efficient  and  Is  not  particularly 
recommended  if  the  number  of  equations  Is  large.  It  Is  being 
mentioned  here  as  It  was  programmed  In  the  preliminary  studies  and 
may  be  used  to  compute  the  Inverse  of  a  matrix. 

a,  Restrictions: 

The  diagonal  elements  must  be  non-zero. 

b.  Input  Format: 

First  card  contains  N,  the  number  of  equations.  Then 
the  elements  of  matrix  A  and  vector  C  are  read  In  in  the 


fol lowing  order: 


RM  SOLUTION  OF  SIMULTANEOUS  EQUATIONS 
RM  BY  INVERSE  MATRIX  METHODS 

RM  650  BUFFTRAN  SOURCE  PROGRAM 

DIMEN  A(12»12) *8(12*12) *C<12) 

1  X<  12) 

12  READ*N 

RM  N  MUST  BE  LESS  THAN  13 

M=N  +  1 

DO  1  1=1. N 
DO  1  J=1«N 
READ*  A ( I  * J ) 

I F  (  I  -  J )  2  *3 . 2 

2  B(I*J)=0.0 
GO  TO  1 

3  B(  I  *  J )  a  1  •  0 

I  CONTINUE 

DO  15  J=1*N 
15  READ.  C(J) 

DO  4  I = 1 *N 
R=A ( I . 1 ) 

DO  13  L=1.N 
B(I.L)-B(I.L)/R 

13  A! I *L ) “A( I *L) /R 
DO  4  K=  1  * N 

P"  A  (  K  *  I ) 

7  I F  (  K" I  )  8 1 4  *8 

tt  DO  14  L=1.N 

R(K*L)=B(K*L) *  B( 1 »L)*P 

14  A(K*L)=A(K»L) *  A ( I *L ) *P 

4  CONTINUE 

IF  < M-2 • 0)10*11. 10 
10  DO  9  I = 1 *N 

DO  9  J=1*N 
9  PUNCH.  BU.J) 

II  DO  5  1  =  1.  N 
X(  I  )=0.0 
DO  5  J=  1 . N 

5  X( I )=B( I *J)*C< J)+X( I > 

DO  61  =  1  *N 

6  PUNCH . X ( I ) 

GO  TO  12 


•  •  < 
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c. 


A,,r  Al ,2*  ’  Al,n-1*  A1,n'  A2,|*  A2,2’ 

A2,n-J*  A2,n . An,l’  An,2’  ’  An,n-r 

An,r,’  C1’  C2 . Cn* 


c.  Output  Format : 

If  we  denote  the  elements  of  the  Inverse  matrix  of 

2 

A  by  B  then  the  first  N  cards  contain  the  elements 
'  iJ 

of  the  Inverse  In  the  following  order: 


B 

B 

B 


,n-l 

.  B.  . 

B  , 

,  B 

’  1." 

2,1 

Bn.l 

’  8nl2’ 

.  B 

2,2* 
n , n- 1  ’ 


The  last  N  cards  contain  the  unknown  X^  In  the  following 

order:  X. ,  Xt,  , . ,  X  . 

I  i  n 

d.  Running  Time: 

The  program  takes  about  one  minute  and  10  seconds  on  the 
IBM  650  for  N  equal  to  four. 

Iterative  Methods 
1 .  Sn I  del  Mefhod  of  Iteration 


a.  Restrictions: 

Matrix  A  has  to  be  symmetric  and  positive  definite.  A 
set  of  normal  equations  has  these  propert ies,  and  therefore 
It  Is  possible  to  use  this  method  to  solve  the  normal  matrix 
of  the  aerotr i angulat Ion  problem. 

b.  Input  Format: 

To  use  this  program  matrix  A  and  the  constant,  vector  C 
should  have  already  been  stored  in  the  computer.  P  oper  storage 
can  be  achieved  by  using  the  program  for  formation  of  normal 
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RM  SIEDEL  S  METHOD  OF  ITERATION  FOR  A  SYMMETRIC 
RM  MATRIX  TO  BE  USED  ALONG  AND  AFTER  THE 
RM  PROGRAM  FOR  FORMATION  OF  NORMAL  EQUATIONS 
RM  650  BUFFTRAN  SOURCE  PROGRAM 

DIMEN  A ( IB  •  19 ) »Y( 1 B  »  ,X(18) 

RM  N  MUST  BE  LESS  THAN  19 
READ.N.R.P 
L=N+ 1  * 

DO  3  I ~ 1 • N 
3  Y  <  I  )  =  A  (  I  ♦  L  ) 

IF(R-5. 15.6.6 
6  DO  8  I = 1 • N 

8  READ  *  X ( I ) 

GO  TO  20 

5  DO  »  I“liN 

9  X(I),=  0. 

20  M= 1 • 

H  =  0. 

10  DO  12  K=1.N 
G»0. 

DO  11  1=1. N 

11  G*G+A(K.I ) *X ( I ) 

Q*(Y(K)-G)/A(K»K) 

X  (  K.  1  *X  (  K )  +Q 
IF(K-M) 12.13.12 

13  lh(U-P> 14.14. 12 

14  M=M+1 

I F ( M-N  >12.12.17 

12  CONTINUE 
H  =  H+1 

GO  TO  10 

17  DO  18  1=1. N 

18  PUNCH. X( I) 

PUNCH. H 
STOP 

END 


Page  61 


equations  described  below. 

This  program  itself  reads  at  least  one  card  on  which  there 
should  be  three  values  designated  by  N,  ft,  and  P  respectively. 

N  Is  the  order  of  matrix  A  .  If  R  5.0,  It  will  be  assumed 
that  there  are  no  initial  approximations  available  to  the 
unknowns,  and  the  program  will  start  from  Initial  approximations 
of  zero  to  the  unknowns  and  no  more  cards  will  be  read.  If 
R  5*0,  It  will  be  assumed  that  there  are  some  initial 
approximations  available  to  the  unknowns,  and  the  program  will 
read  N  cards  containing  the  approximations  to  X ^ ,  X^,  ...  X^ 
respectively.  P  should  be  In  the  order  of  magnitude  of  the 
maximum  error  permissible  In  any  of  the  unknown. 

c.  Output  Format: 

l hen  the  error  in  all  the  unknowns  Is  in  the  order  of 
magnitude  of  P  (described  above),  N+l  cards  will  be  printed 
out.  The  first  N  cords  will  contain  the  unknowns  Xj,  X^,  ... 

X^  respectively.  The  last  card  will  contain  the  number  of 
iterations  the  program  has  gone  through, 

d.  Running  Time: 

Running  time  depends  entirely  on  the  matrix  A  ,  the 
initial  approx imat ions  used,  and  the  accuracy  desired. 

This  program  is,  tn  general,  efficient  only  for  Improving 
the  results  of  a  direct  method  of  solution. 

Stlefei  Method  of  Iteration  with  Plrac  Density  Function 
a.  Restrictions: 

Matrix  A  Has  to  be  symmetric  and  positive  definite.  It 


minim,  w  ir  hi  f * 
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RM  5TIEFEL  METHOD  OF  ITERATION  WITH  DIRAC  DENSITY  FUNCTION 
RM  TO  BE  USED  ALONG  AND  AFTER  NORMAL  EQUATION  PROGRAM 
RM  650  BUFFTRAN  SOURCE  PROGRAM 

DIMEN  A< 18*19 ) »X( 16 ) *DPX(  18)  »R 
1  (18! *DPR( 1 8 ) *  AR ( 18) 

RM  N  MUST  BE  LESS  THAN  19 
READ.N.NUM 
L  =  N+ 1 
Q=0. 

CP*1. 

DO  6  I = 1 ♦ N 
X( I ) =0. 

dpx  ( n=o. 

R ( 1 ) =A( I »L ) 

6  DPR ( 1 ) =0 • 

3  H"0  • 

2  H  =  H  + 1 

DO  7  I = 1 , N 
AR ( I )=0. 

DO  7  J  =  1  «  N 

7  AR( 1)=AR( I >+A( 1  .J)*R< J) 

C  =  0. 

DO  8  I  =  1  .  N 

8  C«C+AR(I)#RII) 

CA-0« 

DO  9  I  =  1  *  N 

9  CA=  CA  +  AR ( 1 ) * AR ( 1 ) 

P  =  <C*0) /CP 

QaC A/C-P 

DO  12  I = 1 »N 

DXS ( R ( 1 )+P*DPX ( I ) ) /Q 

X( I  )=X<  I ) +DX 

DPX (  I  )  aDX 

DR  = ( P*DPR ( I )-AR( I )  ) / 0 
R(  I  )=R(  I)+DR 
12  DPR (I) -DR 

CP=>C 

I F ( H-NUM ) 2 ♦ A  * A 

4  DO  5  I=1»N 

5  PUNCH  »  X (  I ) »K(  1  ) 

STOP 

END 
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should  be  noted  that  a  set  of  normal  equations  has  these 
properties. 

b.  Input  Format: 

To  use  this  program,  matrix  A  and  the  constant  vector 
C  should  have  already  been  stored  in  the  computer.  Proper 
storage  can  be  achieved  by  using  the  program  for  formation 
of  normal  equations  described  below. 

This  program  Itself  reads  one  card  containing  the  two 
values  N  and  NUM  respectively.  N  Is  the  order  of  matrix  A  , 

NUM  is  the  number  of  iterations  desired  -  It  should  In  general 
be  a  multiple  of  N. 

c.  Output  Format: 

When  tho  program  has  gone  through  the  desired  number  of 
Iterations,  It  will  punch  N  cards,  each  containing  two  values. 

The  first  value  on  the  1th  card  will  be  the  unknown  x  ,  and 
the  second  value  will  be  the  residual  associated  with  the 
I  th  eqoat ion. 

d.  Running  Time: 

Running  time  depends  on  the  accuracy  desired  and, 
consequently,  on  the  number  of  Iterations  for  four  equations. 

Each  set  of  four  iterations  takes  about  thirty  seconds  on  the 
IBM  650,  For  four  equations,  at  least  eight  Iterations  were 
necessary  to  obtain  an  accuracy  of  10’^. 

3 .  Stlefel  Method  of  Iteration  Using  Chebvshev  and  Dirac  Density  Funtlons 
This  method  is  essentially  the  same  as  program  8,2  discussed 
above,  the  only  difference  being  that  it  computes  some  initial  approx- 
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Imatlons  to  the  unknowns  before  starting  that  program.  This  process 
consists  of  three  separate  programs  which  are  to  be  used  one  after 
the  other.  Data  common  to  more  than  one  program  is  permanently 
stored  and  carried  on  from  one  program  to  another  by  proper 
dimensioning . 

4.  Determination  of  the  Smallest  and  Largest  Eigenvalues  of  Matrix  A 
This  program  uses  a  very  crude  method  for  the  determination  of 
the  largest  eigenvalue  of  matrix  A  and  then  takes  an  assigned 
fraction  of  It  as  the  smallest  eigenvalue, 
a.  Restrictions 

Matrix  A  has  to  be  symmetric  and  positive  definite. 

I).  Input  format: 

To  use  this  program,  matrix  A  and  the  constant  vector  C 
should  have  already  been  stored  In  the  computer.  Proper  storage 
can  be  achieved  by  using  the  program  for  fo  matlon  of  normal 
equations,  described  below. 

This  program  Itself  reads  one  card  containing  two  values 
denoted  by  N  and  C  respectively.  N  Is  the  order  ot  matrix  A  . 

C  Is  the  fraction  of  the  largest  eigenvalue  that  is  supposed 
to  be  considered  as  the  smallest  eigenvalue,  A  value  of 
C  =  n,  4  or  0.02  usually  can  give  a  reasonable  approximation 
to  the  smallest  eigenvalue.  The  better  this  approximation 
is,  the  better  the  final  results  of  solution  will  be- 
c.  Output  Format: 

This  program  prints  out  one  card  containing  the  computed 

maximum  and  minimum  eigenvalues  respectively. 
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RM  MAXIMUM  MINIMUM  ROUTINE 

RM  PROGRAM  FOR  DETERMINAT ION  OF 

RM  THE  GREATEST  AND  SMALLEST 

RM  EIGENVALUES  OF  A  GIVEN  MATRIX 

RM  TO  BE  USED  AFTER  THE  PROGRAM  EUR  NORMAL 
RM  EQUATION  FORMATION  AND  BEFORE  STIEFEL 

RM  650  BUFFTRAN  SOURCE  PROGRAM 

DIMEN  A ( IB, 19) »X( 18)  »Y(  IB) , 

1  NN '  1)  , M  A  X  A ( 1 ) »MI NA ( 1 ) 

RM  N  MUST  BE  LESS  THAN  19 

RIEAD  »N  *  C 
MAX-0, 

DO  1  I  “  ]  *  N 
M IN-O, 

DO  2  J -  1 » N 
2  MIN«MIN+A( 1 ,J) 

I F ( MAX-M I N ) 5 , 1  *  1 
5  M  A  X  =  M I  N 

1  CONTINUE 

MIN=C*MAX 
M  A  X  A ( 1 ) “MAX 
MI  MAI  1 ) - m I M 
NN(  1  )  »M 
PUNCH ,MAX«M1N 
STOP 
END 


d.  Running  Time: 

Running  time  Is  around  five  seconds  on  the  IBM  650 
for  four  equations. 

Stlefel  Method  of  Iteration  with  Chebyshev  Density  Function 

This  method  of  Iteration  by  Itself  is  not  very  efficient  and 
Is  used  only  to  compute  Initial  approximations  for  the  Stlefel 
method  with  Dirac  Density  Function  which  will  be  explained  later. 

a.  Restrictions: 

Same  as  above. 

b.  Input  Format: 

This  program  reads  only  one  card  containing  the  number 
of  Iterations  desired.  This  number  should  In  general  be  a 
multiple  of  N.  The  matrix  A  and  vector  C  should  already 
have  been  stored.  Proper  storage  will  be  achieved  If  the 
program  above  has  been  used  correctly. 

c.  Output  Format: 

This  program  does  not  have  any  output.  Relevant  Inform¬ 
ation  is  stored  In  memory  to  be  used  by  the  program  next 
descr 1  bed . 

d.  Running  Time: 

An  accurate  running  time  Is  not  available.  But,  it  Is 
expected  that  for  N  equal  to  four,  using  four  iterations, 
about  thirty  seconds  would  be  required  on  the  IBM  650. 
Stlefel  Method  of  Itorat ion  with  Dirac  Density  Function 

Except  for  the  Information  obtained  and  used  from  the  pro- 
•iram  above,  this  pronram  is  the  same  as  program  C,2. 
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RM  STIEFEL  METHOD  OF  ITERATION  WITH  CHEBYSHEV 

RM  DENSITY  FUNCTION  TO  BE  USED  AFTER  THE  MAX 

RM  MIN  ROUTINE  AND  BEFORE  STIFFEL  DIRAC  DENSITY 
RM  650  BUFFTRAN  SOURCE  PROGRAM 
DIMEN  A< 18*19) »X( 18) * Y ( 18)  ♦ 

1  NN ( 1 ) tMAXA ( 1 ) »MI NA  1 1 )  » 

2  DPXl  18) tR( 18) 

RM  NUM  MUST  BE  LESS  THAN  19 
READ  »NUM 
M  A  X  =  M  A  X  A  (  1  ) 

M  1  N  =  M  1  N  A  (  1  ) 

N  =  NN(  1  I 
L  =  N  +  1 

MMMsMAX-MIM 
MPMbMAX+M I N 

BAR“MPM+?*SQRT(MAX*MIN) 

BARDSBAR/MMM 
RARD1  =  !JARD 
RABD=MMM/I3AR 
RAUDl=RADD 
DO  1  1  =  1  *  N 
Y ( 1  ) = A  I  l.L) 

R( I ) =Y( I ) 

DPXl 1 ) =  < I) /MPM 

1  xmwM’xm 

CSIINP=MI>M/MMM 

C05HN-1. 

1-1=0. 

6  IT=tl+ 1 

DO  3  I  1  * N 
YESTI =0. 

DO  6  J  =1  *  N 

5  YEST1-YESTI+A1 I *  J  1  *  X  (  J  ) 

3  R( I  )=Y( I I-YESTI 

CSHNM=COSHN 
COSHN=CSHNP 
I3ARD]=BARD1«BARD 
RABD1=RABD1*RABD 
CSHNP- I BARD1+RABD1 ) / 2 
DO  2  I = 1 .N 

DPXl I ) *  1 4*C0SMN*R (  I ) /MMM+CSHNM 
1  #DPX ( I ) ) /CSHNP 

2  XI  M  =  X(  I  )  +  D P X  (  I  ) 

IF ( H-NUM ) 6  »  7  »  7 

7  STOP 
END 
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RM  STIEFEL  ITERATION  WITH  DIRAC  DENSITY 

RM  TO  BE  USED  AFTER  THE  STIEFEL  WITH 

RM  CHEBYSHEV  DENSITY  FUNCTION 

RM  650  BUFFTRAN  SOURCE  PROGRAM 

DIMEN  A! 18.19) .X( 18) ,Y( 18) ,NNt 

1  1) ,DPX( 18) »RI 18) »DPR( 18)  , 

2  AR ( 18) 

RM  NUM  MUST  BE  LESS  THAN  19 

READ  .NUM 
N"NN ( 1 ) 

L  =  N+1 
(3=0. 

CP=1. 

DO  6  IM.N 
DPR ( I ) =0. 

DPX ( I ) =0. 

YEST I =0. 

DO  11  J=>  1  .  N 

11  YESTI=YESTI+A( I » J ) *X l J ) 

6  R ( I )=Y( I ) -YEST I 

3  H  =  0  • 

2  H=H+ 1 

DO  7  I « 1 .  N 
AR ( I )e0. 

DO  7  J= 1 »N 

7  AR< I )=AR( !)+A( I .J)*R< J) 

C  =  0. 

DO  8  I = 1 » N 
0  (>C+AR<!)*R(I> 

CA=0. 

DO  9  I=1»N 

9  CA=CA+AR ( I ) *AR (1) 

P*< C*Q) /CP 

Q=C A/C-P 

DO  12  I  =  1  *  N 

D  X  =  <  R  ( I )+P*DPX( 1)1/0 

X ( l ) =X ( I ) +DX 

DPX ( [ ) =DX 

DR= (P*DPR( I ) “AR ( I ) ) /U 
R ( I ) =R ( I ) +DR 
12  DPR ( I ) =DR 

CP»C 

I F ( H-NUM 12.4*4 

4  DO  5  1=1»N 

5  PUNCH  »  X (  I  )  .  R (  I ) 

STOP 

END 
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v-  a.  Restrictions: 

t ' 

I,  Same  as  a.  above 

f-  b.  Input  Format: 

j  This  program  reads  only  one  card  containing  the  number 

\  of  iterations  desired.  This  number  should  In  general  be  a 

s  '  multiple  of  N.  Matrix  A  arid  vector  C  would  already 

;  have  been  stored  properly  If  the  two  previous  programs  have 

been  run  correctly. 

c.  Output  Format: 

When  the  program  has  gone  through  the  desired  number  of 
Iterations,  It  will  punch  M  cards  containing  two  values.  The 
I'  •  first  value  on  the  Ith  card  will  be  the  unknown  and  the 

■'  second  value  will  bo  the  residual  associated  with  the  Ith 

equation. 

d.  Running  Time: 

1.  Running  time  Is  about  thirty  seconds  on  the  IBM  6!>0 

“  '  for  M  equal  to  four  and  four  Iterations. 

L  0 ,  Formation  of  Normal  liquations 

:  This  program  can  be  used  to  form  0  set  of  normal  equations  from  a 

redundant  or  non- redundant  sat  of  equations.  It  will  store  the  resultant 
set  of  normal  equations  In  the  proper  storage  space  required  by  several 
of  the  previously  described  programs. 

We  shall  consider  the  set  of  possibly  redundant  number  of  equations 


In  the  following  form: 


.Upp,  liWIU 
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RM  FORMATION  OFNORMAL  EQUATIONS 

RM  650  BUFFTRAN  SOURCE  PROGRAM 

RM  N  MUST  BE  LESS  THAN  19 

RM  M  MUST  BE  LESS  THAN  26 

DIMEN  A( 18*19) »B125. 19) 
READtNfM 
L  =  N+ 1 

DO  1  1 = 1 *M 
DO  1  J  =  1  *  l. 

1  RE AD  »B  < 1 »J) 

DO  5  I = 1 • N 
DO  5  J= 1  * L 
A(  I  * J)=0. 

DO  5  K- 1 *M 

5  A  (  I  *  J  )  -  A I  I  .  JHOUt  J)*U(K.I  ) 
GO  TO  CONSW 
END 
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a.  Input  Format: 

This  program  first  reads  one  card  containing  N  arid  M, 
then  reads  M  x  (Nil)  cards  containing  the  following  values 
respect  1 vely : 


Yr  bi,2 . Bi,N-r  bi,n'  cr  Yr  02,2’  ••• 

02  ,N"  I'  02,N‘  C2 . BM,r  BH ,  2  ’  ’  3M,N-|' 

0M,N*  CM* 

b.  Output  Format: 

This  program  does  not  have  any  output.  It  Just  stores 
the  resultant  matrix  and  constant  vector  in  the  following 
way.  If  we  denote  the  resultant  set  of  equations  by  out 
previous  notation  of  matrix  A  and  vector  C  ,  then  the 

element  A  will  be  In  the  location  specified  by  A(|,J) 

1  *J 

In  Fortran.  The  component  C|  of  vector  C  will  be  in  location 
specified  by  A(l,N+l)  In  Fortran. 

c.  Running  Time: 


This  program  takes  roughly  20  seconds  on  the  IBM  6£0 
when  both  M  and  N  are  four. 

E .  Remarks 


I.  All  the  programs  discussed  above  were  written  solely  for  exper¬ 
imental  purposes.  They  were  essentially  Intended  to  be  used  with 


a  small  number  of  equations;  and,  therefore,  their  Input  and  output 
formats  are  In  no  way  claimed  to  be  the  most  efficient  possible. 

2.  It  should  be  noted  that  In  IBM  650  Bufftran,  which  was  the 
coding  language  used,  all  quantities  are  considered  to  be  floating 
point  quantities.  Therefore,  It  might  be  necessary  to  change  some 
of  the  notations  In  the  programs  If  they  are  to  be  run  on  a 
standard  Fortran  system. 

Conclusions 

1.  The  Jordan  D 1  agonal  1 zat ion  and  the  Grout  methods,  which 
are  quite  similar,  appear  to  be  the  simplest  and  fastest  direct 
methods  for  solution  of  simultaneous  equations. 

2.  The  Inverse  matrix  method  Is  not  very  efficient  because  of 
Its  length  and  Its  Inaccuracy. 

3.  The  single  step  and  total  step  Iterations  arc  not  very 
efficient  because  of  the  restrictions  for  convergence. 

4.  The  Scldcl  method  of  Iteration  for  symmetric  and  positive 
definite  matrices  converges  'apidly  to  a  high  accuracy  when 

a  good  Initial  approximation  Is  used.  However,  convergence  is 
rather  slow  when  poorer  Initial  approximations  are  used  although 
convergence  Is  guaranteed.  Therefore,  In  cases  where  the  set  of 
simultaneous  equations  have  been  solved  by  a  direct  method  and  the 
results  are  in  error  due  to  round-off,  the  Seidel  method  may  be 
used  to  Improve  the  results. 

5.  The  Stlefel  method  with  Dirac's  Density  Function  converges 
quite  rapidly  for  all  matrices  tested,  and  after  2N  Iterations,  N 
being  the  number  of  equations  and  In  this  case  4  or  8,  the  maximum 
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relative  error  was  In  the  order  of  lo"^. 

6.  The  Stiefel  method  with  Chebyshev  Density  Function  proved 
to  give  an  initial  approximation  for  the  Stiefel  method  with 
Dirac's  Density  and  eventually  gave  a  better  result  than  when  an 
initial  approximation  of  zero  was  used.  But,  this  Improvement  still 
was  in  the  order  of  10  ^  which  suggests  that  the  Stiefel  method 
with  Chebyshev  Density  might  not  be  necessary  at  all. 

7.  Although  the  accuracy  achieved  by  a  direct  method  might  take 
a  longer  time  by  an  iterative  method,  In  general  iterative  methods 
have  an  advantage  as  far  as  storage  space  is  concerned.  This  Is 
because  Iterative  methods  In  general  do  not  operate  on  the 
matrix  of  simultaneous  equations,  i.e.,  they  do  not  change  It. 
Whereas  direct  methods  operate  on  the  matrix,  and  non-zero  elements 
might  be  created  where  initially  there  were  zeros  which  did  not 
require  storage  allocation. 

8.  The  accuracy  achieved  by  the  Jordan  Diagonal izat Ion  method 
definitely  requires  a  considerably  longer  time  to  be  achieved  by 
the  St lef e!  method. 

A  Modified  Segment  3  has  been  written  to  be  used  when  the 
higher  accuracy  possible  with  the  Stiefel  method  is  not  required 
and  24  or  less  photographs  are  involved  so  that  computer  time  may 
be  conserved.  The  realized  reduction  in  running  time  has  been 
essential  to  the  execution  of  this  project.  The  Modified  Segment 
3  also  represents  a  practical  implementation  of  the  solution  of  the 
aerotr iangulat ion  problem  normal  matrix. 
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VII.  MATHEMATICAL  DESCRIPTION  OF  THE  PROGRAM 
A,  PROGRAM  FUNCTION 

Given  a  certain  amount  of  photographic,  exposure  station,  and 
ground  point  Information,  the  program  finds  a  simultaneous  solution 
to  the  orientation  and  location  of  the  exposure  stations  with  regard 
to  geocentric  axes.  This  simultaneous  solution  Is  based  on  a  least 
square  adjustment  of  the  random  errors  associated  with  the  measured 
data. 

8 .  INPUT  DATA 

The  Input  data  consists  of  parts  of  the  following: 

1.  The  measured  x  and  y  photographic  coordinates  of  the 
Images  of  ground  points. 

2.  The  geographic  coordinates  (latitude,  longitude,  elevation) 
of  ground  points  and  exposure  stations. 

j.  Estimates  of  the  geographic  coordinates  and  orientation 
of  the  exposure  stations. 

4.  Focal  length  of  the  camera,  major  and  minor  axes  of  the 
ellipsoid  defining  the  geoid. 

C,  THE  EQUATIONS 

They  can  be  divided  Into  two  groups: 

I  .  Ground  point  equations 

These  equations  express  either  the: 

a.  Intersection  of  rays  corresponding  to  the  images 
of  a  ground  point;  or 

b.  enforcement,  when  they  are  known,  of  the  geographic 
coordinates  of  a  ground  point. 
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2 .  Exposure  station  equations 

These  equations  enforce,  when  they  are  known,  the 
geographic  coordinates  of  an  exposure  station. 

0.  UNKNOWNS 

These  are  of  two  types: 

1 .  Exposure  station  unknowns 

These  are  the  only  unknowns  required  for  the  solution  of 
the  pho t oc) r amine tr Ic  problem.  There  are  six  unknowns  associated 
with  every  exposure  station.  Three  arc  corrections  to  the 
estimated  geocentric  coordinates  of  the  perspective  center  (dRx, 
dRy,  dR>;) ,  and  the  remaining  three  are  the  components  along  the 
geocentric  axis  of  a  differential  rotation  vector  dT  which  Is 
the  correction  to  the  orientation  of  the  camera  (dT^,  cIT^,  dT^) . 

2 .  Residual  unknowns 

These  unknowns  are  not  part  of  the  final  solution  of  the 
photogrammetr Ic  problem.  Their  introduction  Is,  however, 
essential  for  the  solution  of  the  set  of  simultaneous 
equations.  The  reason  Is  the  following: 

In  general,  the  available  data  is  redundant  In  that  It 
gives  more  Information  that  necessary  to  solve  the  photogrammetr Ic 
problem.  Were  It  nor  for  errors  In  the  data,  any  part  of  It 
sufficient  to  solve  the  problem  could  be  used  and  the  rest 
discarded.  We  consider,  therefore,  that  the  true  value  A  of  a 
quantity  is  its  measured  value  Aq  plus  an  unknown  error  dA,  so 

A  =■  A  +  dA 
o 


that 


In  the  case  of  random  errors,  the  error  term  dA  has 

associated  with  It  a  standard  deviation  which  expresses 

the  probable  accuracy  of  Aq,  If  the  dA's  were  known,  the 

problem  would  be  solvable  exactly,  except  for  round-off. 

On  the  other  hand,  there  Is,  In  general,  an  Infinite 

number  of  values  for  the  dA's  such  that  the  simultaneous 

equations  are  satisfied.  The  problem  Is,  therefore,  to  find 

that  set  of  dA's  which  gives  the  most  probable  solution. 

I 

This  Is  done  by  th'c  general  method  of  least  squares  which 
makes  the  weighted  sum  square  of  the  dA's  a  minimum. 

The  residual  unknowns  Introduced  In  the  equations  are  the 
fol lowing: 


a.  d I  -  associated  with  the  x  photographic  coordinate  of 
an  Image  point 

'b.  dl  -  associated  with  the  y  photographic  coordinate  of 
^  an  image  point 

c.  d0  -  associated  with  the  known  latitude  of  a  ground 
point  or  an  exposure  station 

cl.  d,v  -  associated  with  the  longitude  of  a  ground  point 
or  an  exposure  slat  Ion 

e.  dh  -  associated  with  the  elevation  of  a  ground  point 
or  an  exposure  station 

NOTAT I  ON 

1 .  Geometr  I  c  Q.uant  1 1 1  es 


The  following  defines  basic  geometrical  quantities  of 
the  photogrammetr Ic  system. 


For  more  details  on  the  method  of  least  squares,  see  I.  and  2. 
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t 

f  - 

a. 

A  letter  with  a  bar  over  It  represents  a  vector. 

The  superscript  o  Indicate!  a  unit  vector.  Thus 
ii”  Is  a  unit  vector  along  V. 

b. 

Oj,  0? ,  0^  represent  3  consecutive  perspective  centers. 

c. 

B  =  00  =  a  vector  Joining  0„  to  0  In  the  direction 

1  '  1  from  02  to  Oj.  1  1 

d. 

B,  =  00  =  a  vector  joining  0  to  0  In  the  direction 

1  i  1  from  03  to  02.  J  1 

e . 

A  =■  vector  joining  an  exposure  station  to  a  ground  point. 

; 

f . 

A  ,  A.,  A,  «  vectors  originating  from  0( ,  0  0  to  the 

i  same  ground  point.  ■* 

r 

g- 

^1X'  ^1Y’  ^1Z  "  9eocentrlc  components  of  A. 

h. 

.0.0.0  oo  o 

A1X’  A1Y'  A|Z  or  a] X '  a i y •  or  aiz  are  geocentric  com¬ 

1  ■ 

ponents  of  A°  They  are  also  geocentric  direction  cosines 
of  A,.  1 

i . 

J. 

A)°,  Aj°,  A|°  or  a|°,  a^,  or  a^°  are  photographic  direction 

cosines  of  A^ .  They  are  also  the  photographic  components 
of  A®. 

R  is  the  vector  Joining  the  geocentric  origin  to  a  per¬ 
spective  center.  11^  refers  to  Oj. 

- 

2 .  Dor i ved  duant i t ies 

The 

following  expressions  define  certain  derived  quantities  where 

r- 

In  a  dot 

(•)  or  a  cross  (x)  between  two  vectors  represents  the  dot 

and  cross  products  respectively  of  those  vectors. 


a. 

f\  X 

A* 

8 

X 

1 

/V'  X 

K: 

‘"Z6 

1 

i_ 

'"5 

ll 

>£*> 

ij 

‘SCALAR. 

CP 

S  ' 

6  • 

1 

r- 

16 
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3/V  C 
v  ■  K*  ^ 


r  _ 

Tjf”  *  (  B  v  ft* ) 

*-  is  1  A 

£fr-  •  CB,  x  A5! 
*-  i  a 

Cj.&_  •  ("B  /  A  3 ) 

r*  i 
^  2& 

S&  •  ^  *  *y 


•  C,6  '  ^  *' 

V  7^  *  ^*,  *  Cja 

1.  c  lj%  1 

1  ^  >  (A"t  x  C  ): 

<e  rf-  4  a« 

1  __  _ 

«/=  i-  ■ (A*  *  c^‘ 

3  CJfl 

V,  --  a, 


Fi,  ■  qr[  l  c;.  *Vl] 

m,  -  Z^f  [5  6; .  a;  )  Q,  *  ^(c.,6,  A-l  j 

V  -  k[(^Al^;  ^  +_^  (J-‘  ^ 

«l'  ^  (r>.9  * 

7--i[Wa7^*  Alpj 
*v  cirp1  c^4  0-5  ^ 

*V  -  (*•  •  A^  C^  +  Cjjc^  A”  )] 

p r  rk  [-  ft "  *1) '  ft  'ft 1  ^ ft  ■  ft')  z <« 


Mj  1  c  1,2,3  J  “  1,2,3  are  the  elements  of  the  trans- 
■*  formation  matrix  [Mj  from  photographic  to  geocentric 
axes 


I  vector  joining  a  perspective  center  to  the  Image  1 


1 


t 


scalar  of  1 


I  =  1,2,3  j  =  1,2,3  are  the  elements  of  the  matrix 


Ij  defined  by: 


1*1-  sf  [&]-!£]  [Af] 


where 


o-*, 


The  superscript  (')  applied  to  the  elements  of  a  t 
matrix  Indicates  that  the  matrix  corresponds,  to  the 
second  photograph  In  a  strip,  while  the  superscript  (") 
Indicates  the  third  photograph  In  a  strip,  and  the 
absence  of  a  (')  or  (")  indicates  the  first  photograph 
of  a  strip. 

a  semi-major  axis  of  the  geo  id 
c  eccentricity  of  the  geoid 
Unknowns 

a.  Photographic  Unknowns 

These  have  been  defined  in  paragraph  0.  The  subscript 
1,2,  or  3  applied  to  dR^,  dR^,  or  dRz  indicates  that  these 
quantities  refer  to  the  first,  second,  or  third  exposure 
station  considered  in  a  particular  equation. 

b.  Residual  Unknowns 

Subscripts  1,2,  or  3  added  to  dl^  and  dl^  refer  respect¬ 
ively  to  the  images  of  the  ground  point  under  consideration 
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fc-L 
»  .. 


ic 


In  the  f I i  at ,  second,  and  third  photograph  considered 
In  a  particular  equation. 

Further  notation  will  be  introduced  as  needed. 


F.  LISTING  OF  THE  EQUATIONS 


For  the  derivation  of  the  coefficients  of  the  exposure  stations 
unknowns,  see: 

"Sequential  Presentation  of  Analytical  Aerotr iangul at  ion 
by  the  Revised  Direct  Geodetic  Restraint  Method" 

ERDL ; 

and  for  the  derivation  of  the  coefficients  of  the  residual  unknowns,  see: 
"Analytical  Aerial  Trlan  gulatlon  Error  Analysis  and 
Application  of  Compensating  Equations  to  the  General 
Block  Tr langu lot  Ion  and  Adjustment  Program" 

Phase  II  Interim  Technical  Report 
GIMRADA  Project  No,  8-35-11-520 

In  the  following,  the  condition  enforced  by  every  equation  Is  briefly 
described,  The  equations  being  generally  long,  their  listing  will  be  in 
the  fol low  I ny  way :  If,  for  example,  an  equation  Involves  two  successive 
camera  stations  and  two  Image  points,  the  unknowns  will  be: 


dTlx' 

dT  . 

'y 

dTlz' 

dRlx- 

dR,y 

dR,Z 

dV 

dT2y  ’ 

dV 

dR2x’ 

dV 

dR2z 

d,lx-  d,!y 
dl2x’  dl2y 


The  listings  will  give  the  coefficients  of  the  above  unknowns  and  the  con¬ 
stant  term  as  It  appears  on  the  right  hand  side  of  the  equation. 


EQUATION  I  -  PASS  POINT  EQUATION 


Two  camera  stations  and  the  corresponding  Images  of  a  ground 
point  are  Involved. 

The  condition  of  Intersection  of  the  rays  A,  and  A^,  from  the 
perspective  centers  to  the  corresponding  Images,  Is  expressed. 

Coef f I c lents 


<rr«  •• 

v.tA;  -  £,»*  ^ 

dT.y  : 

v,  CA,*  *  C.,a*  W 

dTlt  : 

v,  CA,*  »  C-n  ^  j 

d«.K  1 

d^,y  *. 

C-*6V 

dRu  "• 

C-n&r 

d  Tt*  1 

A£  «  C  ,a  ^ 

dTa.t  • 

K  C|B  i'y 

dT«: 

*  C,«  1  g 

dRi*  : 

-  c.;0k 

d  ! 

-  c.;,y 

dlut 

V,  CCia*  dr,,  ■*  C  l8T  4  C.  ,at  t:  3,  ) 

dXiy  1 

V,  XC.J**  "t  1 1,4  C.,*r  't  i*.  *>  C-"#*  -t  3 a. J 

dXy  ? 

'^vOCin.jt,,  4  C.,ay  t  u  4  Ctae  "tj,  ) 

: 

'A.f f*-,6x  dlii.  ^  C,ay+lx  ♦-  C >a*  t s i. 

CONSTANT  T£HM  “  B,„  C  -  Bit  C,6V  -  6  ,*  C.,»,  \ 
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2.  EQUATION  2  -  DIFFERENTIAL  SCALE  RESTRAINT 


Three  camera  stations  and  the  corresponding  Images  of  a  ground 
point  are  Involved. 

The  condition  of  Intersection  at  a  point  of  the  three  rays 
A|,  A*,  A,,  from  the  three  perspective  centers  to  the  correspond¬ 
ing  three  Image  points,  *s  expressed. 

Coef  f  1  dents 


dT 


.  .  fcs,  * ca;  •  ai  )  -  [-  cs,  •  a;)*  c  a;  •  A4»n  c(B„ 

»x  •  -  1  — - - — — - - - - - — — — - — 

J1 


drIY  :  - 


ft©,*  AMyCA,*  •  At*  1  -I-IS.-A,*)-  2 .at  t A  *  ■  A£  )]  ClBv 


} 


d  Tia: 


.  _  fee,  H  a,® )g(A;  A^)  ce,  •  a;)  *  s-a.,  ca:  .*,.•)]  c.lBi. 


t 


<£,»>*■ 


^  ^*|  >(  • 


.  ..  CA,‘-C,a)K 

tc,.)* 

CA,*  * 

>* 


4H„,  -W'£l»lt 


dTjji  - 


[-ca,  •  a* )  r  7.6 *.<.*;  -ca, *  A; )* 


} 


cat*  a;  \  -  CC8f  *i)  *  7-&u  caj •  a,*)  1  c 


/.ox 
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4Ttv  :  CA,»  -A£  )]  C,BV  -  C,jt  w  A  * )  y  7 

L  Cc.t8>ft  J 

C  &  «,  *<  At*  )Y  -  £  L  B  *,  •  Ai*  )  *  2.64  (.A  t.*  •  A»  7  J  £* 


C-CB,  ■  A>)  A*  )]C.ltE  -CQ,>A^7g  } 

CClS^  j 

CBi  *  V  >r  -  [C§t-  As*)+  la,'  CAi  ■  A.‘i]  C^ea 


Cc^*- 

dft**; 

(*Aj  •«.  c.*ja  )x 

LA,*  *•  C.,* )  * 

^,a)L 

d  1 

(,A»  *  C ijg  )y 

ca;  -  c.,,7 Y 

<C„,‘  - 

dRjjtl 

C  A3“  »  C  jj#  )  j. 

C A*  »  c,* )  E 

CC.M)fc 

dT8x; 

LCafc-  A,*)  v  la; 

CA^*  •  At*  A  5  7*  CAJ  A»*) 

<-c„  )*• 

d  r3y  ; 

C  ( 6 ». •  a 3  )  * 

la;  A,*  7]Ctfiv  -CBtAA|7YCAi  A3*  7 

tc.  «.»>*■ 

4T3fc; 

C  ceV’  a,*  )*  i.<s; 

<-*C  ■  Aj* ) 7  -(.§***,•>£  C/V  AJ) 

Ct«)‘ 

dR.3„ : 

C.AV  *  C-lS^X 

CCt.j’* 
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dftSY : 

■ 

d  Iik  1 
dXty  : 

dX  txi 
ciXlY : 
d  tSK: 
di3v : 


CA%*  *  5  tB  )  g, 

CC**)* 

-  M  w  t  |,  -  M  It  *  LI  ~  M  IX  t  31 

~  f'Ajx  *  i*.  ~  M  v(  *  t*.  ~  M  ia  ^  s*. 

CN^  f.N  j,y  “  M  it  +  4  N  x»  - 

CM  ix  "Mix)  ■*•  (  May  -  M  tY)t  u.  -*•  “ 

♦  w*.v  *-f>u  +  Mfc«t3' 


■•fill  Qx  ~  6  iy  Qy  -8|*Oa 

+  0ty  L-;y  «■  6  4.£  l-'i  J 


■  Na.*)  tai 
Mu)tn 


CONSTANT  TERM 


C^XH  *-* 
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3 •  EQUATIONS  3  and  4  -  COMPLETE  GROUND  CONTROL  EQUATIONS 

One  camera  station  and  the  corresponding  Image  of  a  ground 
point  are  Involved,  where  the  latitude^,  longitude  >,  and 
elevation  h  are  known. 

Def Ini t Ions 

A.,  Is  a  vector  to  be  defined  as  follows: 

Associate  a  component  of  A,,,  equal  to  0.7  with  the  coordinate 
that  has  the  smallest  A  component  In  absolute  value.  Associate  a 
component  of  A, ,  equal  to  0.0  with  the  coordinate  that,  has  the  Inter¬ 
mediate-  si zed  A  component  In  absolute  value  and  associate  a  component 

-  A  - 

Of  A.,  equal  to  -0.7  v^wlth  the  coordinate  that  has  the  largest  A 

component  In  absolute  value.  A„iM  and  AMMare  such  that  their 

absolute  values  are  the  smallest  and  the  largest  respectively  of  the 

absolute  values  of  the  components  of  A.  A.,  Is  perpendicular  to  A. 

A;,  Is  obtained  by  normalizing  Aol  .  A.*fc  a  ;*  7f •,  . 

Equation  3  expresses  the  condition  of  Intersection  of  the  ray  A 

with  the  vector  AJ,  passing  through  the  ground  point.  Equation  h 

expresses  the  condition  of  Intersection  of  the  ray  A  with  the  voct..-r 

passing  through  the  ground  point. 

Equations  3  and  4  condlsered  together  express  the  condition  that 

the  ray  A  passes  through  the  ground  point. 

a_, _ Equation  3  Coefficients 

:  -  vtc^oax 


Page  b6 


- 


V*A 


d*n  : 


A  a  * 


-a°wea; 


dlv  : 


“  A“a)Y  -  Anv  A0ix  ) 


^  Aotxt  „  *  Aoty  t  tj  +  Aott-t  3l ) 
vtc  CAatJl"t  ,j.  +  s.t,  + 


3  C  l  -  e.1)  siw  4  cosx 

Cl-  «*•  5.1  M*-  4  1  3/t 


+  K  Sin  ^  cos  a 


&  Cl-  a*-)  SIN  4  Sin  A 

♦  h  5,u 4 >, 


„ .  f  a  c.  i  -  e* )  cos  </>  ,  , 

■•"Aqijl  —  -  ■  ♦  h  COS  £ 

(.!-«*•  Sin  *•  4>  5/*" 


-  Ao\y 


a  tor.  4  sin  A 

f'l“crSlN14 

6  COS  4  COS  A 
(  Simv  4_ 


*■  W  CCS  <f>  Sin  'A 


T  b  COS  4  CjOS  A 


A*w  COS  4  CCS  A  ♦  Aom  ^  4  SlN  *  ♦  *\xa  *IN  4 


CONSTANT  TE.R.M  --  [-  D*A*^  -  Dy  A^y  -  0,  A*fc(t  ] 


b.  Equation  4  Coefficients 


v  U.  ^at-X 


*lTv  t  ~  v< e, 


•  "  VKS. 


*«.*  : 


d  Ry  1 


«*  • 


lAi.)*-  4-  (Aiy^1 
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~  V  IjC.  (.  Aon  ^  +  *  ^Q|»  il  ) 


^  •  ~  ^ic.  ^  A0(a  t  lt  ♦  A^iy  t tt.  r  A  om ^  si.  1 


d  <f) 


ei(,  1  ^  Sum  ^  <"«S  A 

(,  »-  e.tsiN,'4i)  ir<* 


*■  Ia  Sin  4  A 


A(  1-cMSlti  4  S|N  Tv 

- Z _ +  h  SIN  4  SIN  *k 

C 1-  c1,  s  inv  4  5  ait~ 


.  f  a  Ri-e*-  )  cas  <t 

Al«  - - - - T  ^ 


[_  c  l-cl  Sin  1  if)  )J,U 


<3  Ca\  <f)  S'  •*  "A 

Ti-a*-  Sin*-  4 


C.5S  4 


+•  h  coi  4  SlN  a 


-  a;,x 


Cx>S  4  C-®*  A  ,  . 

...  T - - -  +  V\  C.OS  4  005  a 

1  -  *-*■  sin«-  4 
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:  -  A*0(X  Ci3S  4  *os  *  -  A£iy  cos  4  sin  a  -  AJ,*  Sm  4 

CONSTANT  TtRM  =  "D*  Aj,m  ♦  Dy  A^v 
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4.  EQUATION  5  -  HORIZONTAL  GROUND  CONTROL  POINT 

Two  camera  stations  and  the  corresponding  ground  point  images 
are  involved  where  the  1 ong 1 tude  i\  I s  known. 

The  condition  that  the  common  perpendicular  between  rays  A, 
and  A^  Intersecting  ray  A,  lies  In  the  plane  of  equal  longitude, 

Is  expressed. 

Auxl 1 iarv  Def Ini t Ions 
«(,  a  3iN  T. 

/Q  «  cos 
Coef f I c lents 

dTlx  :  [<  A?*  Ox  -  fi  <  Af*  Ox  -  a,*;*  )  3 

dTlv  •,  [?CCA*K  Gy  ♦  A,  A"*  )  ~  /3  O-y^ 

dT14  ■.  l<f.  ATx  f,a-<stA?v)  6,  .  a,  A,;  )  J 

■  (.  («  A%  -  )  Lx  *  <3 

dR,y  ;  [  «  a;*  -  fl  A;v  )  Ly  -  <3  J 

d«ia  :  t<X  a;  -a  A*y  )  le  j 

dT^  :  [«.  A,\  -|3ArY)?x7 

dTtv  :  [c<  a;x  -  (Sa;y  )?y  ] 

dT^  ;  CC<AT*  -SATyjT,  ] 

d*„  :  [Ka:a  -  tfA-)  <-U7 
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f 

T; 

[<.<  ATx 

-  (3At%  )C-Ly>l 

l. 

ii* . 

6; 

t- 

•• 

CK 

-  /S  A.y  i  i-m  i  1 

r, 

d  1,^  • 

a.coc  t„ 

-<3  1*  )♦  C*A7»  -<3A%  KM„ttl  +  M,ytw  ■*•  K«tv; 

jfif  * 

<lTty 

<3,  C°<  V 

-  (3t  *  (A  -  /SAJ,  )  CM,*  t,t  ►  M,y  txvf  ^1’itja.J 

r 

dX  ; 

-  ^3  A,%  )  (N  (X  t  U  ■*  ♦  M  |«t;|i  ) 

= 

dl„  i 

CA  A,; 

-  (3  A,y  )  (  N.^tt  4  Niv*f.i.  1  Ni*1?!-5 

; 

d  Tv  • 

4  (*,x  ' 

f  d,  A*,  1  A.  L  Riy  4  A.A^) 

Constant  TE.RM  ; 


(3,  (  f'.,y  t  -<(*,*  f  A,,) 
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5.  EQUATION  6  -  HORIZONTAL  GROUND  CONTROL  POINT 

Two  camera  stations  and  the  corresponding  Images  of  a  ground 
point  arc  involved  where  the  longitude  tv  and  latttude  are  known. 

The  condition  that  the  point  of  Intersection  with  the  ray  A,, 
of  the  common  perpendicular  to  rays  A,  andAt ,  lies  In  the  prime 
vertical  plane  coorespondlng  to  4  and  *a  .  The  prime  vertical 
plane  being  the  plane  containing  the  normal  to  the  geold  at  the 
ground  point  and  tangent  to  the  circle  of  latitude  equal  to  <f> 
passing  through  the  ground  point. 

Auxiliary  Definitions 
If  coa  >  s  1  n  4 

=  Sits)  7\  SlW  $ 

£  a  -  COS  4 

5  _  ac11-  CO&4 

'mr*'  si**  4 

Coefficients 

dTtx  ;  (T  A;G*.  £  CAS,  G*  -  a,  A*tl  «■  C  CAf.G*  «•  a,  AMI 
dr,y  ;  (.1  (.  Af*  G-,  +  6,  Af,!  t  $  A,"  C-.N  .  fe  (Af,  G,  -6,  A,*»  )  1 

dT,t  U  <*»<»«  *  <3,  AM  ♦  <S  CAS,  G»  ♦  <S,  Af,  \  -  C  A%  G,  ] 

:  C  C  If  A,’, 

:  (C  $  A\  r  ^ A,*t  ■*  €  A,*,  )  lv  r  P  ") 

dfc„  ;  (U  Akt  J  ^  rtA-,)  L,*tl 


,rr> 
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W _ ! 


dT**  *• 

Cuas,* 

*  Art 

♦  t  a,; 

1  til 

dtM  : 

Cc  x  a;. 

f  <f  A,\ 

♦  t  a:c 

S  1»,  ] 

«nv.  ; 

ICl  A,V 

+  Aiv 

*  €  a;, 

)?«! 

dfttx  '• 

Cc «  a:. 

+  i'ATt 

«-  e.  a(*. 

>  c-L,n 

dR-tv 

C  os  A;* 

f  £  A,V 

)  C-Ly)  1 

r 

dRt*  • 

La  A,\  "Mr*  -  C  A,“*  )  C-Lj.l  ] 

i 

i 

r 

d£,K  : 

<3,  C5  t(i  ♦  £  tb  ♦  £  t„  )  +  C )A,\  *■  t  A,*y  ♦  €-A,\)  C  *•  M,Ttt,  *  M,t  tv 

dr1Y  : 

+  sa r,  *«  a«j  (M,**,* 

dT„ 

a  AJk  «•  <P  A,\  *  6  A ,•*  UN«t'w  *  N„  t'u  *  N„t;.  ) 

r 

dX^y  ■ 

Cl  a;„  *  JATv  «•  tAr»/  c  N,x  t/fc  -  N,y  t'tt  *  M 

7 

d  ^  • 

Ca%  x  C  KlX  ♦  A,\  S  +  S<N  A  co!j  <i>  (.  «•  <3,  Aft  1 

: 

„  x  COSl(|)  -  ■SiNt-4u-fc*-Slw'^ 

♦  sim  CfCn  ■*  Aj* )  -  - ~  ~ 

(l-e‘  sin’’*!  ** 

}i 

d  A 

-Sin  Tv  Sin  $  CtC(x  ♦  A,\)  +  A  SI  N  $  C  *•  A,  Aft  ) 

constant  tehm 


AlfcW  f^„*Alt)  -  CCR,,  ♦A,,  )] 
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6.  EQUATION  6A  -  HORIZONTAL  GROUND  POINT 

It  Involves  two  camera  stations  and  the  corresponding  two 
Images  of  a  ground  point  whose  latitude  4  *s  known. 

It  expresses  the  condition  the  the  point  of  Intersection  of 
the  ray  A,  ,  with  the  common  perpendicular  to  rays  A,  and  A^,  lies 


on  the  cone  of  equal 

latitude  4  • 

Auxllarv  Deflations 

,  r  ae-'SiN^ 

'J  4 

V  *.  ^(«,K  *  ) 

=  1C  R,x  ♦  a,  a;x  ) 

£  ~  R  |  y  ♦  A  iy  ) 

-  2(  R(v  f  a,  A,%  1 

,  ,  _  _  ^  C.  R.,*  ♦  Alt 

1-t )  n  _  cRi,  ♦  a,  a*.  ♦  t ) 

tan1-  4 

tan1  4 

a  =  f  a»  ‘  f  a;y 

+  u>  A*# 

Cocf f Iclents 

dT„  •.  -  a,A,V  1  *  Ul)CA\  G*  <  d,A,y  1] 

dTlv  :  ^  *  CA‘X  G  v  *  d,  A,;  )  ♦  a;,  Gy  ♦  us  t  A,*,  Gy  -  <a.A,*x  )  ] 

dTlt  •.  [  VC  A\  Gj  -  d.A.v  )  *  i  CA,*y  G«  «-  d,  A*  I  4  usA\  G,7 

<JR(X  C  c  w  a*,  ♦  &  A,*y  +  u.  a;,  >  Lk  -  V  ] 

dRiy  •  A  ,V  ♦  ^  A*,  ♦  W  Af,  ILy  +  f  ] 

dft,.  (.  CW  A\  ♦  i  A‘  *UI  «•«] 
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dr^  : 

CC  *  A?*  <■  *  A,\ 

*  W  A  **  /  P«,  ] 

d  T«  ’* 

U*a;,  +  J'a;, 

♦  ui  a;,  >  "Py  1 

4T„  • 

U^a;*  -  4a,% 

+  os  A,V  }  'P*  1 

d  Rt*  • 

Uva;,  ♦  *a:„ 

os  A  r»  ^  (-  4* )  1 

d  R*,y 

L  CvAJV  +  4  A,\ 

*  ui  A^4  \  C-  4>< )  ] 

dRte  : 

U  v  a;„  *  t  a,\ 

os  Ar,  )  C-4a)  7 

JX(K  : 

a,  <.♦  t„  ♦  Jtu 

♦  wtj,)*  -fl.  (.M,xti,  ♦*  MlY  «-  M  ) 

d  tlY 

a.c*  ♦  «rt„ 

*■  <-»>  t'M.J  ♦  n.  (_  tAiy  tlu  f  M  t  lt 

;  JlCN„  t,',  ♦  Nrtti  ♦  M1#  t',,  » 

dXw  :  C  Ni*  t,'t  -  N,v  ■t'ut  *  1 

(Se1-  CoS 

d  <fc  :  -—  tan  £  (.  (  ♦  tan  *■  fi  )  ♦  — - — — - 

*•  Cl-C*-  SIN*-  <£)  '** 

cR,«  *  A|*  *  1 1  -  -  CRn  ♦  A,v)1 


CONSTANT  TURM 


TAN  1 
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7.  EQUATION  7  -  VERTICAL  GROUND  CONTROL  POIHT 

Two  camera  stations  and  the  corresponding  Images  of  a  ground 
point  are  involved,  where  the  elevation h  Is  known. 

The  condition  that  the  point  W  ,  Intersection  of  ray  A, 
with  the  common  perpendicular  to  rays  A,  and  A,  ,  lies  on  a  sphere 
centered  on  the  Z  geoenti'ic  axes  and  having  the  elevation  h  along 
the  circle  of  latitude  4  •  If  the  latitude  4  of  the  ground 
point  Is  not  given  It  Is  taken  as  that  of  point  H  . 


Auxiliary  Definitions 
d 


A 


Slrt  *-4 


a 


*•*  ♦  A i r. 

m  h 


e 


a*-  h 


^i*  * 

a*  h 

<a«t 1  si  n  4 
("  i-  eT*  sim  11  4 


Coefficients 


dr,,  :  [C)/A.;<1aI  ♦  ©cA,t,  G*  -  d,A,\)  •  ctA.VG*  ♦«.<»!] 

d  T,y  :  [  ■?(.  A,V  *  «,A,\  )  ♦  ©  A^G,,  «■  <  CA\GX  -  <3,  A»  )] 

dT,t  {,4  -<S,A*  1  *  6<A\  Gt  -d.A*.  )  ♦  <A,\G,] 
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dR,n  ■• 

[Cf  ^r»  *  «  K  +<A\  \  Lx  *  J] 

; 

( Cf  a;x  ♦  a  a:n  ♦  <  a;,  )  l>,  ♦  e  ] 

d  f?  (l  '. 

[,  t  ”?  A**  ♦■  A  AJ|  *■  4  A*|  )  L  j_  »•  i  j 

dT^  : 

[Oa;»  *  6a:v  *  <  a,\  )  t»H  ] 

dT„ 

[C?  A*  +  0A?,  -  <  A,\  )1\  J 

d  T\*  ! 

COAST4  6  a;,  «■  <  A,v  )  ?«  ] 

d  RZ)V 

[(U;  *6A,\  «■  A  A,*c)  (-U)J 

dKa.x 

UlA,;  ♦  ® A|*  *  1  A,;  K-Lyil 

dRta 

[OA.V  .  6Ar,  *  <  A,*,  )  C-L,)l 

dl„  : 

a,  tdtt,  ♦  <rt„)  ♦  OiA\*8a%  ♦  <a;,> 

tM|Xt„  »  M„  t„  *  t„) 

dI.,M 

d(  C>\  ♦  6  t|,i  •*  ^  t  ji )  ♦  ^  •'l.  A,\  «■  ©  A%  t  £A,*  ) 

*d|Y '•ll.  "*  ^i»  ^-JlT 

dll*  J 

<■1  A,;  ♦  AA*  ♦  l  A\)  C  K,^;,  .  -  Nlt  tj,  ) 

d  XM 

C  A,*  *■  6  A^  *■  <  A\  )  (.  N  lx  t,u  ♦  N|^  tjj,  *•  N  k  t,it  ) 

dk 

-  1 

tCR***' A,*)  *•('?,, +CRl#*All’)  J 
A  ♦  h 


CONSTANT  T£RM  = 


4  { 
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8.  EQUATION  56  -  HORIZONTAL  GROUND  POINT  APPEARING  ON  ONE  PHOTOGRAPH 
One  camera  station  and  the  corresponding  Image  of  a  ground  point 
are  Involved  where  the  longitude  and  the  latitude  $  of  the  ground 
point  are  known. 

The  condition  that  the  ray  A  Intersects  the  normal  N  to  the 
geold  at  the  point  of  latitude  £  and  longltu  e  >s  Is  expressed. 

Contrary  to  the  other  equations  the  derivation  of  al 1  the 
coefficients  of  equation  56  Is  to  be  found  In  (  ). 

Derived  Quantities  and  Auxilary  Definitions 

Some  of  the  vectors  used  can  best  be  defined  by  the  sketch  shown 
below: 


.0  =  exposure  station 

N  “  Vector  normal  to  the  geoid  at  the  point  of  latitude  4>  and  long- 
1  tude  7s  so  that 

*  cos  4  cos  >- 
Ny  *  cos^  sin  ■> 
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H  =  Arbitrary  point  on  R  taken  usually  at  an  altitude  g  equal  to 
the  estimated  lowest  altitude  of  the  photographed  area, 
g  can  be  taken  as  zero. 

P  =  Plane  perpendicular  to  N  passing  through  H 
h'=  Orthogonal  projection  of  0  on  P 
GX,  GY,  GZ  =  Geocentric  axes 
I  =  UH 


«  n; 

*  q  N\ 

-  . - 

Mftt1  5>n‘'  $ 

* 

* 

L  N 

4  i-e1  4 

L  -  <S  Cl  ~  3 

* 

'J  1  - « •  mm*-  4> 

R  «  GO 

E  n  OH  u  X  -  R 
k  m  mT'  ^  - 1  +  (E,ir)  n* 

S  -  Fo  -  -  (E,  N*)  N’ 

J  *»  GH'  0  R  -  S 
Ct  -  A*  x  K* 

— 

-  -JT.(K#x  E) 

V,  =  “.(E  XV) 

w  .(—-■]  Cl  t £;  N-  )  N- 

'  K,  '  IX 

ct  •  n*  )  c;  *  (c;  n*  1 1  J 


Page  99 


r 

c*  ' 

fi' 
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c 

<* 

A 

t 

* 

e 


=  COS  4  Cos  > 

=  -  5tN  4 

=  -  CjCJS  4  StN  >* 


-  Sl*J  ^  S  I  M  7> 
COS  4 


<SU-e1)^' 

(.i -e*  sin*-  4  ) 3/l- 


4-  Cj  oC' 


a  i ' 

vTcT- e.-*-  sin1- 4 ) 


v 


■ati-cM  i  ' 
Ci-a*  sins*-  4 \ 


<•  q  l' 


a  <J  _____ 

•J  (.  |-  e’  Slr4’  4‘) 

dCI-a1-)  t' 

cf-e1,  S(hiv  4) v*- 


•f* 


'J 


£' 


Cocf f  Iclents 


arx 

v*.  c  *„* 

x  I’) 

<1  Ty 

v*  ca; 

«  2*  ■> 

dT4 

v*  c  ak* 

x  Ci 

'A 

dR  y 
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:  >/» 

diK  :  vfc  cc;xt()  «■  c.;v  tu  «•  ci«  t„; 

dXv  ;  vfc  C  t  ♦  c£,  t„_  «■  C*.  t R 1 
d  4  ;  -«V/(-ll  Wy-  ew,  t*'w;  *rvj'v4t'w 

d*  ;  -  /9  w*  -  i  ^  ♦  (9'v/;  t  «f'Vy 

CONSTANT  TERM  =  -(_£.■  C.e*  ) 
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EQUATION  70  -  EQUI-ELEVATION  GROUND  POINT 


Two  exposure  stations  and  the  corresponding  Images  of  two 
ground  points  having  the  same  elevations  are  Involved,  each 
ground  point  having  an  Imaqe  1 n  each  photoaraoh. 

Superscripts  (l)  and  (2)  will  refer  to  the  1st  and  2nd  ground 
points  respectively.  It  expresses  the  condition  that  the  point 
Intersection  of  ray  A,"'  with  the  common  perpendicular  to  rays  A,0' 
and  C  has  an  elevation  equal  to  point  defined  similarly  to 
but  for  the  second  ground  point. 


e.  2.  (,  An  ••  A  ft  ) 


tT*  -  Z.  (  ft  ,x  *  A  ) 

\'f  -  Z  ♦  A,\  S 

Vjj  - - -  C  ^ie  +  1 

V*  =■  -HR1X‘A^1 

=  -l  LVZiy  -  A  ft  ) 

V,  *  -  J_  (.«,,*  *  Aft  ) 

i-e*- 
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dT;K  :  CA7V  <i>,- *  <A»  6*  l''  ’ 

■h  [  A*x£»kv*  ♦<■*?*  Gx.-A|«)VY  <-<Ar«  Or  *-A«,  } 

dT,y  1  £  [(A?*  6y  +  A„)  VK  ♦  Aiy  Gy  +  (.  Af,  Gy  -  A,*}  Vm  V 

*  (,^Ai\  Gy  ♦  A,,)  v*  ■*  A,\  6y  Vy  ■♦•  (.A(*»  Cy  -  A  ,*)  "Vg,"]  | 

d r(<  •.  (C  ia^g*  -  a.yI't;  •*•  CAt*  g»  *  a,J  v,  *  A*ma«  v„  )( 

w  ^ 

«■  [CA^G,  -AlVlvK  Y-  CAfy  C,  *A„JVy  +  A,*, o, .  v„  1  ( 

[i  4"*  ♦  A*y  Vy  ■*  A^V^  ^  *-»  ] 

*  U  a:*  *»  -  Aryvy  *  a,v v, )  u  i w  t  k„  1 

dR.y  •  {[  CA^Y*  t  A,V  Tv  ♦  Ar,  r.)  Ly]0' 

*  UA.V  V*  ♦  A*  Vy  t  Ar.  v*  )  Ly  1tW  1*  Ky  \ 

dR„  ••  { [CA ;,v,  V  A.-yTy  »-  Ar,va  ^  LBr’: 

*C.<Ayin  4.  A*y  Vy  ♦  A,m*  vB>i-BlW  *  k2  i 

dT^,  :  {  C  C  K  V*  *  A,*y  V,  f  A\  vr*  1^1° 

*  t (  Afx  ♦  A \y  Vy  *  A,’t  VjJ'Ph  ]*■  ] 

drty  :  {[CA^V*  *  A?,  Vy  ♦  A,;  V^'Pyf' 

f  [  <  V*  +  A,*y  Vy  ♦•  A\  V„1  'Pyl  **  ^ 

iTtg_  Af*  V*  ♦  A'yV^  f  A,,  Til'P^l 

*  t  CA,V  v*  «•  Ary  Vy  «■  AT*  V,  w  } 
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dftu( 

tUA,;v»  - 

*  Yy  ■* 

■  A\r-m)  (-LOl^ 

♦  UA%  v* 

*  A  *y  Vy  i 

•  a;,  v€  j  c-lO  1 (W  \ 

d^tt  ’ 

L  (  at*  y^ 

*  A*  Yy 

♦  A,V  VO  C-U)^" 

•*  CCA.V  v* 

+  A*y  Vy 

*  a,v  v«)  c-LYn*'  1 

(ua:*  v„ 

4-  Aiy  Yy 

4-  A,V  Y^l  C-Ljj)  ] 

*UA*'r* 

♦  A*y  vr 

-  A,*,  vx*  )  C-L,)]*  ] 

■ 

a1:1  cyn  t„ 

■*  Y,  t ,,  1 U  ♦  C-Y*  A,V  +  Yy  A,y 

'T.a;,)1' 

<Mlxt„  + 

MlV  ti<  ♦ 

*«*•*> w 

dX,,1"  ’• 

a?  lrK  tlfc  «■  V,  +  <-y*a:x  *x,  a 

*y  4  Y»  A*,  } 

^«t,u  * 

lAiy  ■*•».». 

^  M.4  t„J1' 

drj’  : 

<■  y,  a:h  - 

Yy  A,\  ♦  Yt  Aft,»UCNmte  +N„tw  ■♦ 

N.»  *j.iW 

dX^'  : 

it*a,v  4- 

Yy  A*y  4 

T,  A\l'[>  <.H,»  tIL  «■  M,y  ttu  «• 

N,» 

<Jt,yW  • 

d,  W,  t  ,| 

1  ■*■  Vy  -t  u 

>  Vx  'ty,  )  +  C^v  AT,  4  '/'y  A  1*1 

♦  v/»  a,;  ) (v 

W 

tMwt„  ♦  MlYt  u  ♦ 


dtlvW  d,^  (.  Vx  tlU  «•  Jy  ta  ♦  \/t  CVK  Atx  +  S  AIV  *  ''ji  A,\\  ^ 

Miy'ttj.  *  ^i*  fc-vi.' 

dx^  ■•  t  v„  a;,  +  s  a;,  w,  Aft) (ll)  CN,*t„  *■  KJnrt*  f  t„)  w 

dlu1"  ;  tv*  A,\  *  V  y  Ar,  +  v/*  A;*)^  C.N^-t.u’-  N,y  - 
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Constant  term  -  -  [R1kkx  ^  (A,k  <  A,?'1  ^ 


(n  ->*. 


K,^  +  <-  A^y  )  *"  "  (■  A|y  ) 


W  ,  »- 


‘I'f 


*  *„  It,  WA^^-CA,®^] 
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1 0 •  EQUATIONS  lOj.,  m?,  10?  -  BLOCK  ADJUSTMENT 

Four  exposure  stations  and  the  corresponding  Images  of  a 
ground  point  on  each  photograph  are  Involved.  Considered  simul¬ 
taneously  they  express  the  condition  of  Intersection  at  one  point 
of  the  four  rays  A^A^A*  ,  Aj  from  the  four  perspective 
centers  to  the  corresponding  four  Images. 


a.  Eauatlon  1Q1  Coefficients 

ar,K  : 

Arx  Cix 

dT,v  : 

A‘x  Gy  ♦  O,  A(il 

dT,z 

A;x  <1*  -a,  Aft 

«.*  : 

AT,  u**  ' 

dR,y  : 

A  lx  i-y 

d  (4 1  e  ! 

a*h  l£ 

a  r„  • 

A'y.  Py 

artv  : 

A|*h  Ay 

dTy,i  Am 

d I?  t*  •  -  A,*  L  n 

JdtY  ;  ■  A,x  Ly 

d  -  AJx  L  ( 


W^W.ftJSIfUrtWW  WlMJW  I4MIW 
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dX„«  :  -  a.4**.  -  At,*  (M»*  +  Nl*  *  Mt«  U*  ) 

di*  -  a;**  in*  -t/„4  ♦  n,;  i'u*  ^  n,*  t;*  > 
dXV  :  -  A,;4  CN*  4  Nt*  t'u*  4  Nt;  t'lt*  ) 

CONSTANT  TE.fi  M  "•  <*tn  4  A,*  -A-i*  -  A,„ 


b.  Equation  10  2  Coefflc Lents 


<*rIA  : 

Aim  0»  '  d,  A, 

dT,,  • 

A,y  Oy 

dT,„  •• 

a;m  c4  ,  a,  A,v 

A*  y  L  H 

dR„  •• 

-  <  A,;  Ly  4  .  ) 

d  A  is  I 

A,\  L  * 

dT«.*  « 

Aim  “P, 

dT&y  : 

drlt  : 

A  nr  "P* 

dftjj; 

-  A*y  L* 
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-A;tW 

d  dL1i 

-  A?,  L* 

-  C  A,*  *  £>  x  '  ^i*  A,V  *  l 

ATj  : 

-  A?y  *  G  v 

dr,;  ' 

-CAr/G**  ♦d*A,\*i 

d«.J  ■■ 

-  AV  L/ 

dftiy 

••  C  A  [y  *  Ly  ♦  1  1 

dR,{ 

.  A*  *  L  * 
u* 

dT.;  : 

.  A  0  •  'p* 

“  rx 

drt; 

A*  *  V  * 

“  Mir  K,< 

dTt;  ■ 

-  aW 

dRj  ' 

A„  *  U* 

dRty 

A  •  ♦  1  * 

“  i  Y  *-  y 

dR«  : 

A*  *  L  * 

A(y  L| 

dl'ix  : 

«S,t«,  ♦  A,*t 

<  M„ t „  ^ 

H,y  t  u 

* 

A  L,y 

<9,  til  «■  A*y 

Mi* 

v  M,» 

A|  y  ^  N ,  „  t  „ 

* 

*  rs„ 
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d  Xty 

A?y  INut,;  ♦  N,y  tu  + 

taj 

d  Xix 

-  a, *  't  2^  -  A,y*  C  (M,x  t.  „* 

*■  M(y  t  t?  ♦  M(t  t  jt  ) 

ax*  : 

-a;  t  **  -  A,‘y*  £  M.J  t  ,i  * 

M|<<  ^*#L  *  ^31.^ 

diti  : 

-  a:/  (  n,<  t'„*  N(*  t;/ 

*  Nlt-  t;-> 

d  1  ty 

-  av  («,:  t;t4  *  n,/  tU* 

r  n(;  i 

CONSTANT-  T'EHH  -  R.*  *  A,y4  -  RIT  -- A,t 


dTix  A,  ,  Ok  ♦  <3,  A,\ 

d  T,y  A,t  G  f  -  d,  A,* 

dr„  1  A,;  Gi 

dRlx  a;.  L* 

d  P^iy  A,t  (~y 

dR,«  ■  A\  L,  ♦' 

dTu  A,\  VA 

4Ttt  A,a  'Py 
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dTt<  A  ‘t  'Pjj 

:  _  Af*  L* 

dR„  :  -  A;t  Ly 

~  A  ,,  L  j 


*r*  : 

-  CAr,  *  Gk  -  4"  A,\ 

dT,;  : 

-  (A%*  C.Y‘  -  a,*  A,*/ 

dr„*  : 

-  Ar/c,/ 

dfc,: 

•  AT*4  l/ 

dR/y  •■ 

-  a,*/  l  ; 

dR/a4 

-  CAT,4  L4  *i  j 

dTw 

-  a;,4  V* 

d  I'm4 

-  A,*,  4  9* 

tr*;  ■■  -  a;,*  ?Y 

d^M4  :  a(;4l; 

4*t?  :  A,*,*  Ly* 

d(R  t*  A,|*  L* 

In 


<3,  tv  +  A,;  C  M„t„  *  M,y  tu  *  t,J 
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dXiy  :  3,  tst  A|i  C  M,x  1 1*,  t- 


a,; 

CM« 

t;  - 

Niv’tt.i  *■  A1  n 

t;,. 

,1 

dX«  ■■ 

A,% 

<■  N«t;A  - 

N,t  tu  '  N- 

.  * 

•JT.-) 

ax,;  : 

-  a,* 

tv 

-  A**  4 

1  c  m*  t,,*  *  m.; 

*  m,; 

tl.‘> 

ax,;  • 

-  a* 

-  A  “  1 

*  c.  m,;  t.4*  * 

tat 

«■  m.; 

t,t) 

diii  : 

"  A|  j 

,*  c 

m.x 

«  *.  M  *  + ' « 
fNiy  '-t.i 

<«■ 

M  * 

t'M 

ax«  ■ 

-  A,  j  4  ( 

<  t 

/  *  „  M  *  4-  '  •  . 

It  *  N<t  -tt 

►  n,; 

*. >  * 
-  IX 

) 

C.ONSTAMT  TERM  =  R,/  *  A,,4  -R,,  "A,, 
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11 .  EQUATIONS  8.  Q.  10  -  EXPOSURE  STATION  OF  KNOWN  ♦  .tn  and  A. 

They  merely  express  dR  ,  dR  ,  and  dR  In  terms  of  the  unknown 

x  y  * 

residuals  d$,  d*.  and  dh. 
a.  Equation  8  Coefficients 

dTK 
d7\ 
dra 

d  ft, 
dRt 
d  4 
d  7k 
d  h 

constant  term  =  o 

b.  Enuation  9  Coefficients 
dT*  :  O 

dTy 


-Vx 


[caKt-e,'lCl-stSiiMt4i)  *  h  ]  E>in4>C£>S7\ 
[  3(1  -  e*  2»iM  *-4  )  L4  h  ]  C£S  4>  SIN  Ti 


-  C.OS  4  ccS  A 


o 


i  jl  I  ILL1,  d  b  M I  k%l  I  ‘■W  j/WiM  MitfbrViWj  f#J  Wii.  SU^iS*  (H»  rt  f)»n«P 
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dTa 

:  o 

a 

*> 

* 

Q 

dfcv 

:  l 

dRa 

c> 

d  <P 

,  -VL 

[  a  C.  i-  e1-  )(.  l  -  e*  Sin*  $  ) 

In  ]  5 »N  $  SloJ  y 

d  A 

T  -{L 

-  [  <s»  c<  )  ►  h  J 

COS  <t>  Cos  A 

dll 

-  COS  4  &<N  > 

CON 

ST  ANT  TE.H.M  =  Q 

C  . 

Equation  10  Coefficients 

6Tk 

O 

d  Tv 

0 

<*r» 

o 

dR  h, 

0 

dRv, 

o 

*** 

1 

d  <$> 

-  [  <3  (  1  -«*■)  Cl  -«**  SIN1-  <fc)  ^ 

*  h  ]  COS 
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d  7s  0 

dW  -  S  in  4 

CONSTANT  TE.RM  =  O 
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12.  EQUATION  1)  -  EXPOSURE  STATION  OP  KNOWN  A 

The  perspective  center  Is  forced  to  lie  In  the  plane  of 
equal  longitude7'. 

Coef f Icients 


drx  : 

o 

ar, 

O 

<j  ra 

o 

dR* 

S(N  A 

•d 

-  CCS  A 

dRa 

Q 

d  A 

OOT,  A  G, 

CONSTANT  TER.NI  =  CO^^R, 


5iN  A  R* 
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1 3 .  EQUAT I  ON  12  -  EXPOSURE  STATION  OF  KNOWN  LATITUDE  <6  AND  LONGITUDE  a 

The  perspective  center  Is  required  to  lie  In  the  prime  vertical 
plane. 

Auxl I iarv  Def 1 nl t Ions 

X  c  coa  A  sin  4 

S  =  SIN  X  S!N$ 

£  =  -  COS  4 

^  da  StN  4  cos  4 

O  =  1 TZZZ  -  1 . Z  - - 

Tt-e*-  5  IN  ‘  4 


Coef f  icients 

arx 

o 

cJT,  • 

O 

*!• 

-a 

0 

4R* 

a 

i 

<JRa  : 

e 

d  X 

-  $  ft*  t  a  Ry 

d  4 

CoS  X  COS  4  *■  St*  ^  cos  4  *  S.IN4  Rt 

t  Cos4  4-  SIN4  4  U- «'  Sin’'4) 

-  da  - - - - - - 

c  l  -a. 1  St  N 1|"  $  ) 


r  C«,  ) 


Constant  term 


3  -  [SRxrfR, 
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1 4 .  EQUATION  12A  -  EXPOSURE  STATION  OP  KNOWN  LATITUDE  cfr 

The  perspective  center  Is  forced  to  lie  on  the  cone  of  equal 
latitude  4> • 

Auxiliary  Definitions 


<3«*"  Slid  4 


ri-e"*  sim*“  4 


Coef  f ! c ' 

lents 

dV* 

0> 

dt^ 

0 

dT, 

o 

dR, 

2Rk 

<iU.y 

2.R., 

-  ;.(.R4  t) 

TAN  4  4 

6<f> 

UR,  ♦  t)  r  R, 

TAN  4  4 


Sim  4  cos  4 


aetco*  4 

U-c.4  mnMT*'4 


<•*«  ►  t  )  *“ 


constant  tsj*.n\ 


tan4-  4 


L  *  *"i  ] 
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15.  EQUATION  M  -  EXPOSURE  STATION  OF  KNOWN  ELEVATION 

It  enforces  the  perspective  center  to  He  on  a  sphere 
defined  similarly  to  Equation  ?• 

AuxIHarv  Definitions 


ca 

n  -  - 

a.*-  SlM1,  $ 

ae'Sifv)  & 

t  =  - . -L— 

'fl-e’-sini  *■ 

Coefficients 


dr, 

d  "Ty 

dTE 


o 

o 

o 

fin 

t  «  flrt) 


dR^ 


(  A  ^  h  ) 


Re  *t 
(.n  *  n  j 


dh  -  l 


_  C  R,  R>,  <-  CAf 


Co  Mutant  term 


=  {_  [  vA  t  In 


t.  n  •*  v» ) 


♦t)1] 
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VIII.  OPERATING  INSTRUCTIONS 


A.  Loading  Seouence 

The  program,  control  cards,  and  input  data  are  written  on  mag¬ 
netic  tape  in  the  following  sequence: 

1.  ID  Control  Card  (monitor  accounting) 

2.  XEQ,  Control  Card 

3.  Chain  I  Control  Card 

4.  Segment  I  Binary  Deck 

5 .  Chain  2  Control  Cord 

6.  Segment  2  Binary  Deck 

7.  Chain  23  Control  Card 

8.  Segment  73  Binary  Deck 

9.  Chain  3  Control  Cord 

10.  Segment  3  Binary  Deck 

1 1 .  Chain  4  Control  Card 

12.  Segment  4  Binary  Dock 

13-  Data  Control  Card 

14.  Data  Deck  -  Binary  Coded  Decimal 

a.  J  -  Job  Card 

b.  P  -  Parameter  Card 

c.  E  -  Exposure  Stations 

d.  G  -  Ground  Control  Points 

e.  S  -  Stop  Card 


B. 


Page  120 


Magnetic  Tape  Usage 

In  addition  to  the  tape  units  used  by  the  Fortran  Monitor 


System,  the  following  tapes  are  used  for  temporary  storage  by  the 
object  program. 

A4  -  Segments  1  and  2 

B1  -  Condition  Equation  Sets 

B2  -  Segments  23,  3,  and  4 

B3  -  Exposure  Station  unit  variance,  ground  point  data, 

VECDEX,  MATDEX 

All  output  Is  put  on  tape  A3  for  off-line  printing  and  punching. 

When  the  run  Is  completed  either  by  norma)  exit  or  after  error 
exit,  control  Is  returned  to  FMS. 

C.  Innut  Format  Guide 

1 .  Summary 

Figure  I  displays  the  eight  types  of  "IBM"  cards  on  which 
the  Input  data  for  a  block  or  strip  tr langul at  Ion  must  be  punched.  The 
card  are  then  arranged  in  a  dock,  as  follows: 

Job  Card 

Parameter  Card 

All  Exposure  Station  Card  sets 

All  Ground  Point  Card  sets 

Stop  Card 

The  following  paragraphs  discuss  the  contents  of  each 
card  type  in  detail,  referring  to  successive  lines  of  Figure  I. 

2.  Job  Card  (Line  l) 

The  computer  Identifies  the  job  card  by  tho  letter  J 
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In  column  1.  The  numbers  of  the  series  and  test,  type  of  photography,  and 
number  of  flights  are  read  and  stored  by  the  computer.  However,  they 
have  no  significance  for  tho  actual  tr langulat Ion  process,  and  may  be 
left  blank. 

N  Is  the  number  of  exposure  stations  or  photographs.  After 
all  cards  have  been  read,  N  and  the  number  of  ground  points  (appearing  In 
columns  29-35)  will  be  compared  with  the  actual  tallies,  and  an  error  stop 
will  occur  If  too  few  or  too  many  exposure  stations  or  ground  point  cards 
wore  read. 

The  unit  standard  deviation  ot  the  if>  and  \  for  exposure 
stations  and  ground  points  Is  defined  in  column  37-39  In  thousandths 
of  seconds  of  arc. 

Tho  datum  g  gives  the  altitude  In  feet,  relative  to  the 
given  spheroid,  of  the  surface  from  which  all  altitudes  arc  given.  If 
columns  hO-50  are  left  blank,  the  program  will  assume  that  all  altitudes 
In  tho  Input  data  are  corrected  to  the  reference  spheroid. 

Tho  focal  length  f  will  be  taken  os  the  "z"  coordinate, 

In  the  photographic  model,  corresponding  to  t he  measured  values  of  x  and 
y  which  give  the  position  of  a  ground  point  on  a  photograph,  lnus,  x,  y, 
and  f  must  hove  the  same  dimension,  e.g.  all  inches  or  all  centimeters. 

The  negative  sign  of  f  Is  necessary  for  ce.  rect  orlerMtlon. 

The  unit  standard  deviation  of  the  altitude  (h)  used  for 
the  exposure  stations  and  ground  points  Is  specified  In  columns  63-68 
in  thousandths  of  feet. 

The  unit  standard  deviation  to  be  used  for  the  photographic 
x  and  y  coordinates  is  defined  In  columns  70-72  In  microns. 
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Cl umns  73-80  of  all  eords  are  not  read  by  the  computer, 
and  the  divisions  shown  are  only  a  suggestion.  Another  possibility  Is  to  use 
them  to  give  consecutive  numbers  to  every  card  of  a  deck,  so  that  Its 
original  order  can  be  restored  In  case  of  accident. 
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J  Card  Format 


In 

-3 

sec  x  10 

Col. 

37 

39 

(4> .  * ) 

ft  x  10“ 3 

63 

68 

(h) 

microns 

70 

72 

(x,y) 

F.  Card  Format 

/r  In 

sec 

Col . 

i 

12 

(<M 

rf 

sec  x  10 

15 

17 

(h 

n' 

ft 

18 

24 

(h) 

a' 

ft  x  IQ-3 

26 

28 

<h) 

<f 

sec 

29 

35 

(M 

see  x  10  3 

37 

39 

(M 

G  Card  Format 

First  Card 

1  in 

sec  x  10  3 

Col . 

26 

28 

(<f>) 

' i 

ft.  x  !0“3 

29 

35 

(h) 

sec  x  10  3 

37 

39 

U) 

R  Care!  Format 

anyth Incj  in 

microns 

Col . 

2 

6 

'1  . 

microns 

7 

13 

1  X 

lV  iy 

n  „ 

microns 

15 

17 

microns 

18 

24 

2x 

26 

28 

'f  2y 
''  3* 

''  3y 
•'  lx* 

'  »v* 

*  i  v< 

2x 

microns 

microns 

29 

35 

microns 

37 

39 

microns 

40 

46 

nlcrons 

48 

50 

microns 

51 

58 

G  * 

*y 

<f _  * 

3x 

r .  * 

microns 

60 

62 

microns 

63 

68 

microns 

70 

72 

3y 

PHYSICAL  UNITS  OF  INPUT  DATA 
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3.  Parameter  Card  (line  2) 
Format 


Column  1 : 
Columns  3*12; 

Columns  14-23! 
Columns  25-32: 

Column  34: 

Columns  36-39: 
Columns  41-43: 
Columns  45-55: 


Columns  57-62: 


Columns  64-69: 


P 

Integral  part  of  1024a,  where  a  -  semi-major  axis 
of  reference  ellipsoid  In  feet  (standard  value  of 
a  «  20,925,832.16) 

Integral  part  of  1024b,  seml-inajor  axis  (standard 

value  of  b  =  20,854,892.02) 

35  2  2  2  2 

Integral  part  of  2  c  ,  where  e  a  I  -  b  /a 

(2^  «  34,359,738,368)  (standard  value  of 

»  .006,768,658,0) 

number  of  Iterations  before  resection  and  break¬ 
point  stop  (pressing  the  START  button  will  Initiate 
another  set  of  Iterations)  (standard  value  4) 
uslnj  3 

number  of  allowable  cycles  of  Simultaneous  Solution, 
assuming  continued  but  slow  convergence, 
cycles  of  Stlofo)  method,  In  each  batch,  hetv.'een 
convergence  tests  (standard  value  6  x  no,  of  photos) 
convergence  criterion  f,  where  Simultaneous  Solution 
Is  assumed  correct  If  for  every  normal  equation  con* 
stant  term  k(,  the  corresponding  rojlduol  r^  satisfies 
r  j  |  <1  f  |  kj  |  ;  f  must  be  greater  than  .000,000,05 
In  any  case  (decimal  point  and  up  to  10  digits) 
(standard  value  .000,005) 

scale  factor  for  dT  coefficients;  coefficients  are 
divided  by  this  factor  before  normalization  of  the 
condition  equation,  so  that  the  dT  values  from  the 
Simultaneous  Solution  must  be  divided  by  It  also 
(division  occurs  before  printing)  (factor  Is  an 
Integer)  (standard  value  l) 

ratio  c  for  computing  artificial  lower  bound  for 
Chebyshev  portion  of  Simultaneous  Solution  (upper 
bound  b  Is  computed  from  normal  matrix,  Chebyshev 
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method  is  then  set  for  eigenvalues  In  the  range 
cb  to  b)  (decimal  point  and  up  to  5  digits) 

(standard  value  .02) 

Columns  71 "72:  number  of  cycles  n  of  Chebyshev  methed;  for  a 

c 

normal  accuracy  a, 

.  .  cosh" *  a 

n&  ■  greatest  Integer  in  — 1  .  - 

cosh"  l+c 
I -c 

(standard  value  50) 

General 

Any  digits  punched  in  the  columns  listed  will  be  taken  as 
new  values,  lo  be  substituted  for  the  standard  cases.  Notice  that  all  values 
except  the  convergence  criterion  f  and  the  ratio  c  are  pure  integers  with 
no  sign  or  decimal  point,  while  f  and  c  each  contain  a  decimal  point  and 
may  (but  usually  will  not)  have  digits  or  zeros  to  the  loft  of  the  decimal 
point.  No  spaces  should  bo  left  between  digits. 

4.  Exposure  Station  Cards  (lines  3  and  4) 

A  pair  of  cards,  each  with  the  letter  E  In  column  1, 
must  be  made  up  for  each  exposure  station.  Columns  2  through  5  of  the 
first  card  (and  the  second,  if  desired)  contain  the  number  of  the  photo¬ 
graphs.  Photographs  must  be  numbered  consecut tvely ,  and  In  the  same 
direction  on  each  strip  of  a  block.  The  number  of  the  first  photograph 
may  be  I  or  larger,  but  not  zero. 

Column  6  contains  the  exposure  station  type  letter  as 
taken  from  F igure  3 • 

Columns  7-13  contain  the  Integer  portion  of  the  standard 
deviation  used  for  the  exposure  station  latitude  in  seconds  of  arc- 
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Columns  15-17  contain  the  fractional  portion  of  the  standard 
deviation  used  for  the  exposure  station  latitude  c|»  In  thousandths  of 
seconds  of  arc. 

Columns  l8-2<f  define  the  integer  portion  of  the  standard 
deviation  used  for  the  exposure  station  altitude  (h)  in  feet. 

Columns  26-28  define  the  fractional  portion  of  the  standard 
deviation  used  for  the  exposure  station  altitude  (h)  in  thousandths  of 
feet . 

Columns  29-35  contain  the  integer  portion  of  the  standard 
deviation  for  the  exposure  station  longitude  A  given  In  seconds  of  arc. 

Columns  37-39  define  the  fractional  portion  of  the  standard 
deviation  tor  the  exposure  station  long  I tude r)  given  In  thousandths  of 
seconds  of  arc. 

Columns  ho-50,  51-62,  and  63-72  contain  such  data  as  Is 
consistent  with  the  typo;  If  a  value  Is  not  known.  It  Is  left  blank.  Note 
that  latitude,  <|’  ,  ranges  from  -90H00I00.000"  to  t90°00l00.000",  whllo 
longitude  ranges  from  - I80°00' 00, 000"  (West)  to  +180°00,00.000"  (East). 
Zeros  for  minutes  and  seconds  must  be  punched  on  the  card;  these  are 
"trailing  zeros"  which  establish  the  degrecs-mlnutes-seconds  relationship. 
Loading  zeros,  however,  need  not  be  punched;  for  Instance,  -00°03' 10,567" 
could  have  simply  -JI0  anywhere  In  the  field  preceding  the  decimal  point. 
Minus  signs,  if  used,  must  not  be  separated  from  the  leading  digit,  nor  may 
digits  be  separated  by  blank  columns  within  a  field. 

Note  that  decimal  points  may  appear  only  In  certain  columns 
(14,  25,  36,  47,  59,  69).  If  desired,  decimal  points  may  be  gang-punched 
in  all  cards  In  these  columns,  even  where  (as  in  N  on  the  Job  Card)  the 
"decimals"  have  no  relation  to  the  "integer"  to  the  left- 
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The  reader  should  also  remember  to  always  write  three  figures 
in  the  decimal  positions  if  they  are  in  fact  part  of  a  number  to  the 
left,  since  (for  all  angles  and  elevations)  conversion  Involves  division 
of  the  number  in  these  columns,  taken  as  an  integer,  by  1000. 

The  second  card  for  each  exposure  station  contains  the 
estimated  heading,  tilt,  and  geographic  location  of  the  photograph's 
perspective  center.  None  of  this  information  may  be  omitted.  If  <1'  ,  ^  , 
and/or  h  differ  from  values  on  the  preceding  card,  the  "estimates"  will  be 
used  in  the  preparation  of  the  orientation  matrix  (M) . 

5.  Ground  Point  Cards  (Lines  5,6,7,  and  8) 

Each  ground  point  Is  represented  In  the  Input  dock  by  two 
or  three  cards  containing  a  G  In  column  1.  The  number  of  the  point  appears 
on  ttic  first  card;  It  will  bo  used  to  identify  the  point  In  the  output. 

The  order  In  which  ground  points  arn  processed,  for  Instance  In  forming 
condition  oguaV, ions ,  Is  the  same  os  their  order  In  the  Input  dock. 

Column  6  of  the  first  card  shows  the  type  of  point  as 
taken  from  F  igurc  3 • 

The  conjunction  point  number,  and  number  of  first  photo¬ 
graph  of  conjunction  point,  are  taken  (if  present)  as  a  requent  for  Condition 
Equation  70.  This  equation  is  not  otherwise  computed.  It  Involves  two 
ground  points  appearing  on  the  same  pair  of  photographs.  The  specified 
conjunction  point  must  fol low  (but  not  necessarily  be  next)  the  current 
point  In  the  Input  deck.  The  pair  of  photographs  must,  of  course,  be 
consecutive  on  a  particular  strip;  they  could,  however,  be  in  different 
locations  with  respect  to  the  ground  points,  e,g.,  as  photos  2  and  3  of 
strip  2  for  the  first  point,  and  photos  I  and  2  of  strip  1  for  the  second 
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point.  The  cards  Tor  the  specified  conjunction  point  do  not  have  any 
special  marks  of  Identification;  If  numbers  appear  In  columns  7-13  and 
1 5-24  of  the  first  card  for  that  point,  they  must  refer  to  a  second 
application  of  Equation  70  with  a  third  point  further  on  In  the  deck. 

Columns  26-28  contain  the  standard  deviation  for  the 
ground  point  latitude  given  In  thousandths  of  seconds  of  arc  ( ), 

Columns  29-35  define  the  standard  deviation  for  the  ground 
point  altitude  (h)  given  In  thousandths  of  feet. 

Columns  37-39  contain  the  standard  deviation  for  the 
ground  point  longitude  (A  )  given  In  thousandths  of  seconds  of  arc. 

INFORMATION  CODE  FOR  GROUND  POINTS  AND  EXPOSURE  STATIONS 


Known  Information  Letter  Code 

latitude,  longitude,  elevation  A 

latitude,  longitude,  elevation  D 

(not  of  equal  quality) 

latitude  and  longitude  C 

latitude  and  elevation  D 

latitude  E 

longitude  and  elevation  F 

longitude  G 

elevation  H 

none  K 

computed  In  resection  but  not  used  for  photo  L 

orientation 


Figure  3-  Type  Codes 
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Geographic  coordinates  appear  If  and  only  If  the  type 
letter  of  column  6  implied  that  some  or  all  of  them  were  known. 

The  second  card  for  a  ground  point  contains  an  Identifying 
"1"  (for  strip  1)  in  column  6.  X^  and  y^  are  the  measurements  (normally 
in  inches)  taken  from  the  photograph  whose  number  is  given  in  columns 
2-5,  X^,  y^  and  x^.  Yj  31-0  left  blank,  or  contain  measurements  of  the 
same  point  from  the  photographs  which  follow  in  the  strip,  i.e.,  If 
tho  first  photo  Is  Ho.  7.  X£  anfl  Y?  must  bo  taken  from  photo  No.  8, 
and  x, ,  Vj  from  No.  9*  If  a  point  happens  to  lie  at  (0,0)  on  the  photo- 
graph,  at  least  one  zero  must  be  punched  In  each  position  for  x  and  y. 

If  x.j  and  y^  arc  blank,  the  computer  will  not  look  for  x^  and  y^. 

Tho  sign,  digits,  and  decimal  point  In  columns  7-15,  18-24, 
etc.  may  bo  separated  by  blanks  In  tho  case  of  x-y  measurements.  Plus 
signs  need  not  appear  (this  also  applies  to  all  signed  data,  on  exposure 
station  cards  as  well). 

A  third  card,  Identified  by  the  G  In  column  I  and  a  2  in 
column  6,  appears  only  If  data  for  a  second  strip  Is  available.  Photo  I 
of  strip  2  must  then  overlap  photo  I  of  strip  I,  and  so  forth.  Unless  the 
ground  point  Is  of  type  A  or  C,  data  from  a  single  photo  In  the  second 
strip  cannot  be  used  and  should  be  left  out  ot  avoid  stalling  the  computer 
In  the  middle  of  the  Condition  Equation  former. 

Strip  I  may  lie  either  above  or  below  strip  2  for  any  ground 
point,  the  only  difference  being  that  condition  equations  5,  6,  6A,  and/or 
7  are  computed  pnly  for  the  first  two  photos  of  the  first  strip. 

Each  set  of  ground  point  cards  must  be  terminated  with 
an  R  card.  This  card  specified  the  standard  deviations  to  be  used  for  the 


Page  130 


n  and  y  coordinates  as  defined  by  the  following  format: 
Column  I  R 


11-13 

photograph  1 

strip  1 

15-17 

photograph  1 

strip  1 

22-2/+ 

photograph  2 

strip  i 

26-28 

photograph  2 

strip  1 

33-35 

photograph  3 

strip  1 

37-39 

photograph  3 

strip  1 

44-46 

photograph  1 

strip  2 

48-50 

photograph  1 

strip  2 

56-58 

photograph  2 

strip  2 

60-62 

photograph  2 

strip  2 

66-68 

photograph  3 

strip  2 

70-72 

photograph  3 

strip  2 

6.  Stop  Card  (not  shown) 

This  card  contains  on  S  In  column  I.  Any  other  data  on  the 
card  Is  ignored. 

D,  Output  format  Guide 

I .  Cycle  Output  * 

On  every  cyclo  including  the  last,  the  corrected  values  of 
(H)  and  R  dr©  printed  in  tins  following  format: 


title  line: 

Photo 

(No.) 

data  1 lnes: 

(1) 

R 

X 

M  M  M 

1 1  12  12 

(2) 

R 

y 

M21  M22  M23 

(3) 

R 

2 

M3I  M32  M33 
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The  first  series  of  this  output  shows  the  Initial  values 
prepared  from  estimates  and  contains,  In  addition,  the  type  letter, 
radian  values  of  phi  and  lamda,  and  altitude  on  the  title  line.  Matrix 
values  are  printed  In  floating  point,  as  a  signed  power  of  10  followed 
by  a  signed  mantissa  of  eight  digits.  For  instance,  the  numbers  0-0987, 

0.987  and  -98,785,1*32,000.0  are  printed  as: 

-I  98700000 
0  98700000 

and  1 1  ~S)8765432  . 

2.  Final  Output 

After  cycle  output  on  the  last  cycle,  the  heading,  swing, 
azimuth,  and  tilt  are  computed  for  each  photograph,  and  geographic 
coordinates  are  computed  for  all  exposure  stations  and  ground  points. 

Angular  data  arc  printed  as  degrees,  minutes,  and  seconds,  the  three 
numbers  being  separated  by  spaces:  for  Instance,  )h6°8' 40. 355"  Is 
printed  W»6  8<t0.355>  The  first  list  Is  for  exposure  stations: 
title  line:  PHOTO (No,) 
data  line:  H  s  a  t  t  t  d>  A  h 

*  y 

Thic  is  followed  by  the  list  of  ground  points: 

title  line:  GROUND  P0!NT(No.)  TYPE (t)  KNOWN  (‘I  )  (*)  (h) 

data  lines:  CALCULATED  ij>  A  h  Ax  1  y  z 

In  the  title  line,  following  "KNOWN",  is  the  control  data  as  given  in  the 
Input  on  the  first  card  for  this  ground  point;  If  any  values  wore  not 
given,  they  are  replaced  by  single  x's.  The  calculated  values  of  f  ,  *  , 
and  h  are  invariably  obtained  as  follows.  First,  values  of  x,  y,  and  z 
computed  for  each  pair  of  photographs  on  which  the  point  appears.  These  are 
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averaged,  and  the  average  values  used  to  compute  </,  ,  a  ,  and  h.  The 
differences  (in  absolute  value)  between  the  average  values  of  x  ,  y  ,  and 
z,  and  the  values  which  are  furthest  out  of  line  (the  maximum  deviates)  ore? 
printed  as  Ax,  Ay,  and  Az. 

3.  Matrix  Solut Ion 

On  each  cy  cle  Including  the  last,  the  solution  vector  to 
the  normal  matrix  Is  printed  In  floating  point,  Its  elements  In  the  order 

dT|x'  dT|y*  dTlz*  dKix*  dRly’  dRlz’  dT2x’  0tC*  thr0ugh  dV 
b.  Condition  Equation  Coefficients 

During  Segment  2  the  coefficients  of  each  condition  equation 
for  each  exposure  station  or  ground  point  will  be  printed  as  soon  as  they 
are  formed  and  have  boon  treated  by  the  weight  former; 

title  lino:  COEFFS.  OF  COND.  EQ.  CALLED  FROM  (nnnnn) 

PH 01  OS  (I)  (J)  (k)  (I) 

data  lines:  (Coefficients  of  dT(^,  dT(,  ,,,  ,  dRjz>  dTjx. 

« • i  ,  dR ,  ,  •••  ,  d R ,  ,  d .  ,  d ,  ,  ■ .  .  d  , 

’  |z’  lz  lx’  Iz  ex’ 

d  ,  dV>  ,  d\  ,  dh  and  constant  term,  all  in 
oy 

floating  point) 

The  octal  number  (nnnnn)  in  the  title  line  refers  to  a  line  In  the  Condition 
Equation  Former  program.  A  list  of  these  numbers  can  be  made  up  by  inspect¬ 
ing  the  list  produced  by  FMS  on  compilation  of  Segment  2,  and  noting  all 
references  to  equation  former  subroutines. 

5,  Normal  Matrix  and  Constant  Vector 

A  suitable  title  line  will  be  printed,  followed  by  the 

36(N+K)  elements  of  the  normal  matrix,  the  6N  elements  of  the  constant 

2 

vector,  and  the  sum  of  the  squares  of  the  constant  vector  (k  ).  The 
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matrix  elements  are  printed  In  36’word  blocks,  beginning  with  the  N 
ma I n-d 1  agonal  blocks  (in  the  same  order  as  the  terms  of  the  vectors  - 
see  h.  Matrix  Solution  above)  followed  by  the  K  non-vanishing  blocks 
In  the  upper  triangle  of  the  matrix,  In  order  row  by  row. 

6.  "Data  Processing"  Output 

Segment  I,  the  following  data  will  be  printed  for  each 
ground  point: 

GROUND  l*T ,  (No.)  (Type)  (0  radians)  (  radians)  (li) 

(X)  (V)  (Octal  Equation  70  code  word) 

(Octal  Strip  I  Index  word)  (a^)  (oj  )  (a^)  (a^)  ...  (n^) 

(Octal  Strip  2  Index  word)  (n^  ) . (a^  ) 

Of  this  data,  .  ,  •  ,  and  h  are  printed  only  If  given;  X ,  Y ,  and  Z  are 

computed  and  printed  only  If  i  ,  ,  and  h  are  given;  and  the  strip  2 

data  line  Is  omitted  If  there  Is  no  Strip  1  Input. 

A t  completion  ol  Segment  I,  a  single  line  on  the  attached  printer 
(lives  N,  36(NiK),  and  T,  l.e.,  the  number  of  photos,  total  matrix  storage, 
and  number  of  photo-pair  linkages. 

I'- .  Error  Routine 

This  program  Is  designed  In  accordance  with  Standard  FMS  operating 
procedures.  Accordingly,  all  error  holts  hove  been  converted  to  a  transfer 
to  an  error  diagnostic  routine:  "ERROR",  The  error  routine  prints  the  con¬ 
tents  of  the  following  AC,  MQ,  SI,  KEYS,  XR|,  XR2,  XR4,  when  control 
reaches  the  TSX  ERROR,  4  instruction.  The  contents  of  the  AC  and  MQ  ore 
printed  both  In  octal  and  decimal.  The  complement  of  the  contents  of  XR)  ,  XR2, 
XR4  are  also  printed.  The  decrement  of  the  SI  contains  the  complement  of 
the  address  of  the  instiuctlon  which  colled  the  error  routine. 

following  this  is  a  dump  of  the  entire  core  which,  joined  with  the 
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knowledge  of  the  location  of  the  error  cell,  gives  maximum  offline  de¬ 
bugging  Information. 

The  relocation  constant  used  by  the  MIT-FMS  system  Is  144  octal.  Thus, 
If  for  example  the  line  printed  by  the  error  routine  gives  as  contents  of 
vhe  SI  07447700 . .  the  location  of  the  TSX  ERROR,  4  Instruction  Is 

BOl -1 44»3 135- 

An  error  stop  In  Segment  2  called  CEALR  prints  before  the  error 
inutlne  ERROR  Is  called  the  following  Information:  SEGMENT  2  ALARM 
Q  &  P  bits,  AC,  MQ,  XRI,  XR2 ,  XR4. 

The  content  of  XR4  Is  the  rolocatod  complement  of  the  location  con¬ 
taining  the  last  executed  TSX  Instruction.  Thus,  If  XR4C 72437,  the 
location  of  tho  TSX  Instruction  Is  !>34 1-144^1  75- 

A  stop  In  data  reading  could  mean  one  of  tho  following:  bad  Job  card, 
bad  parameter  card,  no  slop  card,  wrong  number  of  E  cards,  wrong  numbei 
of  ground  points-  A  stop  In  SEGMENT  2  ALARM  (CEALR)  Is  explained  f refl¬ 
uent  I  y  by  remarks  In  the  coding. 
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_IX.  RESULTS  AND  CONCLUSIONS 


The  Introduction  of  the  least  squares  method  into  the  analytical 
aerotr langu lat Ion  problem  for  the  adjustment  of  random  uiiurs  gives, 
according  to  the  theory  of  probability,  the  most  probable  solution.  This 
solution  consists  of  two  parts:  the  first  is  formed  by  the  location 
and  orientation  of  the  camera  stations  and  the  second  Is  formed  by  the 
adjusted  Imago  points  coordinates,  the  adjusted  ground  control  coordinates, 
and,  If  any,  the  ad  lusted  exposure  stations  control  Information. 

Thu  program  In  Its  present  form  computes  only  the  first  part  of  the 
solution.  Tho  reason  for  tills  is  that  to  compute  the  residuals  as  given 
by  the  least  squares  solution  requires  additional  computer  time  and 
storage  area  which  can  be  Justified  only  If  high  accuracy  Is  desired.  A 
discrepancy  of  say  10  microns  in  the  Image  coordinates  of  a  ground  point 
Is  equivalent  to  a  discrepancy  of  0.3  m.  or  I  ft.  on  the  ground  for  a 
nominal  scale  of  ~  “oogT  '  Is  *n  lh'b  laboratory  that  errors  of 

this  order  of  magnitude  can  be  tolerated  as  long  as  systematic  errors 
corrections  are  not  more  precise  than  they  arc  at  the  present  time.  The 
pho  togr  amnie  tr  Ic.  solution  given  by  the  program  is  then  based  on  the 
exposure  stations  location  and  orientation  as  given  by  the  least  squares* 
adjustment  and  on  the  unadjusted  Image  point  coordinates  and  control 
Informat  Ion. 

To  draw  a  general  conclusion  and  evaluation  of  the  program  requires 
that  a  general  testing  plan,  including  fict  tious  data  and  specifically 
selected  real  data,  be  carried  out.  The  tests  conducted  by  this  laboratory 
fall  short  of  a  conclusive  testing  pl^n,  expecially  as  concerning  the  real 
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data.  The  roason  for  this  Is  that  no  informal  ion  other  than  the  coordinates 
of  the  ground  control  points  was  available. 

Dufoiu  a  short  deScr i ipicn  and  analysis  of  the  tests  conducted  by  this 
laboratory  arc  presented,  the  units  of  tho  residuals  will  be  discussed. 

The  coefficients  of  the  residuals  In  the  condition  equations  are  scaled 
so  that  tho  residuals  of  the  Image  point  coordinates  are  in  microns.  The 
residuals  of  angles  and  for  ground  and  exposure  station  control  are  in 
units  of  0.0206  seconds  (scaling  by  lO^),  the  residual  of  h  Is  in  units 
of  feet,  the  weights  are  all  dimensionless  numbers,  ratios  of  a  reference 
variance  to  tho  respective  variances  of  the  quantities  to  be  weighted. 

A ,  FICTITIOUS  DATA 
1  ■  Tost  ‘i 


No. 

of  photos: 

6 

St  r  1 

1  ps : 

1 

No. 

of  ground  points: 

55 

0  . 

Method  of  solution! 

S 1 1  e  f  c  1 

No.  of  Iterations 

h 

No,  of  cycles  per  Iteration: 

2 

Weigh  is : 

a  1  1  one 

Pesu 1 ts 

tir. 

Pt ,  Type  A  |.  sec  . 

A'  sc-c. 

Ah  ft. 

101 

1)  -0.003 

-0.028 

503 

A  -0.003 

-0.003 

0.21 1 

307 

H 

0,388 

1 1 1 

A  -0.002 

0.001 

-0. 1  1 1 

5M 

F 

-0.001 

0.117 

where  &  ■  calculated  -  input 

b 

Method  of  solution: 

5tlefe! 

No.  of  1 tera 1 1 ons : 

4 

No.  of  cycles  per  Iteration: 

4 

Weights : 

a  1 1  one 
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Results 


Gr .  PT. 

AiL 

A  A. 

AL. 

101 

0 

0.000 

0.1 ''7 

503 

A 

0.000 

-0.001 

0.080 

307 

H 

0.318 

in 

A 

0.00! 

1 . 000 

0.075 

511 

F 

-0.001 

O.OS'i 

To 

transform  M  ani! 

1  A*  1  n  feet , 

a  rough 

est  (mate  Is 

see.  1 

ft.  for /v|.  and 

Too  scc-  c 

1  ft.  X 

COS  <|l  for  , 

2 .  Spec, | at  Test,, Data 


No.  of  photos: 

S  t  r i ps : 

No.  of  ground  points: 

6 

2 

1 

Mot hod  of  solution: 

Jordan 

diagonal i xat ion 

No.  of  Iterations: 

1 

Weights : 

a  1  1  one 

Result  s 

A  <)> 

JAA 

A.V  , 

Gt.  Pt. 

-0.001 

-0.001 

0.004 

Ex.  St .  1 

-0.001 

0.000 

0.002 

■>. 

-0.001 

0.000 

0.002 

3 

-0.001 

0,000 

0.001 

1, 

-0.001 

0.000 

0.003 

5 

-0.001 

0.000 

0.003 

6 

-0,001 

0.000 

0.002 

The  discrepancies  A|  A 

a!  Hn 

the  above  results  arc  not  beyond 

machine  truncation  errors  associated 

with  similar  errors  In  the 

Input  data,  nor  ore  they  beyond  the 

convergence 

limit  reached  In 

the 

last  iteration.  A  comparison 

of  the 

elements  of 

the  or ientat Ion 

matrices  of  the  exposure  stations  as 

compulcd  by 

the  program  arid 

as 

computed  bv  an  Independent  program  which  ge'ncrated  the  fictitious 
data  2  showed  an  agreement  of  4  digits  on  the  average.  This  corresponds 
to  a  dlscrepancey  in  the  angular  orientation  of  the  camera  station 
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-b 

of  the  order  of  magnitude  of  10  radians.  The  solution  for  the  Incre¬ 
mental  angular  displacements,  dT,  of  the  camera  stations  are  of  the  order 

-5 

of  magnitude  of  10  radians.  This  corresponds  to  a  discrepancy  on  the 

ground  of  0.2  ft.  (the  scale  of  data  2  Is  )  and  Is  In  agreement 

with  the  A(i>  ,  a\,  and  Ali  obtained.  Note,  however,  that  the  solution  Is 

computed  after  the  solution  for  the  dT's  Is  used  to  correct  the  camera 

orientations  so  that  the  final  discrepancy  can  be  of  the  order  of  0.2  ft. 

or  much  smaller  arid  In  fact  beyond  ilu:  capacity  of  the  computer.  The 

values  of  the  Incremental  displacements  dR's  are  of  the  order  of  magnitude 
~3 

of  10  ft.  and  clearly  do  not  affect  the  accuracy  of  the  solution.  We 
can  therefore  conclude  that,  for  the  fictitious,  data  tests  run  by  this 
laboratory,  the  results  are  satisfactory  and  the  small  discrepancies 
botwoen  Input  control  data  and  output  Is  In  the  range  of  the  truncation 
errors  and  convergence  limits. 

B .  REAL  DATA 

•  •  Test  10 

a.  No,  of  photos:  1 5 

No,  of  strips:  3 

No.  of  ground  pts:  64 

Method  of  solution:  Jordan  diagonal! zat ion 
No.  of  iterations:  4 
Weights :  all  one 

Resul ts  • 


Or.  Pi, 

Im 

A&. 

A* 

AL. 

9204 

A 

-0,025 

-0.004 

4.326 

9107 

A 

0.010 

-0,020 

-3.165 

9^03 

A 

0.055 

0.040 

-1 .638 

9407 

A 

0.01 1 

-0,039 

0.532 

9603 

A 

-0.018 

-0.019 

-2.133 

960b 

A 

-0.026 

0.058 

1 ,982 
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b.  No.  of  photos;  8 

No.  of  strips*.  2 

No.  of  ground  pts:  37 


Method  of  solution:  Jordan  diagonal izatlon 
No.  of  iterations:  4 


Weights : 

a  1 1  one 

ftesul ts: 

Gr.  Pt .  Type 

aJ> 

A  A 

aU 

9403  A 

0.034 

0.031 

0.074 

960?  A 

-0.034 

-0,003 

-0.020 

9606  A 

-0.012 

0.030 

0.080 

9407  A 

appears 

on  one  photo  only 

c.  No.  of  photos: 

8 

No.  of  strips: 

2 

No.  of  ground 

points: 

37 

Method  of  solution: 

St lefel 

No.  of  Iterations: 

6 

No.  of  cycle: 

per  Iteration:  3  and  2 

Rosy  Its:  for 

3  cycles  per 

1 torat Ion 

Gi  ,  Pt.  Typo 

Ml 

A  A 

Mi. 

9403  A 

0.024 

0.067 

5.806 

9603  A 

-0.005 

0.037 

-7.788 

960G  A 

0.020 

-0.004 

2.239 

Rosy  1  ts.:  for 

2  cycles  per 

1  terat Ion 

9403  A 

-0.01 1 

O.IJI 

(>.  1 1? 

9603  A 

0.023 

0.100 

-5.353 

9606  A 

0.060 

-0.053 

-0.997 

2.  Test,  19  HALCON 

No .  of  photos :  2h 

No-  of  s  t  r I ps :  3 

No.  of  ground  points:  140 


a,  Method  of  solution: 
No.  of  iterations: 
We i ghts ■ 


Jordan  d iagona I i zat Ion 

5 

a  1 1  one 
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fteiui  Its: 


Gr.  Pt, 

Ixaa 

ail-i 

M. 

302 

A 

0.005 

0.316 

28.535 

340 

A 

0.224 

-0.070 

-20.862 

223 

A 

-0.055 

0.060 

-  7.405 

104 

A 

-0,066 

-0.194 

-27-321 

152 

A 

-0. 109 

-0.155 

27.119 

all  points 

type  B 

0.00} 

0-310 

28.651 

0.223 

-0  072 

-20.463 

-0.056 

0.06) 

-  7-169 

-0,066 

-0.196 

-27-052 

-0. 1 10 

-0.153 

27.516 

b.  Method  of  solution:  Jordan  d i agonal i zul ion 

No .  of  Iterations:  6  and  no  convergence  weights 

[of  pass  points:  I 
of  control  points:  8 
Ground  control  coordinates:  ,>10000 
Input  "  output  of  2. a 


!kiiLUi: 


Gr.  Pt. 

Im 

/\<\> 

A,1 

Alt 

302 

A 

-0  009 

0.005 

0.286 

340 

A 

0.015 

-0.002 

-0.561 

223 

A 

0.013 

0.  006 

5  341 

104 

A 

-0.033 

-0-004 

0,346 

152 

A 

0.008 

-0.023 

1.279 

The  solution  for  dT's  and  dR's  shaved  that  after  the  second  Iteration 

they  alternate  In  sign  without  change  In  order  of  magnitude: 

max imum  dT -v  0. 01  radian 
maximum  dR  ~  1000  ft, 

The  results  of  I  Indicate  that  the  Stiefel  method  of  solution  compares 
very  poorly  with  the  Jordan  method.  Wo  will  disregard  it  in  the  following 
as  It  Introduces  significant  errors  foreign  to  the  photogrammetr Ic 
problem. 

The  results  in  I, a.  indicate  that  the  d i screpenc I es  M>  and  A\  when 
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transferred  Into  equivalent  feet  are  generally  of  the  same  order  of 
magnitude  asAkwith  the  maximum  being  approx  I  male! y  5  ft.  These 
results  must  be  considered  33  satisfactory  considering  that  the  systematic 
errors  have  not  been  corrected  and  that  Information  Is  not  available 
on  the  probable  accuracy  of  the  Input  data.  The  orders  of  magnitude 
cf  the  dT's  and  dR's  of  the  last  iteration  are  lO-5  radians  and  1  ft. 

respectively,  and  the  scale  of  the  model  is  20,000* 

The  results  of  l.b.  Indicates  smaller  A«.\> '  s,  A  A.'s,  and  All's  for 
the  same  ground  control  as  l.o.  They  show  further  that  then's  are 
significantly  smaller  than  the  A'ii's  and/SA's  In  equivalent  feet.  This 
merely  reflects  the  facts  that  the  data  In  l.b.  Is  only  slightly  re¬ 
dundant  os  regard  to  altitude  Information.  The  lower  degree  of  control 
redundancy  l«  A.b.  as  compared  to  A. a.  would  also  be  the  reason  for 
smol lor  A^'s  and  AA's  in  l.b. 

In  2. a.  tho  discrepancies  reach  30  ft.  Vihon  one  considers  that 
with  strong  control,  such  os  in  2.b.,  the  A i's.&AU,  andA^'s  become 
very  small  In  comparison  with  the  results  of  2. a.,  even  without  converg¬ 
ence,  it  becomes  apparent  that  the  data  contains  significant  errors. 

A. a.  shows  that  with  a  uniform  weighting  It  Is  impossible  to  lit  the 
ground  control  without  a  "stretching"  In  the  horizontal  direction 
reaching  30  ft.,  and  l.b.  shows  that  If  one  enforces  strongly  the  ground 
control  along  with  the  corresponding  image  coordinates,  the  adjustment, 
of  the  remaining  Image  coordinates  In  each  ileiotion  is  so  high  as  to 
prevent  convergence.  The  errors  In  the  data  9"s  ot  two  kinds,  ground 
control  errors  and  Image  coordinates  errors.  As  was  pointed  out  previously, 
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this  laboratory  has  no  knowledge  of  their  order  of  magnitude  or  of 
the  occurrence  of  these  errors.  Something  can  be  said,  however,  about 
atmospheric  refraction  errors.  For  an  altitude  of  30000  ft,,  an  angle 
with  the  vertical  of  45"  (which  takes  account  of  tilt)  and  a  focal 
length  of  1  ft.,  the  atmospheric  refraction  error  on  a  vertical  photo¬ 
graph  Is  about  50  microns  which  In  turn  produces  an  apparent  ground 
displacement  of  about  5  ft.  Other  systematic  errors,  such  as  film 
shrinkage  and  lens  distortion,  Increase  the  probable  error.  This  accounts 
for  at  least  part  of  the  values  of  the  A  •  h  ,  A  a  1  s  ,  andAW's.  The 
remaining  port  would  be  caused  by  errors  in  the  ground  control  coordinates. 
As  the  least  squares  method  Is  not  Intended  to  adjust  for  large  systematic 
errors,  tho  results  of  2. a.  are  not  conclusive.  The  results  of  2.b., 
however,  show  that  the  weighting  does  Influence  significantly  the 
solution  and  con  oven  prevent  convergence  when  It  Is  not  In  accordance 
with  the  actual  accuracy  of  the  data.  The  A(v's  are  seen  to  be  In  general 
of  the  same  order  of  magnitude  os  theAgt's  andAA's  expressed  In  feet. 

This  would  likely  hold  also  for  pass  points  in  which  case  the  general 
method  of  least  squares,  Independently  of  an  accurate  weighting,  would 
present  an  Impiovement  on  provtous  photogrammotr Ic  solutions. 
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X.  RECOMMENDATIONS 


Tho  possibilities  and  flexibility  of  the  least  squares  method  should 
be  studied  further.  In  order  to  do  this,  it  Is  recommended  that  a 
routine  to  compute  the  residuals  or  adjustment  to  tho  Input  data  and 
the  sum  of  the  weighted  squares  of  those  residuals  be  included  in  tho 
program.  For  data  corrected  for  systematic  errors  and  whose  random  errors 
do  noton  tho  average  exceed  a  certain  limit,  the  additional  computer  time 
and  storage  area  required  may  overshadow  the  resulting  gain  in  probable 
accuracy. 

However ,  for  gaining  o  knowledge  about  the  order  of  magnitude  of  the 
errors  contained  in  a  set  of  data  when  not  available  otherwise  as  well 
as  for  obtaining  a  more  accurate  solution  when  important  errors  cannot  be 
corrected  beforehand,  the  calculation  of  the  residuals  would  be  the  only 
ol tornat Ivo. 

Consider  the  following.  A  sot  of  data  with  minimum  ground  control 
but  redundant  as  regards  tho  relative  orientation  of  tho  camera  stations 
Is  run.  At  eoch  Iteration  the  adjustment  of  the  data  as  given  by  the 
least  squares  method  required  to  make  the  rays  Intersect  as  well  as  the 
sum  of  the  Sv|uaies  of  the  residuals  and  thulr  standard  deviation  can  be 
computed.  An  Iteration  Is  reached  where  the  above  sum  and  the  correspond¬ 
ing  standard  deviation  arc  minimum. 

For  truly  unrelated  random  errors,  the  standard  deviation  computed 
as  outlined  above  Is  a  true  standard  deviation  In  the  probability  theory 
sense.  For  errors  of  any  kind  throughout  the  entire  data  Set,  an  estimate 
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is  p  -vided  of  the  order  of  magnitude  of  the  errors  effecting  the  image 
points  coordinates.  A  concentration  of  errors  in  pans  of  the  data  can  be 
discovered  by  isolating  different  combinations  of  m  ■ndant  subsets  of 
the  data  and  comparing  the  corresponding  standard  deviations.  When  minimum 
ground  control  is  used,  no  discrepancy,  beyond  truncation  errors  and 
convergence  limits,  between  input  and  output  values  of  the  ground  control 
coordinates  will  occur,  but  when  obtaining  the  groum  coordinates  of  the 
redundant  ground  control  points  which  were  used  only  as  pass  points, 
discrepancies  will  occur  If  they  cannot  be  accounted  for  by  a  reasonable 
error  in  the  images  coordinates,  they  provide  an  estimate  of  the  accuracy 
of  the  ground  control  information. 

Although  not  recommended  as  sound  practice  but  still  useful  when 
no  other  course  is  possible,  uncorrected  data  can  be  adjusted  using  the 
computed  residuals  given  by  the  last  iteration  before  the  solution  is 
computed.  Finally,  when  systematic  errors  are  removed  from  the  data  and 
information  on  the  standard  deviations  of  the  random  er  ors  is  available, 
the  adjusted  input  date  can  be  used  to  obtain  high  acco  acy  results. 
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